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Introduction 



Motivations et plan de la these 

Si I'astronomie dans le domaine visible est millenaire, I'astronomie 7 n'est que cinquantenaire et 
n'en etait qu'a ses balbutiements il y a peu de temps encore. L' atmosphere terrestre nous protege 
en grande partie de ces photons tres energetiques, ainsi 1' Homme ne peut observer le ciel dans ce 
domaine d'energie que depuis le vingtieme siecle, en plagant ses instruments sur des ballons ou des 
satellites en orbite autour de la Terre, ou en ayant recours a des telescopes a effet Tcherenkov au 
sol pour ce qui est des photons de tres haute energie (E > 30 GeV). L'astronomie 7 renseigne sur 
des phenomenes violents, capables d'accelerer des particules chargees vers des energies extremes : 
noyaux actifs de galaxie, vestiges de supernova, sursauts 7, pulsars, etc. Le Large Area Telescope 
(LAT), a bord du satellite Fermi, appartient a la troisieme generation de telescopes 7 places en orbite. 
Le satellite Fermi a ete lance le 1 1 juin 2008, et le LAT a debute ses observations entre 20 MeV et 
plus de 300 GeV quelques jours apres. Sa surface effective de collection de photons 7 est considera- 
blement meilleure que celle de son predecesseur, EGRET (actif de 1991 a 2000). Le champ de vue du 
LAT est plus large, et sa resolution angulaire plus fine. Mais le LAT ne se resume pas a une evolution 
d'EGRET : il observe le ciel 7 entre 10 et 100 GeV, domaine du spectre electromagnetique largement 
inexplore jusqu'alors. Ainsi, le LAT observe desormais plusieurs centaines de sources, connues ou 
inconnues auparavant, avec un niveau de sensibilite sans precedent. De plus, le LAT balaye le ciel 
presque uniformement toutes les trois heures, a la difference d'EGRET qui operait principalement en 
pointe, ce qui fait du LAT un instrument de choix pour la decouverte de phenomenes transitoires. 

L'astronomie des pulsars est plus jeune encore que I'astronomie 7, le premier pulsar ayant ete 
decouvert en 1967 par J. Bell & A. Hewish a Cambridge. On connait a present 1900 de ces objets 
environ. Les pulsars sont des etoiles a neutrons fortement magnetisees et en rotation rapide, emettant 
un rayonnement sous forme de faisceaux qui balayent le ciel au cours de leur rotation. Pour une vaste 
majorite, les pulsars sont observes dans le domaine radio. Leur etude a ces longueurs d'onde touche 
a une large gamme de domaines de la physique, de la relativite generate a la mesure du champ ma- 
gnetique dans la Galaxie ou la recherche d'ondes gravitationnelles. D'autres pulsars sont egalement 
detectes en optique, en rayons X ou en rayons 7. Cette categoric de pulsars demeure tout de meme 
largement minoritaire : seuls neuf pulsars avaient par exemple ete detectes dans le domaine 7 au de- 
but de I'activite de Fermi, grace a SAS-2, COS-B, COMPTEL, EGRET et AGILE. Ce faible nombre 
de detections a laisse de nombreuses interrogations, notamment en ce qui conceme le mecanisme de 
production de rayonnement 7 dans leur magnetosphere. En outre, ces objets ont en commun d'etre 
relativement jeunes : aucun pulsar milliseconde n'a ete detecte avec certitude par EGRET ou AGILE. 
Des pulsations faiblement significatives ont neanmoins ete observees pour deux de ces objets, PSR 
J0218+4232 et J1824— 2452A. La question de remission dans le domaine 7 par les pulsars milli- 
seconde, etoiles a neutrons en rotation tres rapide, representant 10% des pulsars connus et pouvant 
disposer de reservoirs d'energie rotationnelle colossaux (typiquement, 10^^ a 10^^ erg, ce qui est com- 
parable a certains sursauts 7), restait done en suspens. Les observations du ciel 7 par le LAT offraient 
la possibilite d'eclaircir certaines zones d'ombre : les pulsars milliseconde peuvent-ils accelerer des 
particules chargees vers les tres hautes energies, afin que celles-ci produisent un rayonnement 7 detec- 
table a distance ? Si oui, quelles sont les particularites de ce rayonnement 7, par rapport a remission 
des pulsars milliseconde dans les autres longueurs d'onde d'une part, et par rapport a celui des pulsars 
ordinaires ? 

Cette these, consacree a la recherche de pulsars milliseconde en rayons 7 avec le LAT, est a la 
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jonction de ces deux astronomies. L' observation des pulsars a haute energie renseigne sur les me- 
canismes operant dans la magnetosphere de ces objets aux proprietes singulieres, d'autant plus que 
pour les pulsars detectes en rayons 7, on s'aper9oit qu'une bien plus grande part de I'energie rotation- 
nelle est convertie en rayonnement yplutot que dans le domaine radio. Plus generalement, I'etude de 
remission /des pulsars se place dans le contexte de revolution des etoiles massives. Par ailleurs, les 
premiers resultats du LAT indiquent qu'en majorite, les sources galactiques sont des pulsars, rendant 
ces objets incontournables pour toute etude de sources galactiques. II a egalement ete propose que les 
pulsars proches soient les principaux responsables de remission des electrons et positrons detectes 
autour de la Terre. Ces differentes raisons motivent I'etude des pulsars a haute energie. 

La problematique peut se resumer de fa9on relativement simple. Aucun pulsar milliseconde n'a 
ete fermement detecte dans les donnees d'EGRET et d' AGILE, ce qui laisse a penser que si ces objets 
emettent en rayons 7, ils doivent etre peu lumineux. Par consequent, les photons eventuellement emis 
sont espaces dans le temps, de sorte que la construction de leur courbe de lumiere (c'est-a-dire I'in- 
tensite du signal en fonction de la fraction de tour du pulsar sur lui-meme) requiert une connaissance 
tres precise de leur periode de rotation et de leur orbite, pour les objets appartenant a des systemes 
binaires. Ces differents parametres sont connus via la chronometrie, en particulier dans le domaine 
radio. La decouverte eventuelle de pulsations en rayons 7 pour ces pulsars a rotation tres rapide re- 
pose sur leur chronometrie, en plus des performances du telescope. La datation des pulsars a done ete 
un premier objectif de cette these. Une fois les parametres de chronometrie acquis, la recherche de 
pulsations dans les donnees du LAT a ete la suite logique du travail preparatoire. 

Dans la premiere partie de ce manuscrit nous rappelons un certain nombre de proprietes des 
pulsars, notamment au travers d'une description simplifiee consistant a les assimiler a des dipoles 
magnetiques toumants. Nous introduisons ensuite les pulsars milliseconde, objets etudies dans cette 
these. Leur scenario de formation et certaines de leurs proprietes se distinguent de ceux des pulsars 
communs. Dans un second volet de cette partie introductive nous etablissons un etat des lieux des 
connaissances sur les pulsars dans le domaine 7, du point de vue theorique et observationnel. 

La seconde partie est consacree aux aspects instrumentaux concemant la preparation et 1' analyse 
des donnees du LAT. Dans un premier temps, nous introduisons le concept de chronometrie des pul- 
sars, puis nous decrivons la campagne de chronometrie des pulsars d'interet pour le LAT, dont la 
coordination a ete I'un des objectifs au cours de cette these. Dans un second temps, nous rappelons 
le principe de fonctionnement et les caracteristiques du Large Area Telescope et de ses composants 
principaux. 

Dans la troisieme et demiere partie, nous presentons la recherche de pulsations pour des pulsars 
milliseconde dans les donnees enregistrees par le LAT. Comme nous le verrons, huit de ces objets ont 
ete detectes, les etablissant pour la premiere fois en tant que puissants accelerateurs de particules char- 
gees. Les courbes de lumiere de ces huit pulsars sont presentees, ainsi que les resultats spectraux. Les 
differentes caracteristiques sont interpretees en terme de modele d' emission dans la magnetosphere. 
Quelques pulsars milliseconde semblent etre de bons candidats a la detection par le LAT au cours des 
mois a venir. 



14 



Introduction 



Publications 

Ce manuscrit porte essentiellement sur la decouverte de pulsations pour huit pulsars milliseconde 
en rayons 7. Ce resultat est du aux bonnes performances du LAT, mais aussi au succes de la campagne 
de suivi des pulsars dans les domaines radio et X, a laquelle j'ai activement contribue. Grace a I'expe- 
rience concemant le probleme des calculs de phases rotationnelles acquise au travers de la campagne 
de chronometrie, j'ai egalement contribue a 1' analyse de pulsars jeunes observes par le LAT, Vela et 
PSR J202 1+3651. Plus generalement, I'ensemble des analyses de pulsars emetteurs radio et X avec le 
LAT ont jusqu'a present beneficie de la chronometrie des pulsars. Mes differents travaux ont trouve 
leur concretisation dans les articles suivants : 

- J. T. O'Brien, S. Johnston, M. Kramer, . . ., L. Guillemot, PSR J1410—6132 : a young, energe- 
tic pulsar associated with the EGRET source 3EG J1410-6147, MNRAS 388, 1 (2008) 

- D. A. Smith, L. Guillemot, F. Camilo, et ah. Pulsar timing for the Fermi gamma-ray space 
telescope, A&A 492, 923 (2008) 

- A. A. Abdo, et al. (Fermi-hAJ' Collaboration), Pulsed gamma-rays from the millisecond pulsar 
J0030+0451 with the Fermi Large Area Telescope, ApJ 699, 1171 (2009) 
Corresponding authors : L. Guillemot, T. J. Johnson, M. Kerr 

- A. A. Abdo, et al. (Fermi-hAT Collaboration), A population of gamma-ray millisecond pulsars 
seen with the Fermi Large Area Telescope, Science 325, 848 (2009) 

Corresponding authors : L. Guillemot, T. J. Johnson, M. Kerr, D. A. Smith 

- A. A. Abdo, et al. (Fermi-LAI Collaboration), Pulsed Gamma-rays from PSR J2021+3651 
with the Fermi Large Area Telescope, ApJ 700, 1059 (2009) 

Corresponding authors : D. A. Smith, L. Guillemot, M. Kerr 

- A. A. Abdo, et al. (Fermi-LAI Collaboration), Fermi Large Area Telescope Observations of the 
Vela Pulsar, ApJ 696, 1084 (2009) 
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Resume 

LES pulsars sont des « etoiles a neutrons hautement magnetisees et en rotation rapide ». Dans 
ce premier chapitre nous rappelons les grandes etapes de la decouverte des pulsars. On peut 
les assimiler en premiere approximation a des dipoles magnetiques tournants, entoures d'une 
magnetosphere en co-rotation qui est le siege de remission de rayonnement a differentes energies. 
On observe que la rotation des pulsars sur eux-memes est pseudo-periodique, la rotation ralentissant 
progressivement. L'hypothese du dipole magnetique tournant permet d'estimer la quantite d'energie 
cinetique rotationnelle perdue au cours du temps par les pulsars, ainsi que leur age ou I'intensite 
de leur champ magnetique, a partir de la periode de rotation P et de sa diminution P. On s'aper§oit 
de plus que differents types de pulsars existent, en fonction de P et de P. Dans la deuxieme moitie 
de ce chapitre nous nous concentrons sur une des categories, les pulsars milliseconde, et rappelons 
leur specificite. Ceux-ci ont des periodes de rotation comprises entre 1 et 30 ms, avec des taux de 
ralentissement extremement faibles. lis appartiennent generalement a des systemes binaires, une 
consequence de leur formation et leur evolution. Bien que representant environ 10% des pulsars 
connus, ils ont une grande importance en astronomic et en physique en general, de par leur nature 
extreme. 
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CHAPITRE I. INTRODUCTION AUX PULSARS 



I.l Historique 



Deux ans seulement apres la decouverte du neutron par James Chadwick, en 1932, Walter Baade 
et Fritz Zwicky proposaient que des etoiles a neutrons, c'est-a-dire des etoiles composees presque 
exclusivement de neutrons, pourraient naitre au coeur des supemovae. Et quelques annees plus tard, 
en 1939, Oppenheimer et Volkoff entreprirent meme les premiers modeles de structure des etoiles a 
neutrons (Oppenheimer & Volkoff 1939j ). Malgre ces travaux pionniers, la communaute astronomique 
n'a guere accorde d'interet a la recherche de ces astres ; on n'avait pas entrevu alors que les etoiles a 
neutrons pourraient etre, par exemple, detectees via la recherche de fluctuations rapides dans le signal 
radio. Si bien qu'il a fallu attendre plus de 25 ans avant que Ton ne decouvre une etoile a neutrons... 
par hasard. 



En 1967, Jocelyn Bell et son directeur de these Anthony Hewish scrutaient le ciel avec le radio- 
telescope de Cambridge. Une instrumentation speciale avait ete con^ue dans le but de detecter des 
scintillations breves du signal radio, une signature caracteristique des quasars. Le 6 aout 1967, des 
trains d'ondes radio periodiques furent remarques a 19'' 19*" d'ascension droite et +21° de declinaison. 
Le signal fut confirme les jours suivants, au transit de cette region du ciel au-dessus du telescope. La 
structure des bouffees d'emission radio restait inaccessible cependant, la resolution en temps etant 
trop grossiere. Une nouvelle instrumentation fut alors mise au point pour sonder le signal en de9a de 
la seconde. Le 28 novembre 1967, on s'aper§ut que les trains d'ondes etaient composes de pulsations 
espacees de 1,337 s environ. Cet astre, source de si gnal radio pulse, a ete nommee pulsar. Ce pulsar 
est maintenant connu sous le nom de B 1 9 19+2 Ipj pewish et al. \ 1968| ). 



Le hen entre les pulsars et les etoiles a neutrons avait ete fait quelques temps avant la decouverte 
de B1919+21. Pacini avait en effet postule que la source d'energie alimentant la nebuleuse du Crabe 



etait une etoile a neutrons fortement magnetisee et en rotation rapide (Pacini 1967 ). Independamment, 
Gold introduisait le concept de « pulsars alimentes par leur rotation » : a savoir que des etoiles a neu- 
trons en rotation etaient analogues a des dipoles magnetiques tournant, done perdant de I'energie par 
rayonnement electromagnetique et par emission de particules relativistes ( Gold||1968 1. Ce faisant, la 
reserve d'energie (rotationnelle en majorite) diminuait, se traduisant par une diminution progressive 
de la periode de rotation. Mais il demeurait la possibilite que le signal radio, de periode 1,337 s, pro- 
vienne en realite de I'oscillation resonnante de naines blanches ( |Meltzer & Thome||1966 ^, ce modele 
expliquant naturellement des periodicites de I'ordre de la seconde. Les decouvertes successives en 
1968 d'un pulsar de periode 89 ms au coeur du reste de supernova de Vela ( Large et a/.|1968 ), et d'un 
pulsar de periode 33 ms au coeur de la nebuleuse du Crabe ( [Staelin & Reifenstein|1968| ) etablissaient 
I'identite des pulsars en tant qu'etoiles a neutrons de fagon definitive. 



La convention de nommage des pulsars est la suivante : les pulsars decouverts de 1967 jusqu'au milieu des annees 90 
sont nommes par leur position dans le ciel dans le systeme de coordonnees equatoriales besselien (B1950) precedee d'un 
B : Bxxxx+yy. Ceux decouverts par la suite sont nommes suivant leur position dans le systeme julien (J2000) precedee 
d'un J : Jxxxx+yyyy. Lorsque deux objets ont le meme nom suivant cette convention, ce qui est le cas pour certains pulsars 
des amas globulaires, on les distingue en ajoutant une lettre : par exemple J1824— 2452A. 
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1.2. PULSARS : PROPRIETES GENERALES 



1.2 Pulsars : proprietes generales 

1.2.1 Etoiles a neutrons 

Les etoiles a neutrons sont les residus compacts des explosions de supemovae de type II, engen- 
drees par I'effondrement gravitationnel du coeur des etoiles massives (typiquement > 8 M©) en fin de 
vie, lorsque la pression de radiation ne suffit plus a compenser leur contraction. Tandis qu'un objet tres 
dense se forme au coeur de I'etoile progenitrice, les couches extemes de 1' atmosphere stellaire sont 
ejectees vers le milieu ambiant, et forment par la suite le vestige de supernova {supernova remnant, 
SNR). 

Les modeles de structure des etoiles a neutrons placent typiquement leur masse entre 0,5 et 



2,5 M© (pour une revue detaillee, voir Lattimer & Prakash (2004)). Un certain nombre de masses de 



pulsars ont ete mesurees {cf. Figure 3 de Lattimer & Prakash (2007)), les mesures les plus precises 



provenant de la chronometrie de pulsars binaires. Les masses relevees sont generalement compatibles 
avec 1,4 M©, valeur de masse la plus couramment utilisee pour les etoiles a neutrons. En revanche, 
on ne dispose a ce jour d'aucune mesure precise de rayon d'etoile a neutrons. Les modeles d'equation 
d'etat de la matiere, qui decrivent la relation entre densite et pression au sein de I'etoile, placent 



generalement ce rayon entre 10 et 12 km, comme le montre la Figure I.l (A). La mesure simultanee 
du rayon et de la masse d'une etoile a neutrons placerait des contraintes fortes sur les differents 
modeles. 

L' equation d'etat de la matiere definit le profil de densite radial de I'etoile, et done le moment 
d'inertie / = kMR^. Pour une sphere de densite homogene, onak= |, alors que les modeles predisent 
des valeurs variant entre ,3 et 0,45 lorsque le rapport M/R est compris entre 0,1 et 0,2 M© km~^ 
([Lattimer & Prakash|26oi]). En utilisant A: = 0,4, M = l,4M0eti?= 10km, onobtient/~ lO'^^gcm^, 



valeur que Ton utilisera dans la suite. Pour fixer les idees, la Figure|rT](B) montre quelques exemples 
de relations masse-rayon ou Ton utilise k = 0,4 et Ton fait varier / entre 0,5 Iq et 1 ,5 /q, avec Iq = 10^^ 
g cm^. Bien que les differents modeles d'equations d'etat de la matiere vont au-dela de I'hypothese 



de la sphere de densite uniforme (cf. Figure 6 de Lattimer & Prakash (2001 )), on s'aper9oit que la 



valeur du moment d'inertie de 10^^ g cm^ est tres incertaine, ce que nous garderons a 1' esprit dans 
la suite de ce manuscrit. File pourrait varier fortement d'un pulsar a un autre, en particulier s'il y a 
accretion de matiere au cours de la vie de I'etoile a neutrons, comme c'est a priori le cas pour les 
pulsars milliseconde, comme nous le verrons par la suite. 

1.2.2 Description simplifiee : modele du dipole 

Un certain nombre de proprietes essentielles des pulsars peuvent etre comprises en les assimilant 
a des dipoles magnetiques tournants (voir Figure [L2]). 

Le champ magnetique dipolaire, dont I'axe n'est pas necessairement aligne avec I'axe de rotation, 
induit en toumant un champ electrique E (^Q.xr^xB,ouQ. est la vitesse angulaire et B est le 

champ magnetique regnant a la distance r du dipole. Dans ces conditions, Goldreich & Julian ont 
montre que I'etoile a neutrons ne peut etre entouree de vide : en effet la composante de champ elec- 
trique parallele au champ magnetique a la surface de I'etoile est intense, et arrache des particules 
chargees vers I'exterieur, peuplant ainsi la magnetosphere ( Goldreich & Julian||l969| ). A I'equilibre, 



la magnetosphere est remplie d'une distribution de charges p = —Cl.B/ {lite), et la composante 
du champ electrique est ecrantee. Les particules chargees et les champs electrique et magnetique de 
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(a) Modeles d'equation d'etat de la matiere des etoiles a neutrons 
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Figure I.l - (A) : diagramme masse-rayon pour certains modeles d'equation d'etat de la matiere des 
etoiles a neutrons. Voir |Lattimer & Praka sh (2001 ) et |Lattimer & Prakash| ( 2007] l pour la nomenclature des 
differents modeles representes et des zones d'exclusion. Figure issue de Lattimer & Prakash] ( |2007| l. (B) : 
courbes representant / — ^MR^, pour differentes valeurs du moment d'inertie / autour de Iq = 10*-' g cm^. 
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la magneto sphere entrent en co-rotation avec I'etoile. Un pulsar de periode de rotation P est entoure 
d'un « cylindre de lumiere » imaginaire, de rayon cP/{2tz), qui est en co-rotation avec le pulsar, 



et dont la surface se deplace a la vitesse de lumiere. Comme le montre la Figure 1.2 les lignes de 
champ magnetique comprises dans ce cylindre de lumiere peuvent se refermer normalement (ce sont 
les lignes de champ « fermees »), tandis que celles s'etendant au-dela du cylindre ne peuvent boucler 
(lignes de champ « ouvertes »). Une repartition de charges existe done la oia les lignes de champ 
peuvent se refermer. Au contraire, les lignes de champs ouvertes laissent s'echapper les particules 
chargees au-dessus des poles magnetiques. Par cascades electromagnetiques, les particules chargees 
engendrent remission de photons dont les directions sont delimitees par les demieres lignes de champ 
ouvertes, celles qui interceptent le cylindre de lumiere. Ainsi un faisceau d'emission centre autour de 
I'axe magnetique du pulsar se forme. II balaye I'environnement du pulsar au cours de la rotation. Si le 
faisceau intercepte la trajectoire de la Terre, I'observateur terrestre peut detecter un signal periodique, 
se repetant a chaque rotation du pulsar ( jLorimer & Kramer [2004 1. 



Ce modele donne une interpretation simple de I'origine de remission des pulsars dans le domaine 
radio, comme creee par les particules chargees se depla9ant le long des lignes de champs ouvertes. 
Comme nous le verrons par la suite, remission de haute energie des pulsars n'est pas necessairement 
produite au-dessus des poles. Neanmoins I'idee d'une zone d'emission situee dans la magnetosphere, 
visible de fa§on periodique par I'observateur terrestre si la configuration geometrique le permet, est 
toujours applicable. 
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1.2.3 Periode, ralentissement et quantites derivees 

La description faite ci-dessus de la magneto sphere des pulsars, bien que tres incomplete, indique 
qu'il existe au moins deux sources de perte d'energie pour le pulsar et sa magnetosphere : remission 
de rayonnement dipolaire a basse frequence, due a la rotation du moment magnetique, et I'echappe- 
ment de particules chargees le long des lignes de champ ouvertes, aboutissant a un faisceau d'emission 
le long de I'axe magnetique. L'energie cinetique rotationnelle E = jIO?, ou I est le moment d'inertie 
de I'etoile a neutrons et Q. = In/P la vitesse angulaire, est le principal reservoir d'energie du pulsar. 
La diminution de E avec un taux E = entraine par consequent le ralentissement de la rotation 

du pulsar, avec un taux P = —InCl/Q? > 0. La quantite E est la perte d'energie due au freinage 
(Spin-down power). En fonction des parametres P et P, et en prenant / = 10^^ g cm^, on a : 



^ = 4^'^^-3,95xl03iergs-i(^)(^) ' (LI) 

Ces pulsars dont la periode augmente au cours du temps en raison de leur perte d'energie rota- 
tionnelle forment la categorie des pulsars « alimentes par rotation » (Rotation-Powered). On connait 
a ce jour environ 1900 de ces astres, parmi lesquels B 1919+21, Vela et le Crabe, cites auparavant . 



Le catalogue ATNlj^repertorie ces pulsars, et certaines de leurs proprietes ( Manchester et a/. [2005 ) 
C'est cette categorie de pulsars qui est etudiee dans cette these. Une autre categorie d'objet est celle 
des pulsars « alimentes par accretion » (Accretion-Powered). Le rayonnement de ces demiers est du 
a la chute de matiere provenant d'une etoile compagnon. La forte temperature acquise par la matiere 
dans sa chute cree des points chauds a la surface de I'etoile a neutrons. Ces zones sont emettrices de 
rayons X et sont en rotation avec le pulsar, ce qui cree un signal pulse pour un observateur fixe. On 
comprend cependant que le mecanisme de rayonnement de ces objets est tres different de celui des 
pulsars alimentes par leur rotation. 

Si Ton suppose que la perte d'energie cinetique rotationnelle du pulsar est due exclusivement au 
rayonnement dipolaire magnetique, on a (Jackson|1962) : 



E = -laa 



3c3 



(1.2) 



Dans cette equation, Jx designe le moment dipolaire magnetique du pulsar, et a est Tangle que 



forment I'axe magnetique et I'axe de rotation. D'apres L2 il vient : 



Cette demiere equation est le cas particulier de Q pour une perte d'energie purement di- 
polaire, soit n = 3 d'apres [L3| La quantite n est appelee « indice de freinage ». Cet indice peut se 
mesurer : en effet Cloc Q," implique n = et comme nous le verrons dans la suite, Clet Cl 

sont accessibles via la chronometrie des pulsars. En pratique les mesures de n sont rares, car Q. est 



souvent contamine par des instabilites de la rotation (le timing noise, cf. IIL2.3 1, dont il est necessaire 



Disponible en ligne : http ://www.atnf.csiro.au/research/pulsar/psrcat/expert.html 
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de s'affranchir pour acceder au comportement intrinseque du pulsar. Les mesures disponibles sont 
en contradiction avec I'hypothese n ~ 3, par exemple le pulsar du Crabe : n = 2,51 ±0,01 ( Lyne| 
et al. 1993), ou Vela : n = 1,4±0, 2 ( |Lyne et ar]|1996[ ). II est generalement admis que I'ecart entre 



les indices de freinage mesures et la valeur theorique de 3 est une preuve qu'une part importante de 
I'energie cinetique est dissipee via un vent de particules chargees et de champ magnetique emis par 
les pulsars. 

La relation de proportionnalite Q.oc Q/^ peut s'integrer simplement, a condition que n ^ \, pour 
estimer I'age du pulsar : 



P 

T = r 

{n-\)P 

Dans cette equation, T designe I'age du pulsar, P et P sont les valeurs actuelles de la periode et 
du ralentissement, et Pq est la periode du pulsar a sa naissance. En faisant I'hypothese que la periode 
initiale Pq est negligeable devant la periode actuelle, et que le ralentissement est d'origine dipolaire 
uniquement (soit n = 3), on a : 



La quantite T est nommee « age caracteristique », et represente une estimation raisonnable de I'age 
veritable du pulsar. Le pulsar du Crabe a par exemple un age caracteristique de 1240 ans, alors que 
r association avec la supernova historique indique un age de 950 ans, soit un ecart de 25% environ. 
Pour la quasi-totalite des pulsars, I'age historique est inconnu, ainsi I'age caracteristique T est le seul 
estimateur dont on dispose. 

Enfin, les valeurs de P et P et le modele dipolaire de champ magnetique des pulsars permettent 
d' estimer le champ magnetique regnant en differents endroits de la magnetosphere, et en particulier a 
la surface du pulsar et au cylindre de lumiere. On peut montrer que I'intensite du champ magnetique 
induit par le moment dipolaire a la distance r est donne par : 



Bir) 



' 3/c3 




2 sin^ a 





-PP (1.6) 



Soit Bs le champ magnetique en surface, dit « champ magnetique caracteristique », et B^c le 
champ magnetique au cylindre de lumiere, de rayon Ri = cP/{2Tt). Pour un pulsar canonique, de 
rayon R = 10 km, de moment d'inertie / = 10^^ g cm^, et en choisissant un angle d'inclinaison 
magnetique a egal a 90° (rotateur orthogonal), on a : 

U9 ■ 



Bs = B{R) ~ 3,2 X lO^^G VPP (1.7) 
flTtR^ ^ 



Blc = Bs\^-^] (1.8) 
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II faut souligner que ces deux dernieres equations utilisent une description dipolaire du champ, 
alors que les mesures d' indices de freinage disponibles indiquent un mecanisme de perte d'energie 
rotationnelle plus complexe. Les incertitudes existant sur R et I sont importantes, et enfin une valeur 
de 90° pour I'inclinaison magnetique a ete utilisee, alors que certaines valeurs de a mesurees sont 
tres differentes : citons par exemple le cas de PSR J02 18+4232, pour lequel cet angle est proche de 



0° (Stairs et a/.||1999). II en va de meme pour E et I'age caracteristique T : les parametres sont mal 



connus, ainsi les incertitudes peuvent etre du meme ordre de grandeur que les valeurs elles-memes. 



Mais a defaut d'informations plus precises, les equations |I . 1 [ p3] et pT8] fournis sent des renseignements 
utiles sur les pulsars et leur magneto sphere. 

Remarquons que des calculs alternatifs des proprietes energetiques et magnetiques des pulsars 
existent. Par exemple, Spitkovsky a applique les equations de la MHD a la magnetosphere des etoiles 
a neutrons dans I'hypothese d'un plasma libre et obtenu une solution numerique a ces equations 
( |Spitkovsky|2006[ ). La perte d'energie par rayonnement electromagnetique ainsi calculee est : 



£'spitkovsky = ki (1+^2 sin^ a) (1.9) 
Les coefficients k\ et k2 valent respectivement (1 ± 0,05) et (1 ± 0,1). Le rapport entre ^spitkovsky 



et la perte d'energie cinetique rotationnelle par rayonnement dipolaire {cf. 1.2 1 est done : 



£'spitkovsky 3 1+^2 sin^ a 
^'oipoie 2 sin^ a 

A r inverse, le champ magnetique surfacique calcule numeriquement est proportionnel a 

(l +fc2sin^ a) tandis que le champ magnetique de surface dans I'hypothese du dipole tournant 

est proportionnel a (sin^ a) . Les formules de Spitkovsky fournissent done une description plus 
reahste des proprietes des pulsars, en particulier lorsque a est faible. En effet : pour un angle a 
proche de 0, £'Dip6ie tend vers et le champ magnetique surfacique correspondant diverge tandis que 
les quantites calculees numeriquement prennent des valeurs finies. Neanmoins, Tangle d'inclinaison 
magnetique a est inconnu pour la plupart des pulsars, aussi nous utiliserons les equations [LT| [13] et 



1.8 avec a = 90° dans la suite de ce manuscrit. Notons tout de meme qu'a mesure que des valeurs de 



Tangle a seront obtenues, Temploi de tels resultats deviendra plus approprie. 

1.2.4 Ralentissements apparents et intrinseques 

Les pulsars se deplacent par rapport au systeme solaire. lis sont done animes d'un mouvement ap- 
parent, le mouvement propre, detectable si leur deplacement est suffisamment rapide et leur distance 
suffisamment faible. Dans le plan du ciel, un pulsar d' ascension droite a et de declinaison 5 possede 



un mouvement propre transverse egal a /ij- = y /i^ + ou = Of cos(5) et jUg = 5. 

Une consequence de ce mouvement propre transverse est T effet Shklovskii ( Shklovskii|1970 ) : la 



distance projetee du pulsar par rapport au systeme solaire croit, ce qui a pour consequence le fait que 
le taux de ralentissement apparent P est superieur au taux intrinseque, d'un facteur : 

APshk = 4x/' (Lll) 
ca 
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ou vt est la vitesse transverse, donnee par vt = fir x d et d est la distance de I'objet. L' equation 



LI 1 pent done se reeerire en : 



AftM='i^xP.2,43xlO--(^^) (1.12) 

c \mas yr ^ / \kpe / V ^ / 

Pour obtenir le taux de ralentissement intrinseque Pi du pulsar, il faut done tenir compte de I'effet 



Shklovskii en soustrayant APshk au taux apparent P. Comme le montre I'equation 1.12 I'effet n'est 
visible que pour les pulsars proches et a faible valeur de P. Pour les pulsars concemes, il pent avoir 
des consequences importantes sur les quantites derivees donnees par|LT}|l3]et|L8j Par exemple, PSR 
J0437-4715 (P = 5,757 ms, P = 5,73 x lO^^o^ J = 0, 16 kpc, = 141 mas yr^^) a un taux de 
ralentissement intrinseque Pi = 1, 39 x 10^*^ inferieur a P d'un facteur quatre, diminuant d'autant sa 
perte d'energie due au freinage E. 

Plus generalement, toute acceleration a le long de la ligne de visee conduit I'observateur a observer 
un ralentissement /* different du taux intrinseque d'un facteur AP = aP/c; c'est notamment le cas pour 
les pulsars appartenant a des amas globulaires, soumis a des forces gravitationnelles locales. Pour ces 
objets il faut done prendre les effets d' acceleration locaux dans le calcul de E, T, Bs et B^c- 



1.3 Pulsars milliseconde 



1.3.1 Differents types de pulsar 

On a vu dans les paragraphes precedents que les pulsars perdent de I'energie au cours du temps. 



ce qui se traduit par une diminution de leur periode de rotation. La Figure L3 montre le diagramme P 
- P, representant le taux de ralentissement des pulsars en fonction de leur periode, pour ceux dont on 
dispose d'une mesure de P liable, parmi les quelques 1900 recenses par le catalogue ATNF. 



Comme le montre la Figure L3 la majorite des pulsars connus ont une periode de rotation P 
comprise entre 0,1 et quelques secondes, et un ralentissement P entre 10^^^ et 10^^^. Ce sont les 
pulsars dits « normaux », ou « ordinaires ». On estime qu'ils naissent avec une periode comprise 
entre 14 et 140 ms ( Kramer et q/.|2003 ). La majorite des pulsars normaux est relativement jeune : de 
quelques milliers a quelques dizaines de millions d'annees typiquement. On constate I'existence, dans 
le coin inferieur droit du graphe, d'une zone depeuplee de pulsars, nommee « cimetiere ». Le faible 
nombre de pulsars detectes dans cette gamme de (P, P) amene a penser qu'il doit exister une limite a 
partir de laquelle les pulsars, apres avoir vu leur periode de rotation P croitre et leur ralentissement P 
diminuer au cours de leur vie, cessent d'emettre un rayonnement detectable a grande distance. 

Pour les six pulsars jeunes pour lesquels une mesure de I'indice de freinage n est disponible (cf. 
Kaspi & Helfand (2002) pour une revue des resultats), la Figure L3 montre le comportement de P en 



fonction de P, simule sur 10^ ans, a partir de I'equation differentielle n = = 2 — PP/P^ ou 

Ton a suppose que n est constant. Les parcours ainsi calcules sont quasi-lineaires, avec une pente de 
— 1 si I'indice de freinage est de 3. On constate que les six vecteurs ne sont pas diriges vers la region 
la plus densement peuplee en pulsars, signifiant vraisemblablement que I'indice de freinage n evolue 
au cours du temps afin que les pulsars jeunes viennent peupler cette region. On remarque que I'un des 
pulsars voit son ralentissement P augmenter, alors que ce taux diminue chez les cinq autres pulsars. 
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P(s) 

Figure 1.3 - Ralentissement P en fonction de la periode P pour les pulsars du catalogue ATNF en dehors 
d'amas globulaires et avec P positif. Les valeurs de P sont corrigees de I'effet Shklovskii lorsque possible. 
Les pulsars isoles sont representes par des points tandis que les pulsars appartenant a des systemes multiples 
sont representes entoures d'un cercle. Des lignes d'age caracteristique T constant sont representees, ainsi 
que des lignes de champ magnetique caracteristique Bs constant et des lignes de E constant. La majorite 
des pulsars connus sont les pulsars « ordinaires », au centre de ce graphe. En has a gauche du graphe on 
trouve les pulsars « milliseconde ». En haut a droite, on trouve les « magnetars », des pulsars au champ 
magnetique tres intense. Pour les six pulsars dont I'indice de freinage n = Q.Q./Q? a ete mesure, des traits 
pleins montrent le comportement simule de P en fonction de P sur 10^ ans, dans I'hypothese oil n est 
constant. 



II s'agit du pulsar de Vela, dont I'indice de freinage est Inferleur a 2, ce qui Impllque P > (Lyne 



et al. 1996[ ). De tels objets pourralent etre a I'orlglne des « magnetars », pulsars ayant les champs 



magnetlques les plus Intenses ( |Lorlmer|2009[ ). 

Enfin, on trouve dans la partle Inferieure gauche du graphe les pulsars « milliseconde » (MSPs), 
qui se demarquent de la majorite des pulsars ordinaires par une rotation extremement rapide (P < 
30 ms) et stable (P < 10"^^). Le plus rapide des MSPs connus a ce jour est PSR J1748-2446ad, 
decouvert dans I'amas globulaire Terzan 5, ayant une periode P ~ 1 , 396 ms ( Hessels et al. [2006 ) ! En 
outre, on remarque que les MSPs appartiennent a des systemes multiples pour la plupart, au contraire 
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des pulsars ordinaires. 



1.3.2 Formation 

Le modele du dipole tournant conduit a des ages caracteristiques pour les MSPs de I'ordre de 
10^ a 10^^ annees, soit superieurs a I'age de I'univers dans certains cas. Par consequent, le schema 
d'evolution decrit precedemment pour les pulsars normaux ne permet pas d'expliquer I'existence des 
MSPs. La Figure [13] indique par un cercle les pulsars appartenant a des systemes multiples. Comme 
on I'a deja evoque, il est frappant de constater que moins d'un pourcent des pulsars normaux appar- 
tiennent a des systemes multiples, alors que plus de 80% des pulsars milliseconde sont concemes, et 
tandis que ce taux est bien superieur a 1% pour les etoiles en general. Cette faible valeur s'explique 
par I'importante probabilite de rupture du systeme, du fait de la violence du phenomene de super- 
nova. Pour les systemes qui survivent a la supernova, il a ete propose que si I'etoile compagnon vient 
a deborder du lobe de Roche du pulsar, celui-ci peut acquerir une partie de la masse et du moment 
cinetique, accelerant ainsi sa rotation ( |Alpar et aL\\l9S2) . La decouverte recente d'un pulsar millise- 



conde qui a ete entoure d'un disque d' accretion dans un passe proche tend a confirmer ce scenario 
( Archibald et a/. [2009 1. Pour cette raison, on parle egalement de pulsar « recycles ». 



Pendant la phase d' accretion, il y a production de rayonnement X thermique par la matiere tom- 
bant vers I'etoile a neutrons, ce qui rend le systeme detectable en tant que « binaire a rayons X » 
(X-ray Binary). On en distingue deux categories, selon la masse de I'etoile compagnon : les systemes 
a haute masse (HMXBs) et a faible masse (LMXBs). Dans le premier cas, il se peut que I'etoile com- 
pagnon soit suffisamment massive pour evoluer en supernova, puis en etoile a neutrons ; si le systeme 
survit a la seconde explosion, on obtient un double pulsar. On connait neuf de ces systemes, le pre- 
mier decouvert etant PSR B1913-I-16 ( Hulse & Taylor|1975 1. Plus generalement, la phase d'accretion 



dans les HMXBs peut etre interrompue par I'explosion du compagnon, ce qui produit des pulsars 
« moderement recycles » (mildly recycled). Pour ce qui est des LMXBs, ou cette fois le compagnon 
est peu massif, la periode d'accretion n'est pas interrompue par une supernova, ainsi la phase d'accre- 
tion est beaucoup plus longue. II en resulte generalement un pulsar milliseconde a tres courte periode 
(typiquement, quelques ms), en orbite avec une naine blanche. 

II existe neanmoins des MSPs isoles et a tres courte periode, pour lesquels le recyclage n'a pour- 
tant pas ete interrompu. Dans les amas globulaires, ce phenomene peut s'expliquer par la forte densite 
d' etoiles : les interactions entre etoiles y sont plus nombreuses et les systemes multiples peuvent se 
defaire. En revanche la cause est moins bien connue pour les pulsars du champ galactique. II a ete 
propose que I'etoile compagnon puisse etre evaporee par I'intense vent de particules produit par le 
MSP dRudermang? a/.|1989D. 



1.3.3 Nombre et repartition 

La catalogue ATNF recense 169 pulsar^ ayant une periode inferieure a 30 m^ Parmi ces 169, 
97 appartiennent a des amas globulaires, les autres sont dans le « champ galactique ». En realite, si 
Ton tient compte des pulsars decouverts recemment et dont les decouvertes n'ont pas ete publiees, le 

^Enjuin 2009. 

"^Cette limite superieure sur la periode des pulsars milliseconde est une des valeurs couramment utilisees dans la 
litterature. La condition f < 50 ms est egalement utilisee. La limite de 30 ms, bien qu'arbitraire, presente I'avantage 
d'exclure le pulsar du Crabe (P ~ 33 ms) aux proprietes tres differentes de celles des MSPs, en particulier son taux de 
ralentissement P beaucoup plus grand. 
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nombre total de MSPs connus a ce jour s'eleve a un peu plus de 200. En particulier, bien que quelques 
pulsars soient en attente de confirmation, de nombreux pulsars milliseconde des amas globulaires ne 
sont pas recenses par le catalogue ATNlj^ Les listes des pulsars milliseconde du champ galactique et 
des amas globulaires et certaines de leurs caracteristiques sont donnees en Annexes |A.l 



et 



En moyenne, les MSPs connus sont moins lumineux dans le domaine radio que leurs alter-ego 
normaux ( [Kramer etal.\l99E . Bailes et a/.|1997 ). II en resulte par construction que les MSPs detectes 
dans le champ galactique sont plus proches que les pulsars normaux en moyenne. A cause de leur 
proximite, de nombreux MSPs paraissent etre situes en dehors du plan galactique (soit \b\ > 10°). Les 
amas globulaires se situant dans le halo de la Galaxie, ils apparaissent egalement en dehors du plan, 
en coordonnees galactiques. Pour ces differentes raisons les MSPs sont majoritairement situes a haute 



latitude galactique, comme le montrent les Figures 1.4 (A) et (B). 

Moins d'une centaine de pulsars milliseconde galactiques sont done recenses. Les etudes de po- 
pulations de MSPs emetteurs radio et des LMXBs predisent de s taux de naissance de pulsars milli- 
seconde compris entre 3 et 6,5 x 10^ par an (Ferrario & Wickramasinghe||2007 ; Story et aZ]|2007 



Kiel & Hurley, 2006 ). En prenant I'age de la Galaxie egal a 13,6 milliards d'annees, on obtient une 
population galactique de MSPs comprise entre 40000 et 90000 objets. 



1.3.4 Profil d' emission 



La Figure L5 montre quelques exemples de profils radio de pulsars milliseconde. On constate 
une certaine diversite parmi ces courbes de lumiere, dans la multiplicite des pics, la complexity des 
structures, leur largeur et I'ecartement des pics (pour les cas oh le nombre de pics est superieur a 
un). Les courbes de lumiere aux differentes longueurs d'onde representent les empreintes des pulsars, 
et donnent une idee de la configuration geometrique (inclinaison de I'axe magnetique par rapport 
a I'axe de rotation, angle de visee par rapport a ces axes) que peuvent completer des etudes de la 
polarisation de remission ( Ord et ar||2004 , Stairs et al. 1999, Xilouris et q/.||1998] ). Par exemple, le 
pulsar J1939+2134 presente deux pics, separes d'environ 180°. On peut en deduire que Ton observe 
les deux poles magnetiques au cours d'une seule rotation, ce qui correspond a une configuration ou la 
ligne de visee et I'axe magnetique sont proches, et sont orientes orthogonalement a I'axe de rotation. 

Les profils radio des MSPs sont en regie generale semblables a ceux des pulsars normaux. Comme 
on I'a deja evoque, remission radio est supposee produite par les particules chargees en deplacement 
le long des lignes de champ ouvertes. Le cylindre de lumiere, de rayon cP/ {2n) est beaucoup plus 
petit chez les MSPs, ce qui se traduit par une ouverture de faisceau plus large en moyenne, par rap- 
port aux pulsars normaux. Differents modeles phenomenologiques ont ete proposes pour expliquer 
la morphologic des profils radio des pulsars en general, et la diversite de ces profils, dans un schema 
simple de dipole magnetique tournant. D'apres le modele « coeur et cone », le faisceau radio possede 
un coeur entoure d'un ou plusieurs cones evides, ces differentes composantes etant centrees sur I'axe 
magnetique (Rankin 1983| ). L alternative principale est le modele « a taches », pour lequel le cone 
delimite par les lignes de champ ouvertes est peuple de taches d'emission, distribuees de fa9on non 
uniforme ( |Lyne & Manchester|1988j ). Selon des travaux recents, les donnees sont le mieux reproduites 
par un modele hybride utilisant des taches d'emission dans une structure ay ant un coeur et des cones 
( Karastergiou & Johnston|[2007[ ). On constate cependant d'apres la Figure L5 que certains profils de 
MSPs sont extremes de par leur complexity, par exemple celui de J2124— 3358. Cette complexite est 
parfois attribuee a 1' existence d'un champ multipolaire a proximite de la surface ( Kramer et al. [1997 ). 



cf. http ://www.naic.edu/ pfreire/GCpsr.html 
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Figure 1.4 - (A) : latitude galactique b en fonction de la distance pour les pulsars du catalogue ATNF. (B) : 
Carte du ciel en coordonnees galactiques. Les pulsars normaux sont representes par des points, les pulsars 
milliseconde sont entoures d'un cercle. 
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J1744-1134 J1824-2452A J1939+2134 J2124-3358 



Figure 1.5 - Exemples de profils radio de pulsars milliseconde. Ces profils ont ete enregistres a Nangay, a 
la frequence de 1400 MHz. 



1.3.5 Les pulsars milliseconde en physique et en astronomie 

Pour terminer ce chapitre introductif sur les pulsars en general et les pulsars milliseconde en 
particulier, repla^ons ces objets dans un contexte plus general que I'etude de leur emission en rayons 
7, qui sera abordee dans la suite. L' etude des etoiles a neutrons touche un panel large de domaines de 
la physique fondamentale et de I'astrophysique. De par leur nature extreme, les MSPs ont une place 
toute particuliere dans les differentes decouvertes issues de 1' astronomie des pulsars. Citons quelques 
exemples d'avancees realisees grace a I'etude des pulsars (liste non exhaustive) : 

- Les premieres exoplanetes ont ete decouvertes via la chronometrie du pulsar milliseconde 
B 1257+ 12, en 1990 a Arecibo ( |Wolszczan & Frail|1992D . 



Des tests fins de la theorie de la relativite generale ont pu etre realises grace aux doubles pulsars 



tels que B 1913+16, 1'orbite observee etant remarquablement conforme aux predictions (Hulse 



1994). Un prix Nobel a ete deceme a Hulse et Taylor en 1993 en recompense de ces travaux. 
Plus recemment, 1' observation du systeme double J0737— 3039 dans lequel les deux pulsars 
sont visibles a fourni les contraintes les plus fortes dont on dispose a ce jour sur la theorie de la 
relativite : par exemple la valeur du parametre orbital « post-keplerien » s mesuree est conforme 
a celle predite par la relativite generale, avec une precision de 0,05% ^ Kramer et al. [2006 1 . 



La dispersion du signal radio emanant des pulsars donne la possibilite de cartographier la den- 
site d' electrons dans la Galaxie. De la meme maniere, I'etude de la polarisation des pulsars 
permet de sonder le champ magnetique regnant dans la Galaxie ( Han et q/. [2006 1. 



- Une limite superieure sur la variation temporelle de la constante de gravitation G a pu etre 
determinee, grace au pulsar J0437— 4715 ( Verbiest et a/.|2008 ). 

- La chronometrie des pulsars appartenant a des systemes binaires, pulsars milliseconde pour 



la plupart, a permis la mesure de nombreuses masses d'etoiles a neutrons (Lattimer & Pra 



|kash 2007 1 ). Ces mesures de masses permettent de contraindre I'equation d'etat de la matiere 



super-dense des etoiles a neutrons. D'autres contraintes sont fournies par la detection de pulsars 
milliseconde ultra rapides : le plus rapide connu a present est J1748— 2446ad, dont la frequence 
de rotation est de 716 Hz ( [Hessels et q/.|2006) ). 



Tandis que les telescopes actuels operant du domaine radio aux rayons y tentent de percer cer- 
tains mysteres sur les pulsars - combien sont-ils dans la Galaxie, comment se forment-ils, comment 



evoluent-ils, quels sont leurs mecanismes d' emission - de grands projets tels que LOFAR et SKA (van 
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Leeuwen & Stappers]|2008[ |Smits et 'aL\\2009^ sont a venir. En plus d'augmenter considerablement 



le nombre de pulsars connus (plusieurs milliers de detections sont attendues) ce qui permettra une 
bien meilleure comprehension de ces objets, ils pourraient mettre en application une idee ambitieuse : 
utiliser un reseau de pulsars milliseconde comme d'un detecteur d'ondes gravitationnelles. 
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Les pulsars en rayons 7 
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Resume 

NOUS presentons ici les differents modeles de production de rayonnement 7 par les pulsars, 
ces modeles se repartissant en deux grandes families : Polar Cap (PC) d'une part, et Outer 
Gap (OG) d' autre part. La difference principale entre les differentes theories d'emission 7 est le 
site de production du rayonnement : pres des poles magnetiques ou au contraire en altitude dans la 
magneto sphere. La geometric de remission ainsi que les mecanismes variant d'un modele a 1' autre, 
les proprietes predites sont differentes et le LAT pent permettre de discriminer les modeles. Nous 
nous penchons par la suite sur le cas, specifique a priori, des pulsars milliseconde. Ceux-ci ont 
notamment un champ magnetique bien plus faible que celui des pulsars ordinaires. Les modeles 
theoriques prevoient neanmoins une dizaine de detections de MSPs par le LAT dans sa premiere 
annee. Dans la suite de ce chapitre, nous passons en revue les proprietes des pulsars 7 connus. A ce 
sujet, EGRET a bord de CGRO (Compton Gamma Ray Observatory) a ete la principale source de 
connaissances avant le lancement de Fermi. Des tendances se degagent parmi les profils d'emission 
7 ainsi que les proprietes spectrales des pulsars detectes par EGRET. Un point commun important 
est leur grand E, et c'est sur ce critere qu'a ete basee la campagne de suivi des pulsars par les 
telescopes radio et X, evoquee dans le chapitre suivant. Dans la fin de ce chapitre nous presentons les 
motivations de I'etude des pulsars dans cette gamme d'energie. 
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II. 1 Modeles theoriques 

Les premieres sources 7 identifiees etaient les pulsars du Crabe et de Vela, detectes par les teles- 
copes spatiaux SAS-2 et COS-B dans les annees 70, plusieurs decennies avant le lancement de Fermi. 
Bien que le rayonnement ydes pulsars ait ete etudie depuis cette epoque, le mecanisme d'emission 
est toujours debattu. 

Les modeles existant actuellement ont en commun de necessiter une region faiblement peuplee 
en particules chargees dans la magnetosphere des pulsars. En presence d'un tel deficit de charge par 
rapport a la distribution de charges a I'equilibre dans la magnetosphere, un champ electrique parallele 
^ll au champ magnetique se developpe, et celui-ci peut accelerer des particules chargees vers les 
hautes energies. Ces particules produisent par la suite des photons 7, par rayonnement de courbure, 
rayonnement synchrotron ou diffusion Compton inverse. 

Le disaccord entre les differentes theories porte justement sur la region d' acceleration des parti- 
cules chargees. Historiquement, deux grandes families de modeles se sont confrontees, ceux de la « 
calotte polaire » (Polar Cap) d'une part, pour lesquels I'acceleration de charges se produit au-dessus 
des poles magnetiques, ceux de la « cavite exteme » (Outer Gap) d' autre part pour lesquels les pho- 
tons 7 sont emis en altitude dans la magnetosphere. Des evolutions de ces theories ont recemment 
amene a considerer des zones d'emission intermediaires. Pour plus de details relatifs aux differents 



scenarios d'emission 7, le lecteur pourra par exemple se reporter a |Harding| ( |2007a[ ); |Hirotani| ( |2008[ ). 



II. 1.1 Modeles Polar Cap et Slot Gap 

On a vu dans le chapitre precedent qu'il existe un rayon au-dela duquel les lignes de champ 
magnetique ne peuvent boucler, c'est le rayon du cylindre de lumiere R^c = cP j (2;r). Les dernieres 
lignes de champ fermees delimitent, au niveau de la surface de I'etoile, la calotte polaire. Pour les mo- 
deles Polar Cap (PC), I'acceleration des particules chargees se produit sur quelques rayons stellaires 
au-dessus de cette calotte polaire, le long des lignes de champ magnetique ouvertes. 

Les modeles PC se repartissent eux-memes en differentes families, selon le mode d'approvision- 



nement en particules chargees a accelerer (Arons & Scharlemann 1979; Ruderman & Sutherland 



197 5). Pour les modeles a cavite vide (vacuum gap), des particules chargees sont piegees dans les 
couches superficielles de I'etoile a neutrons, et une zone vide de charge se forme au-dessus de la 
calotte polaire. Dans les modeles SCLF (Space-Charge Limited Flow), des particules chargees sont 
emises a la surface vers la magnetosphere, ou la densite de charges est insuffisante pour annuler le 
champ electrique parallele. Dans tous les cas, la region situee au-dessus du pole presente un deficit de 
charges, induisant I'existence d'un champ electrique intense E\^\^ en son sein, parallele aux lignes de 
champ magnetique ouvertes. Des particules chargees y sont done accelerees parallelement au champ 
magnetique, et emettent des rayons 7 par rayonnement de courbure ou par diffusion Compton inverse. 
A proximite de la surface de I'etoile, le champ magnetique est intense. Par consequent, les photons 7 
creent des paires electron / positron par interaction avec le champ magnetique, dans un zone appelee « 
front de formation de paires » (Pair Formation Front, PFF). Une cascade de paires electron / positron 
se developpe dans cette region de I'espace. Au-dela de I'epaisseur optique du PFF, le champ elec- 
trique parallele est completement ecrante. Ainsi les particules chargees qui s'en echappent ne sont 
plus accelerees, mais peuvent a leur tour emettre du rayonnement, formant le faisceau d'emission 7. 
Des evolutions du modele PC ont con duit a I'elaboration du modele « a fentes » (Slot Gap, SG) 



(Arons 1983 ; Muslimov & Harding 2003 1. A proximite de la surface, la repartition des charges cree un 
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Figure II.l - Schema de la magnetosphere des pulsars, representant les regions d' acceleration des particules 
chargees dans les differents modeles theoriques : calotte polaire (Polar Cap), cavite a fentes (Slot Gap) et 
cavite exteme {Outer Gap). Credit ; Alice Harding. 



champ electrique E\\ ^ 0, qui est ecrante plus haut dans un front de formation de paires. Cependant la 
frontiere de la calotte polaire est parfaitement conductrice, imposant = 0. Par consequent, le champ 
electrique au-dessus de cette frontiere est moins intense qu'au-dessus du pole. Les particules chargees 
ont done besoin d'une plus grande distance d' acceleration pour emettre des photons 7 energetiques 
pouvant induire des cascades de paires. Le PFF est done recourbe au-dessus des bords de la calotte, 



comme indique par la Figure II.2 Cette courbure cree des fentes entre le PFF et les demieres lignes 
de champ magnetique ouvertes, ou le champ magnetique £ | n'est pas ecrante, et ou par consequent 
les particules chargees peuvent accelerer jusqu'a haute altitude. Conduites dans la magnetosphere 
exteme, ces charges emettent des photons 7 par rayonnement de courbure ou synchrotron. 

Dans une approche purement geometrique qui a conduit au modele Two-Pole Caustic (TPC) (IDyks 



& Rudak|2003| ), il a ete montre que les profils d'emission des pulsars 7connus correspondent relative- 



ment bien a un rayonnement produit le long des lignes de champ ouvertes, depuis la surface de I'etoile 
jusqu'a haute altitude, a I'image de ce qui est predit par le modele a fentes. II est done interessant de 
remarquer que celui-ci pourrait servir de support physique a une description empirique. 
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Figure 11.2 - Schema de la region situee au-dessus de la calotte polaire. Le front de formation de paires 
(PFF) est represente au-dessus du pole. L'interstice delimite par le PFF et les dernieres lignes de champ 
magnetique ouvertes est la region d' acceleration des particules chargees dans le modele Slot Gap. Figure 
extraite de Harding & Muslimov (2004[l. 



II. 1.2 Modele Outer Gap 



Les modeles a cavite exteme originaux (Cheng et al. 1986, Romani 1996) partent du principe 



que les charges positives et negatives sont distribuees separe ment dans la magnetosphe re, suivant le 
signe de Q..B, comme dans le modele de Goldreich & Julian ( Goldreich & Julian||l969 ). Des cavites 
pauvres en particules chargees peuvent se former dans la region delimitee par la surface de charge 
nuUe (definie par p oc Q„B = 0) et le cylindre de lumiere, que les charges provenant de la surface de 
I'etoile ne peuvent peupler. Dans ces cavites, la composante de champ electrique parallele au champ 
magnetique est forte, de sorte que les particules chargees qui s'y trouvent sont accelerees vers 
les hautes energies. Des photons 7 sont ainsi emis, et ceux-ci peuvent former des paires electron / 
positron par interaction avec les photons X thermiques provenant de la surface de I'etoile a neutrons. 
Des cascades de paires sont ainsi engendrees : les electrons et positrons sont acceleres, ils emettent des 
photons, qui eux-memes forment des paires, et ainsi de suite. La cavite est done peuplee de particules 
chargees, venant ecranter le champ electrique le long des lignes de champ magnetique et limiter le 
volume de la region d' acceleration. La taille de la cavite ainsi maintenue, ou le champ electrique 
parallele En n'est pas ecrante, definit la geometric de remission. 



II.1.3 Differences attendues : geometries et spectres 

Les trois modeles prevoient des zones d'emission des rayons 7 et des mecanismes de production 
differents. Les deux observables de I'analyse des pulsars, quelle que soit la longueur d'onde consi- 
deree, sont le profil temporel de remission, et la forme spectrale du rayonnement. Les observations 
de pulsars en rayons 7 permettent de confronter les predictions morphologiques et spectrales des 
differents modeles, et par consequent de les discriminer. 



La Figure II. 3 montre des simulations de profils d'emission 7 attendus par les differents modeles. 



pour le pulsar de Vela et pour differentes configurations geometriques (a savoir, a : angle entre I'axe 
magnetique et I'axe de rotation, et C, : angle entre la ligne de visee et I'axe de rotation). Pour le modele 
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PC, r acceleration de charges se produit au-dessus de la calotte polaire, et est delimitee par les lignes 
de champ ouvertes. L' emission radio nait egalement au-dessus des poles. Par consequent le modele 
PC n'autorise que de faibles ecartements radio / 7. Au contraire, des ecartements plus grands sont 
permis pour des emetteurs de type SG et OG, du fait de la zone d'emission 7 decentree. De plus, 
OG et SG prevoient davantage de pulsars 7 non detectables en radio en raison du desalignement des 
faisceaux d'emission aux differentes energies. La fraction de pulsars detectables en radio parmi les 
emetteurs 7 varie done d'une theorie a 1' autre. 

A ces differences intrinseques s'ajoutent des effets relativistes qui modifient la morphologic de 
remission. Le temps de parcours des photons 7 est un premier exemple : des photons emis a une 
certaine altitude dans la magnetosphere et d'autres photons emis plus tot ou plus tard a une alti- 
tude differente peuvent parvenir a I'observateur simultanement. L'etalement en phase de remission 
7 produite en differents endroits de la magnetosphere engendre la formation de caustiques dans les 
courbes de lumiere. A proximite de la surface de I'etoile, le champ gravitationnel est intense ; ainsi 
les champs electrique et magnetique sont distordus et le parcours des rayons lumineux courbes. En al- 
titude, a proximite du cylindre de lumiere, les champs subissent des retards plus ou moins importants 
selon la vitesse de rotation. Ces differents effets jouent sur la direction d'emission des photons 7 voire 



sur la geometric des zones d'emission, modifiant le profil d'emission vu par I'observateur (Grenier & 
Harding|2006l ). 



II est clair d'apres la Figure II.3 que pour un pulsar donne (avec un angle d'inclinaison magne- 
tique a donne), deux observateurs ayant des lignes de visees distinctes par rapport a I'axe de rotation 
peuvent observer des profils et des intensites lumineuses differentes, voire ne pas detecter les pulsars 
si I'orientation des faisceaux est defavorable. Ceci conduit a definir un facteur de correction f^, repre- 
sentant le fait que I'intensite du signal observe sous un certain angle de visee n'est pas necessairement 



representatif de remission rayonnee sur le ciel entier ( Watters et q/. [2009 1 : 



Dans cette expression, h designe le flux en energie emis dans la direction (a, ^, 4>). En pratique, ce 



facteur est inconnu a priori, mais comme le montre la Figure IL3 on peut I'estimer par des hypotheses 
sur la geometric des faisceaux d'emission. Typiquement, ce facteur de correction est d'environ 1 pour 
les modeles SG et OG, compris entre et 0,5 pour le modele PC. 

Du point de vue des spectres d'emission, le modele PC prevoit une coupure super-exponentielle 

(c'est-a-dire en oil Ec est I'energie de coupure etb> 1) du fait du mecanisme de production 

de paires electron / positron par interaction avec le champ magnetique. Pour SG et OG, les coupures 
sont simplement exponentielles (soit b = l). On peut done ecrire le flux de photons attendu de la fa9on 
suivante : 



dN 
dE 



NoE 



(11.2) 



Dans cette expression, E designe I'energie du photon, T est I'indice spectral avant la coupure, Ec 
est I'energie de coupure du spectre et b est le parametre de forme de la coupure. Le terme A^o est un 
facteur de normalisation. La quantite dN /dE represente le nombre de photons emis par le pulsar et 



detecte par unite de temps, d'energie, et de surface de collection. La Figure II.4 montre, pour le pulsar 
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Figure II. 3 - En haut a gauche : courbe de lumiere y du pulsar de Vela vue par EGRET entre 30 MeV 
et 10 GeV. Les trois autres cadrans presentent des simulations de la courbe de lumiere attendue pour Vela, 
dans differents modeles d' emission {Polar Cap, Slot Gap et Outer Gap) et pour differentes configurations 
geometriques (a, oil a est Tangle d'inclinaison de I'axe magnetique par rapport a I'axe de rotation, 
et C, est Tangle que forme la ligne de visee avec Taxe de rotation. Chacun des cadrans est separe en deux 
parties. Celle du haut represente T emission yproduite pour les differentes valeurs de phase rotationnelle 4> 
(ici entre —180° et 180° ou entre 0° et 360°) et d' angle de visee C,. Pour un observateur fixe, on fait une 
coupe dans ce plan avec une valeur de constante, symbolisee par les lignes pointillees. La partie inferieure 
montre la courbe de lumiere predite pour la valeur de C, correspondante. Figure issue de Harding ( 2007b I. 



de Vela, des predictions faites avant lancement du spectre vu par le LAT au bout d'un an d'activite, 
pour une emission de type Polar Cap ou Outer Gap. Contrairement a EGRET, le LAT devrait etre en 
mesure de distinguer entre une coupure simplement ou super-exponentielle. La forme des spectres, et 
de leur coupure en particulier, sont done d'autres discriminateurs disponibles. 

II.1.4 Emission 7 des pulsars milliseconde 

Un certain nombre de differences importantes existent entre les pulsars normaux et les MSPs : 
ces derniers ont un champ magnetique en surface typiquement 10^ fois inferieur, comme le montre la 
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Figure II.4 - Spectres du pulsar de Vela pour une observation simulee d'un an avec le LAT, dans I'hypothese 
des modeles Polar Cap et Outer Gap. Les donnees d'EGRET sont superposees. Credit : Massimiliano 
Razzano. 



Figure [O] Par consequent, la creation de paires electron / positron par interaction des photons avec le 
champ magnetique a proximite de la calotte polaire est peu favorable. Pour les pulsars milliseconde 
on s' attend done a ce que le champ electrique au-dessus des poles magnetiques ne soit pas ecrante, 
a cause de 1' absence de cascades de creations de paires. Selon le modele PC, les electrons sont done 
acceleres depuis la calotte polaire jusqu'a haute altitude et emettent par rayonnement de courbure 
(Harding et al. 2005 p . Le modele prevoit de plus que les electrons puissent etre acceleres vers des 
energies plus grandes que pour les pulsars jeunes du fait de 1' absence de front de formation de paires. 
De ce fait remission de photons peut s'etendre au-dela de 10 GeV, avec une coupure simplement 
exponentielle (soit b — \ dans I'equation II.2), contrairement aux pulsars jeunes dans le modele PC. 
Une synthese de population basee sur la theorie de remission au-dessus de la calotte polaire a predit 
12 detections en rayons 7 de MSPs emetteurs radio par le LAT au bout d'un an (Story et al. 2007[ ). 
Enfin, une consequence de la non-existence de PFF au-dessus des poles est 1' absence de fentes au 
contact des demieres lignes de champ ouvertes, requises par le modele SG pour 1' acceleration de 
particules chargees et remission de rayons 7. 

L'emission 7 des MSPs a egalement ete modelisee dans le scenario OG ( [Zhang & Cheng|[2003 1. 
Dans cette etude on fait I'hypothese d'un champ magnetique multipolaire regnant a proximite de la 
surface de I'etoile a neutrons. Dans la cavite externe, le champ electrique parallele aux lignes de 
champ magnetique permet aux particules chargees de se deplacer, engendrant un courant de retour 
prenant la direction de la surface de I'etoile. Ces particules interagissent avec I'intense champ ma- 
gnetique multipolaire et emettent des rayons X, qui viennent alimenter le reste de la magnetosphere 
et en particulier la cavite externe. L'arrivee de rayonnement X dans la cavite externe recree les condi- 
tions decrites precedemment pour le modele Outer Gap applique aux pulsars normaux : les photons 
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X interagissent avec les photons yemis dans la cavite pour engendrer des cascades de paires electron 
/ positron, qui sont accelerees et emettent des photons 7, etc. L' emission de photons 7 se produisant 
dans la cavite exteme, elle peut done apparaitre decalee du rayonnement radio. Les photons 7 emis 
dans ce modele ne peuvent pas exceder quelques GeV, et la coupure attendue est simplement ex- 
ponentielle. Le nombre de detections de MSPs par le LAT dans sa premiere annee predites dans ce 
modele est egalement de I'ordre d'une dizaine ( Zhang et a/. [2007 1. 

Bien que les modeles ne s'accordent pas sur les sites et les mecanismes de production des photons 
7, on constate que ceux-ci predisent un nombre relativement important de detections de MSPs par le 
LAT. 



II.2 Observations de pulsars en rayons 7 

Avant le lancement de Fermi, I'essentiel des connaissances relatives aux pulsars emetteurs 7 a 
ete acquis grace au telescope EGRET, a bord de CGRO. Entre 1991 et 1997, EGRET a permis la 
detection de six pulsars en rayons 7, parmi lesquels le Crabe et Vela, et quelques candidats serieux, 
comme nous le verrons dans la suite de ce chapitre. En ajoutant la detection du pulsar B 1509— 58 
jusqu'a 10 MeV par le telescope COMPTEL, egalement a bord de CGRO, on comptait sept detections 
certaines de pulsars en rayons 7 apres la mission de CGRO, auxquelles se sont ajoutees recemment 
les deux detections effectuees par AGILE et de nouveaux candidats. 



II.2.1 L'heritage de CGRO 
Pulsars normaux 

L'echantillon de pulsars emetteurs 7 vus par EGRET et COMPTEL a bord de CGRO etait de sept 
objets et quelques candidats. On ne pouvait done pas tirer de conclusions definitives a partir de leurs 
proprietes. Un certain nombre de tendances se degageaient neanmoins. 

Ces sept pulsars sont relativement jeunes et energetiques (Thompson 2004| ). lis sont parmi les 



pulsars qui perdent le plus d'energie cinetique rotationnelle par freinage electromagnetique ; disposant 
ainsi d'un important reservoir a convertir en rayonnement 7. 



La Figure II.5 montre les courbes de lumiere des pulsars d'EGRET et de COMPTEL dans cinq 
bandes d'energie. Le profil d'emission et I'alignement en phase des composantes radio, optiques, X 
et 7 varient fortement d'un pulsar a 1' autre. On remarque par exemple que le pulsar du Crabe presente 
deux pics en alignement etroit dans les differentes longueurs d'onde. Le pulsar de Vela n'a qu'un pic 
visible en radio, mais deux en optique, en X et en rayons 7. Parmi ces sept pulsars, Geminga est le 
seul a n'avoir jamais ete detecte en radio jusqu'a present. Malgre ces differences de comportement en 
fonction de I'energie, on remarque qu'a I'exception de B1509— 58, les pulsars presentent deux pics 
d'emission 7. Dans certains cas, les plus frappants etant ceux de Vela et de Geminga, des photons sont 
egalement detectes dans une zone situee entre les deux pics (le bridge). Ces observations indiquent 
que remission 7 est produite dans une zone large, vraisemblablement evidee. 

Des similitudes se degagent egalement parmi les proprietes spectrales des pulsars d'EGRET. Pour 
un ajustement de leur spectre par une loi de puissance avec coupure exponentielle (soit A^o fi^^e^^^/^o^), 
les indices spectraux F sont typiquement inferieurs a 2. L'energie de coupure Eq quant a elle est ge- 



neralement comprise entre 1 et 4 GeV ( |Thompson|2004 ). On peut ecrire la luminosite 7de la maniere 
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Crab B1 509-58 Vela B1 706-44 B1951+32 Geminga 81055-52 
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Figure II.5 - Courbes de lumiere dans differentes longueurs d'onde pour les sept pulsars y detectes avec 
EGRET et COMPTEL. Figure issue de |Thompsonl ( |2008ap . 



suivante 



(11.3) 



ou h est le flux en energie detecte par le telescope, d est la distance du pulsar et /fj est le facteur de 



correction geometrique defini precedemment {cf. equation II. 1 ). Connaissant la luminosite 7, on peut 
calculer I'efficacite de conversion de la perte d'energie due au freinage en rayonnement 7, t] = Lj/E. 
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En I'absence d'hypothese sur la geometric des faisceaux 7, il apparait que la luminosite Ly(l sr) 



hd^ au-dessus de 100 MeV varie grossierement comme pour les pulsars d'EGRET, soit Tj v£ 
( Arons|1996 ). Cette propriete a ete renforcee depuis, sur la base de considerations theoriques (Walters 



et al. 



20091. L'efficacite de rayonnement diminue done, a mesure que E croit. Etant donne que les 



pulsars ne peuvent convertir plus de 100% de E en rayonnement 7, une autre implication de cette loi 
est que 77 doit saturer a has E. Une autre possibilite est qu'il existe une valeur de seuil pour remission 
7, situee en de9a de 3 x 10^"^ erg/s, valeur de E la plus faible parmi les pulsars detectes par EGRET. 

Pulsars milliseconde, et cas de PSR J0218+4232 

Une recherche de pulsations en rayons 7 a ete effectuee dans un lot de donnees de deux ans 
d' observations d'EGRET, pour un echantillon de 19 pulsars milliseconde ( Fierro et q/.||1995 I. Cette 
etude n'a permis aucune detection de pulsations en rayons 7, y compris pour J0437— 4715, parmi 
les meilleurs candidats car tres proche et de E comparable a ceux des jeunes pulsars 7 detectes par 
EGRElQ Cependant, les limites superieures a 3 C7 sur les flux observes etaient bien superieures aux 
valeurs de flux predites (cf. Table 4 de Fierro et al. ( 1995 )). 

En effectuant 1' analyse de la distribution spatiale des evenements detectes par EGRET autour de 
la source du troisieme catalogue 3EG J0222+4253, il est apparu que la position de la source entre 100 
et 300 MeV etait compatible avec celle du pulsar J02 18+4232, decouvert peu de temps auparavant 
(Kuiper et al. 2000 , [Navarro et al 1995[ ), apres la fin de la mission EGRET. Au-dela de 1 GeV, le 
blazar 3C66A situe a une distance angulaire de 1° du pulsar milliseconde est la contrepartie evidente 
pour remission 7, creant la confusion spatiale. L' analyse temporelle des photons entre 100 MeV et 
1 GeV a quant a elle abouti a un signal de significativite 3,5 o, avec un profil a deux pics separes de 
0,45, analogues a ceux des pulsars normaux. Plus tard, Kuiper et al. (2002) ont reanalyse les donnees 
EGRET conjointement avec des observations aux rayons X du telescope Chandra (Kuiper et al. [2002]). 



Les courbes de lumiere radio, X et 7 ainsi obtenues sont montrees dans la Figure II. 6 L'alignement 
apparent entre les differentes longueurs d'onde leur permit d'accroitre la significativite du signal 7de 
3,5 a 4,9 a, grace a des considerations statistiques basees sur l'alignement apparent entre pics 7et X. 



II.2.2 Observations de pulsars par le telescope AGILE 

Le telescope AGILE {Astro -rivelatore Gamma ad Immagini LEggero) a ete mis en orbite en avril 
2007, soit un peu plus d'un an avant Fermi. II a depuis observe le ciel 7 entre 30 MeV et 30 GeV, 
avec une sensibilite comparable a celle d'EGRET mais un champ de vue plus large (~ 2,5 sr) et une 
meilleure resolution angulaire. AGILE a dans un premier temps observe certains pulsars 7 connus 
prealablement : Vela, Crabe, Geminga et B1706— 44, obtenant ainsi des courbes de lumiere 7 plus de- 
taillees que celles d'EGRET Pellizzoni et al. ( 2009b[ ). Par la suite et pendant les premiers mois d'acti- 
vite de Fermi, de nouvelles detections fermes de pulsars emetteurs 7 ont ete realisees grace a AGILE : 
PSR J2021+3651 ( Halpem et oTpOO^ et J2229+6114 ( [PeUizzoni gFQr]|2009aD , en plus de quelques 
detections marginales : J1016-5857, J1357-6429, J1513-5908, J1524-5625, J2043+2740, et le 
pulsar milliseconde appartenant a I'amas globulaire M28, J1824— 2452A (nous reviendrons sur cette 
detection marginale a la fin de ce manuscrit). Comme pour EGRET, on s'aper§oit que ces pulsars 
comptent parmi ceux qui perdent le plus d'energie cinetique rotationnelle par freinage electromagne- 
tique. 



'En realite, I'effet Shklovskii n'avait pas ete pris en compte. Dans le cas de PSR J0437— 4715, la valeur de E intrin- 
seque est inferieure a la valeur apparente d'un facteur quatre. 
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2.0 



Figure 11.6 - Courbes de lumiere de PSR J0218+4232 a differentes energies, a : profil radio, a 610 MHz. 
b : observations en rayons X entre 0,08 et 10 keV, avec Chandra, c : phasogramme yobtenu avec EGRET, 
entre 100 MeV et 1 GeV. Figure issue de|Kuiper et al.\\2QQ2\. 



Aucun pulsar milliseconde ne figurait done parmi les detections formelles d' EGRET et AGILE 
aux rayons 7 au moment des premieres observations du ciel par le LAT. Cependant, la detection 
probable de PSR J02 18+4232 par EGRET faisait de ce pulsar une cible prioritaire pour le LAT, tout 
comme J1824— 2452A, possible detection d' AGILE. Ces observations suggerent la possibilite que 
d'autres pulsars de cette categoric soient des sources emettrices de rayons 7. 

II.2.3 Quels pulsars pour Fermi ? 

Empiriquement, les neuf pulsars detectes par EGRET, COMPTEL et AGILE ainsi que les dif- 
ferentes detections a faible significativite ont en commun de presen ter des taux de perte d'energie 



par freinage electromagnetique superieurs a 3 x 10-^^ erg/s (Thompson et al. 1999, Pellizzoni et al. 



2009a[ ). Cependant, la geometric des faisceaux d'emissions 7 etant mal connue, la relation entre la 



valeur de E et le flux 7 attendu, done la detectabilite du pulsar, est tres incertaine. Pour cette raison, 
une liste de cibles pour Fermi a ete definie en retenant les pulsars verifiant : 
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E > lO^^erg/s (11.4) 



Cette coupure conduit a une liste de 230 pulsar^ ( [Smith et a/. [2008 1. Les pulsars des amas globu 



laires pouvant avoir une valeur de E apparente exageree a cause d' accelerations dans des potentiels 
gravitationnels locaux, certains ont ete exclus, amenant a une liste de 218 pulsars a suivre a I'aide des 
dilferents radiotelescopes de par le monde, et les telescopes X pour les pulsars non emetteurs radio (au 
nombre de cinq dans cette liste : Geminga, J181 1-1925, J1846-0258, J0540-6919 et J0537-6910). 

Enfin, un critere a ete defini dans le but d'accorder une priorite plus ou moins importante au suivi 
des pulsars. On peut definir un « flux du au freinage » comme le flux d'une source de luminosite E 
situee a une distance d, soit E/d^. Cette quantite represente en quelque sorte le flux maximal pouvant 
etre emis sur I'ensemble du spectre electromagnetique. On a vu precedemment que pour les pulsars 
d'EGRET, la luminosite 7 au-dessus de 100 MeV semble croitre comme V£, soit 77 = Ly/E oc v^. 
On peut done ponderer le flux du au freinage par I'efficacite d'emission en rayons 7 empirique, et 
obtenir un flux 7 effectif : 

^■e<r=l| = ^ (11.5) 

Pour les pulsars ayant des valeurs de E equivalentes, la priorite est done donnee aux plus proches. 
Selon ce classement, les six pulsars d'EGRET se classent parmi les 20 premiers pulsars, interca- 
les avec quelques objets decouverts apres la fin de I'activite du telescope 7. Neanmoins quelques 
pulsars milliseconde pour lesquels EGRET n'avait pas detecte de pulsations, PSR J0034— 0534 et 
J0437— 4715, se classent parmi les dix premiers, soulignant a nouveau le fait que Fgff ne prend pas 
en compte tous les parametres, notamment la geometric d'emission variable. La Figure [1177] montre, 
pour les pulsars recenses par le catalogue ATNF, la valeur de F^^ff, normalisee a la valeur de Vela, en 
fonction de leur periode de rotation. Les pulsars cibles pour Fermi sont indiques par des carres et des 
triangles. Par construction, ils sont generalement distribues aux grandes valeurs de VS/d'^. 



II.3 Pourquoi etudier les pulsars en rayons 7? 

Dans la premiere partie de ce chapitre nous avons vu que differents modeles sont en competition 
pour I'explication de remission 7 des pulsars : Polar Cap, Slot Gap et Outer Gap notamment. Chacun 
d'entre eux propose une region particuliere de la magnetosphere des pulsars comme point de depart 
du rayonnement. On a vu par exemple que les modeles prevoient des coupures spectrales de formes 
differentes, des morphologies d'emission 7differentes et des ecartements plus ou moins importants 
avec la composante radio. Avec un echantillon plus large de pulsars emetteurs 7, on peut etablir des 
tendances dans les profils d'emission ainsi que dans les spectres, et ainsi placer des contraintes sur la 
geometric de remission dans la magnetosphere. L' augmentation du nombre de pulsars connus permet 
done de contraindre, de raffiner voire d'eliminer les theories proposees. 

Les telescopes predecesseurs du LAT, et en particulier EGRET, ont revele que les pulsars emet- 
teurs 7 sont ceux qui perdent le plus d'energie par freinage electromagnetique. Mais EGRET n'a pas 

nouveau, ce sont les valeurs apparentes de P qui ont ete utilisees. 
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Figure 11.7 - Flux de freinage effectif \ft,j(p- pour les pulsars du catalogue ATNF. Les flux sont normalises 
a celui du pulsar de Vela, 3,13 x lO'^ (erg/s) '^^/t/^. Les pulsars de valeur de E comprise entre 10^* erg/s et 
3 X 10^* erg/s sont indiques par des triangles, ceux verifiant £ > 3 x 10^^ erg/s sont montres par des carres 
vides. Les carres pleins sont les pulsars detectes par EGRET, COMPTEL et AGILE, ainsi que les detections 
marginales. 



permis de comprendre combien remission en 7 est un phenomene general pour les pulsars a grand E. 
En particulier, on peut s'attendre a ce que d'autres parametres entrent en jeu, I'age ou I'intensite du 
champ magnetique en surface ou au cylindre de lumiere. Le scenario d' evolution peut egalement avoir 
son importance : existe-t-il des pulsars milliseconde emetteurs de rayons 7? Si c'est le cas, le pro- 
cessus d' emission est-il commun aux populations « jeunes » et « recyclees » ? Ces questions ont leur 
importance pour les syntheses de population ( jGonthier et al. [2002 ; Story et q/. [2007 1, qui cherchent a 
determiner le nombre de pulsars dans la Galaxie et leurs proprietes. Enfin, la non-detection de pulsars 
a grandes valeurs de E est tout aussi importante dans la mesure oil elle renseigne sur la geometric des 
faisceaux aux differentes longueurs d'onde. 

Enfin, si I'etude des pulsars en rayons 7 touche aux modeles d'emission theoriques et a revo- 
lution des pulsars dans un contexte plus general, il est bien evident que 1' identification de pulsars 
7 peut permettre la decouverte et I'etude de sources avoisinantes. De meme, on s'attend a ce que 
les pulsars eloignes, non detectables mais tout de meme emetteurs, contribuent a remission diffuse 
qui est source de bruit pour beaucoup d'etudes. La contribution des pulsars a remission diffuse est 
encore mal connue, mais des modeles predisent une contribution importante, en particulier par les 
pulsars milliseconde (voir par exemple [Faucher-Giguere & Loeb| ( |2009| ); |Wang et al] ( |2005| )). Les 
raisons mentionnees ici sont autant de motivations pour la recherche de pulsars 7 avec le Large Area 
Telescope a bord de Fermi. 
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Partie II 



Instrumentation radio, X et 7 



49 



50 



Chapitre III 

Chronometrie des pulsars 
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Resume 

LA specificite des pulsars est leur emission cadencee, par consequent 1' analyse temporelle est 
centrale dans toute etude de ces objets. Pour construire les courbes de lumiere des pulsars, 
c'est-a-dire le signal emis en fonction de la phase rotationnelle, il faut connaitre la frequence de 
rotation (ou de fa9on equivalente, la periode) en fonction du temps. Connaissant cette frequence, on 
peut calculer les phases pour un lot de dates donne, et ainsi acceder a la morphologic de remission. 
La mesure des parametres de rotation, comme par exemple la frequence /o et ses premieres derivees, 
constitue la « chronometrie » des pulsars, que nous detaillons dans la deuxieme partie de ce chapitre. 
EUe consiste en plusieurs etapes, parmi lesquelles I'enregistrement de « temps d'arrivee » de 
pulsations, le transfert de ces temps vers le barycentre du systeme solaire, et le calcul des residus 
de datation. L'aboutissement de la chronometrie est la construction d'ephemerides, qui permettent 
de calculer la phase rotationnelle dans un intervalle de temps donne. Dans la troisieme partie nous 
revenons sur la campagne de datation des pulsars candidats a la detection en rayons 7 par le LAT. En 
particulier, nous detaillons I'algorithme de construction et de mise a jour des ephemerides construites 
a partir des donnees du radiotelescope de Nan^ay. Les ephemerides ont ete centralisees dans une 
base de donnees accessibles par les pulsaristes de la collaboration Fermi-LAT. Outre la chronometrie, 
les instruments impliques dans la campagne de datation ont egalement permis la verification des 
outils d' analyse temporelle mis en place pour les donnees du LAT, ce que nous presentons dans la 
quatrieme partie du chapitre. 
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III.l Courbes de lumiere de pulsar 

Pour un observateur fixe, les pulsars emettent un signal module par leur frequence de rotation. 
Construire la courbe de lumiere (ou phasogramme) d'un pulsar consiste a decouper le signal accu- 
mule pendant un certain temps d'observation en morceaux correspondant a une rotation de I'objet 
sur lui-meme, puis les empiler. On a vu precedemment que les pulsars ralentissent, que leur taux de 
ralentissement varie, et ainsi de suite. Une fa9on d'exprimer la frequence de rotation / d'un pulsar en 
fonction du temps est alors d'utiliser un developpement en serie de Taylor, que Ton tronque a I'ordre 
approprie : 



/W-/o+/ox(?-ro) + ^x(?-ro)2 + ... (III.l) 

ou /o, /o et /o sont la frequence de rotation et ses premieres derivees a une date de reference Tq. 
Pour un signal continu puis echantillonne, tel qu'en radio, ou un signal discret re§u sous la forme de 
photons espaces dans le temps, on dispose de dates pour lesquelles on souhaite connaitre la phase 



rotationnelle. D'apres I'equation III.1[ cette phase est donnee par : 



4>(0 = 4>o+/ox(?-ro) + ^x(?-ro)2 + ^x(r-ro)3 + ... (111.2) 

Dans cette expression, 4>o est la phase (arbitraire) du pulsar a la date Tq. La partie entiere de 4>(?) 
donne le nombre de tours du pulsar effectues entre les temps Tq et t. La partie fractionnaire indique le 
stade de la rotation, au cours d'un tour. 

Pour construire la courbe de lumiere il faut done connaitre /o et un certain nombre de ses derivees 
avec precision. Le besoin de precision est d'autant plus important que le signal est faible : plusieurs 
semaines d'observations etaient necessaires a EGRET pour detecter certains pulsars en rayons 7, 
ce qui pouvait representer plusieurs millions a milliards de rotations des etoiles sur elles-memes. 
Comme nous le verrons dans la suite de ce chapitre, une connaissance approximative des parametres 
de rotation a pour consequence le decalage progressif du profil ou son etalement, empechant in fine 
la detection du pulsar. 

Pour cette these, consacree a la recherche de pulsars milliseconde emetteurs 7 dans les donnees 
de Fermi, le besoin de precision sur f{t) est done crucial. Deux strategies sont envisageables pour 
rechercher les /o, /o, /o, etc. qui permettront eventuellement la detection en rayons 7 : 

- Une technique consiste a rechercher les parametres optimaux dans le signal lui-meme, en es- 
sayant un certain nombre de valeurs de /o et de ses derivees (et eventuellement d'autres pa- 
rametres, par exemple orbitaux si Ton recherche des systemes binaires), a la recherche de la 
configuration qui optimise la significativite du signal. Cette technique « a I'aveugle » a ete mise 



en oeuvre pour la recherche de pulsars dans les donnees de Fermi, avec succes (Abdo et al. 



2008[|Ziegler et a/.|2008j . Plus generalement, elle est surtout employee dans le domaine radio, 



domaine dans lequel presque tous les pulsars connus ont ete decouverts. 
L' autre methode, mise a I'oeuvre dans le cadre de cette these, est la « chronometrie » des 
pulsars, dans les domaines d'energie oil le signal est plus abondant (essentiellement en radio). 
Cette technique est presentee dans la suite de ce chapitre. 



52 



III.2. PRINCIPE DE LA CHRONOMETRIE 



III.2 Principe de la chronometrie 

Le but de la chronometrie (ou datation) d'un pulsar est la construction d'un modele decrivant 
la rotation de I'etoile en fonction du temps avec precision et dans un intervalle de temps donne. 
Pour cela, on commence par observer le pulsar a plusieurs reprises, et on determine pour chaque 
observation un « temps d'arrivee » (time of arrival, TOA), a savoir, une date a laquelle I'objet etait 
a un stade donne de sa rotation (par exemple, le point oii remission est maximale). Disposant de A'^ 
dates auxquelles la phase du pulsar etait identique, on recherche alors une fonctioii[^qui reproduit au 
mieux la rotation du pulsar. 



III.2.1 Temps d'arrivee de pulsations 

Si la morphologic de remission peut varier fortement d'une rotation a I'autre (iRramer et al. 1999 1, 
I'experience montre que des profils accumules sur une certaine duree sont des empreintes stables. Ces 
profils integres fi peuvent etre decrits par la relation : 



fi{^)=a + b^F{<p') + c{<p) (III.3) 

Dans cette equation, atib sont des constantes, c est une fonction representant la contribution du 
bruit, et F est un « profil ideal » de remission du pulsar, decale de = 0' — ^. L'enregistrement de 
TOAs revient a determiner le dephasage A^, pour chaque observation puis connaissant les dates 
des A'^ observations, en deduire les dates pour lesquelles le pulsar etait au meme point de sa rotation. 

Plusieurs choix sont possibles pour la definition du profil de reference F . Une solution simple est 
un pic de Dirac, auquel on aligne un point arbitraire des observations individuelles /,, par exemple 
le maximum de remission. Avec cette solution on omet cependant la morphologic de remission 
du pulsar. On peut creer le profil F en modelisant une observation par une somme de gaussiennes 
d' amplitudes et de phases centrales differentes. On peut egalement lisser une observation individuelle, 
c'est-a-dire reduire la contribution du bruit. Enfin, on peut creer un profil ideal de maniere iterative, 
en alignant observations entre elles grace a un modele de datation obtenu au prealable, et en les 
additionnant pour ameliorer le rapport signal-sur-bruit. 



III.2.2 Modele de datation 

Les TOAs determines precedemment sont generalement topocentriques, c'est-a-dire enregistres 
par un observatoire en mouvement par rapport a la source de remission, et dans une echelle de temps 
variant de fa9on irreguliere par rapport au temps propre du pulsar. Le barycentre du systeme solaire 
est, en premiere approximation, en mouvement inertiel par rapport aux pulsars de la Galaxie ; on 
choisit done de convertir les dates « topocentriques » en temps « barycentriques » (7]). L' equation 
de transfert reliant les et (7^) est la suivante : 



T =t + ^c-^D + ^RB+ Ae,q - As,Q + Ab (III.4) 

Dans cette equation, Ac represente les corrections d'horloges a apporter entre I'appareil de me- 
sure, qui enregistre des temps dans une certaine echelle, UTC (Temps Universel Coordonne), TAI 

'll n'y a pas necessairement unicite de la solution. 
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{Temps Atomique International) ou TT {Temps Terrestre), et une echelle standard proche du temps 
propre du pulsar, corrigee des irregularites dues a la rotation de la Terre. L' utilisation de 1' echelle 
TT, liee au temps atomique et dont 1' unite est la seconde SI, est preconisee pour cette derniere. Plus 



de details sur les definitions de systemes temps-coordonnees sont disponibles dans | Andersen| ( |1999] ) ; 
Rickman| ( |200T] ); [McCarthy & Petit] ( |2004l ) 



Le membre A© est un terme dispersif. Dans le milieu interstellaire, les rayons lumineux sont 
disperses par les electrons presents sur leur parcours. Le terme dominant est : 



DM 



^ ^ /sSB 



(III.5) 



ou DM est la « mesure de dispersion », correspondant a la densite de colonne d'electrons pour 
une direction donnee, /ssb est la frequence d'observation et ^ = 2,410 x 10 MHz^ cm pc s^ 
est la constante de dispersion (Manchester & Taylor 1977). Notons que pour un observateur situe 
sur Terre, la frequence d'observation effective varie au cours de I'annee, principalement a cause de 
I'effet Doppler du au mouvement de la Terre. La frequence /ssb de I'expression precedente est par 
consequent celle d'un observateur situe au centre de masse du systeme solaire. 

A cause de cet effet dispersif, des photons de frequences differentes mais e mis si multanement 
arrivent en decalage (une illustration de I'effet du terme A/) est montre en Figure IILl I. A I'inverse, 
il est crucial de maitriser cette correction si Ton veut mesurer I'ecart temporel entre des pics emis 
a des energies differentes, ces pics pouvant etre reellement decales dans le temps et done signifier 
des origines differentes dans la magnetosphere du pulsar. Corriger de cet effet revient a passer a 
une observation faite a frequence infinie. Remarquons enfin que la mesure de dispersion DM peut 
varier au cours du temps, pour les cas ou la densite de colonne d'electrons sur la ligne de visee evolue 



(Isaacman & Rankin 1977 , Ramachandran et al. 2006 1. Dans ce cas on peut ajouter a III. 5 des derivees 
temporelles DM, I)M, etc. La Figure 111.2 illustre I'effet d'une variation de la mesure de dispersion. 

Le membre A^ est appele terme de Roemer, et contient le temps de trajet classique de la lumiere, 
entre I'instrument de detection et le centre de masse du systeme solaire. On peut le decomposer de la 
maniere suivante : 



{rBT + rTt)rBP 



(III.6) 



oil rsT est un vecteur pointant du barycentre du systeme solaire vers le centre de la Terre, rj-t 
relie le geocentre et le telescope, et rsp est un vecteur unitaire donnant la direction du pulsar, depuis 
le centre de masse du systeme solaire. Pour calculer ces vecteurs, on a recours aux « ephemerides 
planetaires » telles que DE200 ( |Standish||1990[ ) et DE405 ( |Standish||1998[ ) publiees par le Jet Pro- 
pulsion Laborator^ Ces ephemerides donnent la position des principaux astres du systeme solaire 
en fonction du temps. Le terme de Roemer est largement dominant dans le transfert de dates vers le 
barycentre. Sa magnitude est approximativement (1 U.A.) x cos(/3)/c, ou /3 est la latitude ecliptique 
du pulsar. Pour un pulsar situe dans le plan de I'ecliptique (c'est-a-dire j8 — 0), I'effet est maximal et 
son amplitude vaut environ 500 s. 

Certains pulsars sont animes d'un mouvement apparent dans le plan du ciel, le mouvement propre. 
L' ascension droite et la declinaison varient au cours du temps, modifiant ainsi le vecteur rgp. Si le 



http ://www.jpl. nasa.gov 
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Figure III.l - Effet de la dispersion interstellaire, sur une bande de frequence de 288 MHz centree a 1380 
MHz. Les pulsations arrivent a I'observatoire avec un decalage en temps, variant comme I'inverse du carre 
de la frequence. Figure extraite de Lorimer & Kramer ( 2004| l. 



mouvement du pulsar n'est pas pris en compte, le TOA au barycentre est errone ; I'ecart avec le temps 
d'arrivee « vrai » correspond a I'erreur positionnelle, modulee par I'orbite terrestre. 

Enfin, I'effet de parallaxe de chronometriej^ induit par la variation de la courbure du front d'onde 
de r emis sion du pulsar au cours de I'orbite terrestre, est visible pour certains pulsars proches (I'equa 
tion 



III.6 fait en effet 1' approximation d'une onde plane). Son amplitude est donnee par (1 U.A.) 



X 



cos (j8)/ (2c J), et donne done acces a la distance d du pulsar, avec une grande precision dans certains 
cas ( [Lommen et al.\2006\ |Hotan et a/. 12006] ). 

Les termes A£.0 et A5 © sont des corrections relativistes. Ce dernier, I'effet Shapiro, correspond 
aux retards induits par la deformation de I'espace-temps a proximite des masses du systeme solaire 
( |Shapiro||1964[ ). II est le plus intense pour les photons passant pres du Soleil (I'effet est typiquement 



^Different de la parallaxe classique. 
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Figure 111.2 - Residus de chronometrie pour le pulsar milliseconde B 182 1—24 observe par le radioteles- 
cope de Nangay. Les croix indiquent des TOAs enregistres a 1,4 GHz, les carres sont des observations a 
2 GHz. Pour ce graphe une mesure de dispersion constante a ete utilisee. Les TOAs enregistres aux deux 
frequences d' observations s'ecartent. Dans ce cas il est necessaire de tenir compte d'un e variation de la 
mesure de dispersion, DM, dans le modele de datation. Figure issue de Smith et al. (2008 l 



de 100 /is). La contribution des planetes du systeme solaire varie de quelques ns a quelques centaines 
de ns. L'effet Einstein, represente par le terme ©, correspond a la dilatation de I'espace-temps 



au cours du mouvement de la Terre dans le potentiel gravitationnel du systeme solaire ( ,6 acker & 
Hellings|1986l ). 



Pour finir, le terme A5 correspond aux corrections a prendre en compte pour les pulsars apparte- 
nant a des systemes multiples. Ceux-ci sont animes d'un mouvement orbital, ce qui se traduit en des 
effets Roemer, Shapiro et Einstein locaux, analogues a ce qui concerne le systeme solaire. Plusieurs 
modeles permettent la description de ces orbites, selon leur excentricite ou leur complexite, selon le 
nombre de compagnons ou leur masse. Par exemple, le modele BT (Blandford-Teukolsky) qui est 
adapte aux orbites lentes obeissant a la mecanique newtonienne, DD (Damour-Deruelle) qui est une 
extension du modele BT convenant aux systemes en champ fort, ou encore ELLl, adapte aux orbites 
quasi-circulaires ( Taylor & Weisbergl|1989 1. lis mettent en jeu les parametres Kepleriens classiques : 



periode orbitale P^, demi grand-axe projete x, excentricite e, longitude du periastre to et temps de 
passage au periastre Tq, ainsi que des parametres « post-kepleriens », pour les theories le permettant. 

III.2.3 Ephemeride de pulsar et signatures de residus 

On a vu dans le paragraphe |III. 1 | que Ton peut ecrire la frequence de rotation / sous la forme d'un 
developpement en serie de Taylor tronque, autour d'une date de reference Tq, avec (/o, /o, /o, • • •) la 
frequence de rotation et ses derivees a cette date {cf. equation III.1| ). Construire une ephemeride pour 



un pulsar consiste a rechercher des parametres de rotation (/o, /o, /o, • • ■), des parametres astrome- 
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triques (ascension droite a, declinaison 8, mouvement propre jia = 6c cos(5) et /ig = 5, parallaxe 
n = l/Jkpc)^ et des parametres orbitaux le cas echeant, permettant de reproduire au mieux la rotation 
du pulsar observee au telescope. 

Pour y parvenir, on se donne un lot de TOAs enregistres a I'observatoire, tj, que Ton convertit 



en temps barycentriques, 7], d'apres III.4 Pour chaque 7], on calcule la phase ^{Tj) predite par le 



modele que Ton souhaite ameliorer, a I'aide de I'equation III.2 Une ephemeride optimale minimise 
les « residus » de chronometrie, donnes par : 

R=pxy ^ , (III.7) 

Dans cette expression, o{Ti) designe le poids attribue a 7] (en effet, I'incertitude sur la mesure des 
TOAs etant variable, on peut choisir des ponderations differentes) et n, est I'entier le plus proche de 

En pratique, on debute 1' analyse des temps d'arrivee avec une ephemeride minimale, donnant une 
connaissance approximative de la phase en fonction du temps. De maniere iterative, on optimise les 
parametres du modele un a un, en fonction de la forme des residus de datation. Plus precisement, 
un parametre mal adapte donne lieu a une signature caracteristique dans le graphe des residus. La 
Figure III.3| illustre quelques exemples de signatures, pour le pulsar J0030+0451, observe avec le 



radiotelescope de Nan9ay sur quatre ans. Dans I'exemple (A), la frequence de rotation est faussee de 
10^^ Hz. A cause de cette erreur sur la frequence du pulsar, I'ecart entre phases predites et observees 
grandit lineairement dans le temps. Dans ce graphe on constate un saut d'une rotation (~ 4,9 ms) 
vers MJD 53900, signifiant que le modele de chronometrie differe d'un tour complet par rapport a la 
rotation reelle. La Figure (B) correspond a une erreur sur /o de 5 x 10^^^ Hz^, induisant un ecartement 
en parabole. Les cas (C) et (D) correspondent a des erreurs sur la position du pulsar. Pour le premier 
cas, la position du pulsar dans le modele de chronometrie est erronee de 1 arcseconde. La signature est 
une sinusoide annuelle, due a I'orbite terrestre. Pour les pulsars ayant un mouvement propre, comme 
PSR J0030+0451, on peut observer des residus tels que ceux montres dans I'exemple (D) : ici on n'a 
utilise aucun mouvement propre dans 1' ephemeride. La signature correspondante est une sinusoide 
annuelle dont 1' amplitude augmente lineairement avec le temps. 

Ayant ameliore le modele de fa9on iterative, on parvient a une ephemeride optimale, et un graphe 



de residus semblable a celui montre en Figure IIL4 Ici les residus sont repartis aleatoirement autour 



de 0, et la deviation standard, ici de 4,3 /is, est du meme ordre de grandeur que la precision moyenne 
sur la mesure des TOAs individuels. L' ephemeride correspondante est donnee dans la Table |IILT 



Les residus de datation accessibles pour les pulsars milliseconde sont generalement depourvus 
de bruit, a I'instar de ce qui est presente en Figure [IIL4 pour PSR J0030+0451. La plupart des pul- 



sars milliseconde sont extremement stables et ne produisent pas de glitches ou de timing noise, au 
contraire des pulsars jeunes, beaucoup plus bruites. Les glitches sont des augmentations soudaines 
de la frequence de rotation, probablement provoques par une modification de la structure de I'etoile 
a neutrons. Un quart environ des pulsars suivis dans le cadre de la campagne de chronometrie (cf. 



IL2.3[ ) sont connus pour avoir produit des glitches ( |Melatos et a/.||2008^ . Le timing noise est la va 



riation aleatoire de la frequence de rotation ou d'une de ses derivees, avec une pseudo-periodicite 
comprise entre quelques jours a quelques mois voire annees. Des illustrations de timing noise sont 



donnees dans Hobbs et al. (2006). Ces deux comportements perturbent la precision des modeles de 



datation. Pour les glitches, il faut incrementer la frequence pr6-glitch d'une certaine quantite pour 
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Figure 111.3 - Exemples de signatures de residus caracteristiques, pour le pulsar J0030+045 1 observe a 
Nan9ay. 



pouvoir suivre la rotation a nouveau. Dans le cas du timing noise, on peut lisser les residus grace a un 
grand nombre de derivees de la frequence. Une autre possibilite, offerte par le logiciel TEMP02, est 
d'ajuster la marche aleatoire due au timing noise par un developpement en sinusoidea^ 



III.3 Campagne de datation 



III.3.1 Couverture des pulsars 

Comme on I'a deja evoque, une campagne de chronometrie de pulsars prometteurs a ete mise 



en place dans les mois precedant le lancement de Fermi ( jSmith et al. 2008 1. Des radiotelescopes 
parmi les principaux dans le monde ont participe : le telescope de 64 m de Parkes, en Australie 



([Manchester et a/.]|2001| ), le telescope Lovell de 76 m a Jodrell B ank, en Angleterre ( |Hobbs et al. 
2004), le radiotelescope de 110 m de Green Bank aux Etats-Unis (Kaplan et al. 2005 1, et la station 



^cf. fonctionnalite FITWAVES de TEMP02 
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Figure 111.4 - Residus de chronometrie optimaux, pour le pulsar J0030+0451 observe au radiotelescope de 
Nan9ay. La deviation standard de ces residus est de 4,3 /is, I'incertitude moyenne sur les TOAs individuels 



etant de 3,5 /is. L'ephemeride correspondante est donnee en Table III.l 



de Nan^ay en France (Theureau et al. 2005| ), telescope meridien equivalent a une antenne de 94 
m {cf. Figure III.5). En outre, des contributions ponctuelles ont ete apportees par les telescopes de 
Westerbork aux Pays-Bas ( |Voute et a/.|2002[ ) et Arecibo, a Porto Rico ( |Dowd etal.\1000\ . 

Les radiotelescopes cites precedemment permettent de couvrir I'essentiel des pulsars connus. Pour 
certains pulsars ne produisant pas ou peu d'emission detectable dans le domaine radio, mais detec- 
tables en rayons X, la chronometrie a ete realisee a I'aide des satellites RXT^et XMlvj^ 



^http ://heasarc. gsfc.nasa.gov/docs/xte/XTE.html 
^http ://xmm.esac.esa.int/ 



Norn J0030+0451 

Intervalle de temps (MJD) 53285,9—54829,7 

Nombre de TOAs 413 

RMS des residus (ils) 4,26 

Ascension droite, a 00:30:27,4301(11) 

Declinaison, S +04:5 1 :39,78(4) 

Frequence, /o (s"') 205,530699274958(19) 

Premiere derivee de la frequence, /o (s^^) — 4,2988(7)x 10^'* 

Mouvement propre en ascension droite. Ha (masyr^') —3,5(18) 

Mouvement propre en declinaison, (masyr^') .... —7(5) 

Parallaxe, TZ (mas) 3,4(5) 

Date de reference de la frequence (MJD) 50984,4 

Date de reference de la position (MJD) 50984,4 

Mesure de dispersion, DM (cm^'pc) 4,33822 

Ephemeride planetaire DE200 

Unites TDB 



Tableau III.l - Exemple d' ephemeride de pulsar, ici pour PSR J0030+0451. Les nombres donnes entre 
parentheses sont les incertitudes aid sur les dernieres decimales citees. 
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Figure III. 5 - Le radiotelescope de Nangay (region Centre) vu du ciel. Le miroir visible en arriere-plan 
(200 m X 40 m) joue le role de premier reflecteur. C'est un miroir plan et orientable, au contraire du second 
reflecteur au premier plan (300 m x 35 m), fixe et de forme spherique. Les miroirs sont separes de 460 m. 
Ce second reflecteur concentre les ondes incidentes vers le point focal, au centre. Le point vise dans le ciel 
se deplace au cours du temps a cause de la rotation de la Terre, et par consequent le point focal se deplace. 
Pour cette raison, les antennes sont disposees dans un « chariot focal » (on peut le distinguer au centre 
de I'image) pouvant se deplacer sur des rails. Grace a ce mouvement, le radiotelescope de Nan9ay peut 
couvrir a tout moment ±5 en ascension droite. Le chariot focal embarque deux recepteurs observant dans 
les gammes de frequence 1,1 a 1,8 GHz et 1,7 a 3,5 GHz, permettant de couvrir la raie de I'hydrogene HI 
et au-dela. L' ensemble forme un telescope meridien, equivalent a une antenne simple de 94 m de diametre. 



Ces instruments ont couvert prioritairement les 218 pulsars consideres comme prometteurs. Ce 
suivi des pulsars candidats a la detection par Fermi est en fait similaire aux travaux realises auparavant, 
dans le contexte de la mission CGRO ( [Arzoumanian et aL\\l994\ [Johnston et a/.|p^95| [D'Amico 



et al. 1996) a la difference que Fermi fonctionne principalement en mode de balayage, et non de 



pointe. Cela signifie que si pour CGRO, la duree et la direction des pointes definissaient la liste des 
pulsars a observer et la frequence des observations necessaire dans un temps limite, un suivi continu 
et homogene est requis pour Fermi qui balaye le ciel en permanence. Une autre difference est bien 
evidemment la difference de performance entre EGRET et le LAT. Le LAT, plus sensible, a acces a 
davantage de sources. 

Les periodes d' observation prolongees favorisent 1' apparition de timing noise dans les residus de 
chronometrie des pulsars jeunes et bruites, generalement a grande valeur de E et done bons candidats 
a la detection en rayons 7. La consequence est que la bonne connaissance de la phase en fonction du 
temps requiert plus de complexity dans le modele de chronometrie. Une illustration de ce besoin de 
complexity est le cas du pulsar du Crabe, quotidiennement observe a I'aide d'un des instruments de la 
station de Jodrell Bank pres de Manchester. Avec ces observations des ephemerides sont construites 
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chaque mois, de telle sorte que leur periode de validite ne couvre que la duree du mois en question^ 
Sauf exceptions, un simple developpement de la frequence de rotation en fonction du temps en /o, 
/o et /o suffit a assurer une precision de I'ephemeride de quelques dizaines de /is. Lorsque la plage 
d' observation du pulsar du Crabe se prolonge, le timing noise pent devenir preponderant et imposer 
r utilisation de termes d'ordre superieur sur la frequence de rotation, voire un developpement en 
sinusoides. La Figure III.6 (A) montre les residus obtenus pour 697 TOAs enregistres conjointement 
a I'aide des radiotelescopes de Jodrell Bank et de Nan^ay, entre le 20 juin 2008 et le 9 avril 2009, 
avec un modele de rotation restreint a /o, /o et /q. Le timing noise domine, 1' excursion en residus 



s'etendant sur 2 ms environ, soit plus de 5% de la periode de rotation du pulsar. Dans la Figure IIL6 



(B) on a lisse ce bruit oscillant par un developpement en harmoniques. Cette fois I'amplitude typique 
des residus de chronometrie est de 50 /is, de I'ordre du millieme de rotation du pulsar. 



La Figure III.7 dresse un bilan de la couverture des pulsars dans le cadre de la mission Fermi, 
jusqu'a mai 2009. Les pulsars definis a priori comme prioritaires sont indiques par des points epais. 
Rappelons que le seuil a E = 10^^ erg/s qui avait ete fixe ne tenait pas compte de I'effet Shklovskii, 
ce qui explique que certains d' entre eux se trouvent en de^a de la ligne pointillee. Les carres vides in- 
diquent les pulsars qui n'etaient pas consideres comme prioritaires, mais pour lesquels des parametres 
de rotation ont ete fournis au moins ponctuellement. Le nombre de ces pulsars est de 500 environ. Les 
objets non couverts pendant la premiere annee de campagne sont indiques par des points fins. 

La centralisation des ephemerides issues de cette campagne de chronometrie dans une base de 
donnees constitue I'un des travaux effectues au cours de cette these. Afin de rendre accessible les 
parametres de rotation aux pulsaristes y de la collaboration Fermi-LAT, nous avons construit une 
base de donnees php/MySQL accessible par internej^ Le mode de fonctionnement est le suivant : en 
amont, le site internet absorbe et conserve les donnees provenant des radiotelescopes et telescopes 
X, a savoir les ephemerides et eventuellement les graphes de residus et les profils de reference. En 
aval, la base de donnees permet a I'utilisateur cote LAT d'effectuer des requetes : pulsar etudie, 
gamme de validite de I'ephemeride souhaitee, observatoire souhaitee, etc. Puis, une fois la requete 
passee, le site propose a I'utilisateur de telecharger I'ephemeride voulue au format initial (propre aux 
telescopes radio/X) ou dans un format special adapte aux outils de la collaboration Fermi. Sur une 
base mensuelle a bimestrielle, cette base de donnees a ete mise a jour de fa9on a disposer d'un suivi 
continu des pulsars. In fine, les ephemerides utilisees dans des publications de Fermi sont mises a 
disposition libre sur un serveur du Fermi Science Support Center (FSSCj^ 

Enfin, soulignons que par le biais de la campagne de chronometrie, presque tous les pulsars mil- 
liseconde du champ galactique ont ete suivis. Des ephemerides contemporaines avec les observations 
de Fermi n'ont pas ete construites pour six des 72 objets de la liste |A.1| seulement ; a savoir PSR 
J1125-6014, J1216-6410, J1618-39, J1629-6902, J1757-5322 et J1933-6211. Cependant, ces 
six pulsars ont des valeurs de E intermediaires et sont situes a des distances relativement grandes 
parmi les MSPs connus. Disposant de chronometrie pour I'essentiel des 72 pulsars milliseconde ga- 
lactiques, nous avons eu 1' opportunity de rechercher des pulsations parmi un echantillon de pulsars 
limite non pas par une hypothese formulee a priori sur les meilleurs candidats a remission, mais par 
la completude des recherches de MSPs conduites dans le domaine radio. 



''cf. http ;//www.jb.man.ac.uk/ pulsar/crab.html 

^http ://www.cenbg.m2p3.fr/ephem/ 

'http ://ferim.gsfc. nasa.gov/ssc/data/access/lat/ephems/ 
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Figure 111.6 - Residus de chronometrie pour le pulsar du Crabe, observe a 697 reprises entre le 20 juin 
2008 et le 9 avril 2009 par les radiotelescopes de Nan9ay et de Jodrell Bank. (A) : seule la frequence de 
rotation du pulsar et ses deux premieres derivees ont ete ajustees. (B) : le modele de datation contient en 
plus un developpement en harmoniques de la phase rotationnelle. 



III.3.2 Analyse des donnees de Nangay 

En plus de la centralisation des ephemerides dans une base de donnees accessible par internet, la 
construction d' ephemerides a partir des donnees radio du telescope de Nan9ay a ete un autre volet de 
la participation a la campagne de chronometrie au cours de cette these. 

Le radiotelescope de Nan§ay a servi a la construction d'ephemerides pour environ 150 pulsars. 
La procedure de mise a jour des ephemerides pour ce grand nombre de pulsars a ete rendue automa- 



tique. La Figure IILS montre les differentes etapes de la procedure automatique d'optimisation des 
parametres de rotation. On dispose initialement de I'ephemeride valide au moment de la mise a jour 
precedente. Si de nouvelles observations ont ete realisees pour le pulsar considere, les profils sont 
alignes et additionnes. A partir des observations, nouvelles et anciennes, et utilisant comme reference 
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Figure III.7 - Couveiture des pulsars dans le cadre de la campagne de datation. Les points epais sur ce 
diagramme P -P sont les pulsars a grand E, definis avant le lancement comme bons candidats a remission 
7. Les Carres vides indiquent les pulsars a valeur de E moindre, pour lesquels des ephemerides ont tout de 
meme ete fournies au cours de la premiere annee. Les points fins sont les pulsars non couverts. 



le profil integre, on calcule les temps d'arrivee. Pour ces TOAs, on calcule les residus de chronometrie 
« pre-fit », correspondant a I'ephemeride initiale, ainsi que les residus « post-fit » pour lesquels un 
ou plusieurs parametres de rephemeride ont ete ajustes. Si les residus « pre-fit » sont meilleurs que 
les residus « post-fit », il n'est pas necessaire de mettre a jour I'ephemeride : les parametres initiaux 
ajustent toujours les donnees de fa9on satisfaisante. Dans le cas contraire, la procedure est repetee 
jusqu'a la convergence des residus de chronometrie. 



III.4 Verification des outils 

Pour r analyse des donnees LAT, une distribution d'outils dedies a ete mise en place, incluant 
des codes consacres au calcul des phases des pulsars, a partir des dates enregistrees au telescope. 
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Figure 111.8 - Processus de construction des ephemerides a partir des observations de pulsars faites a 
Nan^ay. 



Cette distribution, appelee Science Tools, est developpee, maintenue et distribuee par le FSSC 



10 



avec I'aide des membres de la collaboration Fermi. Outre la construction d' ephemerides couvrant 
les observations 7 par Fermi, les telescopes radio et X ont egalement offert la possibilite de tester 
les outils d' analyse de Fermi a I'aide de donnees exterieures, en particulier gtbary, procedure de 
barycentrisation des dates des photons 7, et gtpphase, outils de calcul des phases des pulsars a partir 
des dates barycentrees. Ce travail de verification est detaille dans Smith et ar] ( 2008| ). Ici nous avons 
retenu les deux tests les plus contraignants, realises avec des pulsations geantes du pulsar milliseconde 
B 1937+21 et des TOAs du pulsar milliseconde en orbite binaire J0437— 4715. 



cf. http ://fermi.gsfc.nasa.gov/ssc/data/analysis/software/ 



64 



in.4. VERIFICATION DES OUTILS 



III.4.1 Pulsations geantes de B1937+21 



Le phenomene de pulsations radio geantes (Giant Radio Pulses, GRPs) se manifeste pour quelques 



pulsars, jeunes ou anciens, tels que le pulsar du Crabe ou certains pulsars milliseconde (Johnston & 



Romani 2004, Knight et al. 2006b), et consiste en de breves bouffees d'emission radio, beaucoup 
plus intenses que le niveau d'emission moyen. Le pulsar milliseconde isole B1937+21 (P = 1,56 
ms) est un emetteur de GRPs ( Cognard et a/.|[T996] ). L' emission de GRPs pour ce pulsar consiste en 
des composantes principale et secondaire, decalees de 55 a 70 /is par rapport aux pics standard et 
intervenant dans des fenetres de temps inferieures a 10 /is. 

Entre le 13 mars et le 7 avril 2007, 251 GRPs pour lesquels le rapport signal- sur-bruit etait supe- 
rieur a 30 deviations standard ont ete enregistres au radiotelescope de Nan9ay. Une ephemeride a ete 
produite a I'aide des memes observations. Pour verifier les outils de Fermi, les GRPs ainsi enregistres 
ont ete traites de la maniere suivante : 

- Conversion des TOAs, donnes dans un format de temps specifique a I'observatoire, vers des 
dates MET (Mission Elapsed Time, format des dates de Fermi ainsi que de certaines missions 
spatiales), en 1' occurence le nombre de secondes depuis le Janvier 2001 a minuit (UTC), 
dans I'echelle de Temps Terrestre (TT). 

- Calcul de la position de I'observatoire dans le systeme solaire en fonction du temps et conver- 
sion vers le format specifique des donnees de Fermi. 

- Transfert des temps a I'observatoire, vers le barycentre du systeme solaire, grace a gtbary. II faut 
cependant corriger des effets dispersifs au prealable si besoin est : gtbary ne permet pas cette 
correction puisqu'elle est negligeable pour les photons ydont la frequence est quasi-infinie. 

- Enfin, calcul de la phase du pulsar pour chaque TOA, et ce grace a I'outil gtpphase. 



En realisant les etapes precedentes, on obtient les phasogrammes montres en Figure III.9 A et B. 
Les composantes principale et secondaire des GRPs retardent de respectivement 60,1 et 67,3 /is par 
rapport aux pics d'emission classique, avec des deviations standard de 1,9 et 2,4 /is. Ce resultat, en 
accord avec ceux de Kinkhabwala & Thorsett] ( |2000[ ), demontre que les outils d' analyse de Fermi, 
de la procedure de barycentrisation au calcul des phases, sont precis a I'echelle de quelques /is. Ce- 
pendant, le pulsar B 1937-1-21 est isole, par consequent la procedure de correction des orbites binaires 
n'est pas testee ici. De plus, le faible intervalle de temps considere (trois semaines environ) ne permet 
pas de deceler des erreurs de positionnement de I'observatoire, eventuellement visibles a plus long 
terme. 



III.4.2 Orbites complexes : cas de J0437-4715 

Le pulsar milliseconde J0437— 4715 (P = 5,76 ms) appartient a un systeme binaire. La periode 
de I'orbite est de 5,7 jours environ, et cette o rbite permet la mesure de parametres post-kepleriens 
(van Straten et al. 2001 , Verbiest et a/.|2008 1. A partir de 1' ephemeride pour ce pulsar foumie dans 



Verbiest et a/.| ( |2008] ), de 300 ns de precision, 37 TOAs au telescope de Parkes ont ete simules sur 
un intervalle de temps de 500 jours (de decembre 2003 a mai 2005), grace a la fonctionnalite FAKE 
de TEMP02. La procedure de test est la meme que celle decrite precedemment, a la difference que 
cette fois la procedure permettant a gtpphase de corriger de I'orbite des pulsars binaires est mise a 



I'epreuve. Les 37 phases calculees a I'aide des outils de Fermi sont montrees dans la Figure IILIO En 
moyenne, les TOAs traites dans ce test ont un retard moyen de 320 ns par rapport aux memes TOAs 
analyses avec TEMP02. On en deduit ici que les outils de barycentrisation et de calcul des phases 
gtbary et gtpphase donnent des resultats compatibles avec ceux de I'outil standard TEMP02, pour 
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(A) GRPs de B 1937+21 sur I'ensemble de la rotation 
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Figure 111.9 - Pulsations radio geantes du pulsar millisecondeB 1937+21 enregistrees aNan9ay et analysees 
avec les outils de Fermi (ligne continue, echelle de gauche). Un profil radio standard du pulsar a 1,4 GHz est 
montre en ligne pointillee (echelle de droite). Le pic de GRPs a une largeur d'environ 2 /is, demontrant la 
precision des outils d' analyse. La figure (A) est issue de [Smith et a/. | ( [20081 ). (^) • zoom sur la composante 
principale des pulsations geantes. 
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Figure III. 10 - Temps d'arrivee simules pour le pulsar binaire J0437— 4715, observe a Parkes, et analyses 
avec les outils de Fermi. La phase absolue est definie de telle sorte que le pic radio principal soit centre 
a 0,5. Les TOAs ont un retard moyen de 0,32 jj.s, indiquant que les outils d' analyse ont une precision de 
I'ordre de la ^us. Figure issue de|Smith et a/. | ([2008]). 



un pulsar milliseconde dans un systeme binaire complexe et pour des dates couvrant un intervalle de 
temps de plus d'un an d' observation. 

D'autres tests ont ete realises au cours de cette these, utilisant notamment des GRPs du pulsar 
du Crabe enregistres a Nan9ay, et des observations en rayons X de PSR J02 18+4232 avec XMM. 
Ces autres verifications ont amene a des conclusions identiques. Par consequent, les outils d' analyse 
temporelle distribues par le FSSC foumissent des resultats compatibles avec ceux des outils deja 
existant, a I'echelle de la /is. 
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Chapitre IV 
Le Large Area Telescope a bord de Fermi 
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DANS ce chapitre nous presentons I'lnstrument utilise en majorite au cours de cette these, le Large 
Area Telescope (LAT) a bord du satellite Fermi, lance en juin 2008. Dans un premier temps 
nous faisons un rappel des proprietes de Fermi. Celui-ci embarque deux instruments : le GBM dedie 
a la detection de sursauts 7 a des energies comprises entre 8 keV et 40 MeV, et le LAT, instrument 
principal embarque par l'observatoire, sensible a des photons d'energie entre 20 MeV et 300 GeV. 
Le LAT est lui-meme compose de differents elements dont nous rappelons le fonctionnement dans 
la deuxieme partie de ce chapitre. Ceux-ci sont le bouclier anti-coincidences, element participant 
au rejet des rayons cosmiques, le trajectographe, determinant la direction des photons detectes, et 
le calorimetre qui mesure I'energie qu'ils deposent. Ces differents sous-systemes interagissent pour 
rejeter les particules chargees indesirables, et pour reconstruire la direction et I'energie des photons 
tombant dans le telescope. Nous nous penchons enfin sur les performances du LAT, resolution 
angulaire, resolution en energie ; mais egalement la precision en temps, vitale pour I'etude des 
pulsars milliseconde. 



Le Fermi Gamma-ray Space Telescope (FGST), a I'origine nomme GLAST (pour Gamma-Ray 
Large Area Telescope), est un observatoire dedie a la detection de sources de rayonnement 7, mis en 
orbite autour de la Terre le 11 juin 2008. La mission Fermi succede a CGRO, en activite pendant 



IV.l L'observatoire Fermi 
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Figure IV.l - Le Fermi Gamma-Ray Space Telescope, surmonte du LAT, visible dans la partie superieure 
de I'image. 



les annees 90. A I'image de CGRO qui embarquait quatre instruments, I'observatoire spatial Fermi 
dispose a son bord du GBM (Glast Burst Monitor) ( [Meegan et al. 2007 ), instrument consacre a la 
detection de sursauts yd'energie comprise entre 8 keV et 40 MeV, et du LAT (Large Area Telescope, 
cf. Figure |IV.1[), instrument principal de Fermi dont nous allons rappeler le fonctionnement dans la 



suite de ce chapitre. Le lecteur interesse pourra trouver davantage de details sur le LAT dans Atwood 



et al. (2009 1. Fermi et ses instruments devraient fonctionner de cinq a dix ans. 



Le satellite est place sur une orbite basse, quasi-circulaire, a environ 565 km d'altitude et 25,6° 
d'inclinaison, pour une periode orbitale de \ j ^lf — 95 minutes, avec a = Rt+ 565 km. II fonctionne 



essentiellement en mode de balayage : I'axe du satellite pointe dans une direction a 35° du zenith vers 
le pole de 1' orbite, pendant une orbite. Au cours de la rotation suivante I'axe du satellite pointe a -35° 
du zenith, et ainsi de suite. Le ciel est balaye de fagon quasi-uniforme toutes les deux orbites (soit 
environ trois heures) grace au grand champ de vue du LAT (2,4 sr). A plusieurs reprises au cours des 
60 jours qui ont suivi le lancement Fermi a opere en pointe, vers le pulsar de Vela notamment. Dans 
ce mode le satellite pointe son axe vers la cible desiree et bascule vers une autre direction lorsque la 
Terre s'approche de la ligne de visee. 
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Figure IV.2 - Schema en coupe du LAT, montrant un photon 7 tombant dans le trajectographe et in- 
teragissant avec la matiere pour former une paire electron / positron. Le calorimetre, dans la partie in- 
ferieure, absorbe les particules chargees et mesure leur energie. L' ensemble est protege par le bouclier 
anti-coincidences. 



IV.2 Le Large Area Telescope 

IV.2.1 Principe de fonctionnement 

Le Large Area Telescope (LAT) est un telescope sensible au rayonnement yd' energie comprise 
entre 20 MeV et 300 GeV. Ses trois elements principaux sont le trajectographe, le calorimetre et 
le bouclier anti-coincidences, auxquels il faut ajouter I'electronique associee, enregistrant les infor- 
mations. La technique de detection du rayonnement 7 est la suivante : un photon 7 tombe dans le 
telescope et interagit avec la matiere pour se convertir en une paire electron et positron. Ces particules 
chargees sont detectables par les elements du telescope, qui, en determinant leur trajectoire et leur 
energie, permettent de remonter a la direction et I'energie totale du photon 7 initial. Ce principe de 



fonctionnement est illustre par la Figure IV.2 



Les photons 7 se convertissent en paires en presence de matiere, par consequent 1' atmosphere 
terrestre empeche I'observation au sol dans cette gamme d'energie, d'ou la necessite de satelliser 
le LAT. Place en orbite autour de la Terre, Fermi est cependant soumis a un intense flux de parti- 
cules chargees qui interagissent egalement avec les elements du LAT. Ces interactions peuvent etre 
interpretees a tort comme des detections de photons 7. Pour cette raison, le detecteur est entoure du 
bouclier anti-coincidences charge de rejeter le plus possible de ces particules chargees indesirables. 
Nous reviendrons par la suite sur le rejet du fond. 



IV.2.2 Le bouclier anti-coincidences 

Le but principal du bouclier anti-coincidences (ACD) est le rejet des rayons cosmiques, particules 
chargees provenant de toutes les directions de I'espace. Comme on I'a evoque precedemment elles 
peuvent interagir avec le detecteur, et engendrer de fausses detections de photons. L'objectif fixe pour 
le rejet des particules non 7est de 99,97%. L'ACD est compose de 89 scintillateurs plastiques d'un cm 
d'epaisseur, dits tuiles, et de huit rubans couvrant les interstices. L'ACD recouvre le trajectographe et 
le calorimetre. 

Le principe de fonctionnement est le suivant : les rayons 7 interagissent peu avec I'ACD. On 
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s'attend alors a ce que la detection d'un rayon 7 par le trajectographe coincide avec une absence de 
detection de signal par I'ACD (c'est-a-dire, une « anti-coincidence »). Au contraire, une particule 
chargee tombant sur le telescope depose une partie de son energie dans I'ACD, puis interagit avec le 
trajectographe. Ainsi on s'attend a ce que pour les particules chargees, la detection de signal par les 
deux instruments coincide. 

La nouveaute du systeme ACD embarque par le LAT est qu'il ne s'agit pas d'une structure mo- 
nobloc mais de nombreuses tuiles et rubans, permettant cette fois de connaitre les positions 011 des 
particules ont ete detectees. Pour les rayons yde haute energie, typiquement au-dela de quelques GeV, 
la gerbe engendree par I'interaction dans les instruments peut contenir des photons sortants, diriges 
vers le bouclier anti-coincidences. Ce phenomene est appele back-splash. La probabilite que ces pho- 
tons d' energie moindre (typiquement 100 keV a 1 MeV) interagissent avec I'ACD est faible. S'il y a 
neanmoins interaction, un ACD monobloc sans informations spatiales n'a aucun moyen de distinguer 
ce back-splash de rayonnement cosmique traditionnel. En consequence I'ACD produit systematique- 
ment un signal de self-veto qui empeche la detection du photon initial. Cet effet a rendu EGRET tres 
peu sensible aux photons au-dessus de 10 GeV. L'ACD du LAT est segmente, permettant I'acces a 
I'information spatiale en cas de signal, et reduisant ainsi I'effet de self-veto. II a ete con9u pour que 
seuls 20% des photons de 300 GeV tombant sur le LAT ne soient rejetes en raison de coincidences 
fortuites. 

IV.2.3 Le trajectographe 

Cet element (egalement appele Tracker, TKR) est con9u pour determiner la direction du rayon 7 
tombant dans le telescope. Le principe est ici d' avoir un materiau riche en protons, pour que la pro- 
babilite que des photons / incidents s'y convertissent en paire electron / positron soit maximale. Puis, 
en cas de conversion, de mesurer les points de passage de ces particules chargees dans le detecteur, 
et en joignant les points, en deduire la position oii la conversion a eu lieu (le « vertex ») et Tangle 
d'incidence de ce photon. 

Le TKR est un assemblage de 4 x 4 tours, chacune consistant en un empilement de 19 plateaux. 
Chaque plateau contient deux plans de detecteurs a pistes de silicium (Silicon-Strip Detectors, SSDs). 
disposes aux deux faces et dont les pistes sont orientees dans la meme direction. Les plateaux aux 
extremites ne comportent qu'un seul plan de SSDs. Chaque plateau est oriente perpendiculairement 
a celui qui le precede. De cette maniere, 1' ensemble se repartit en deux orientations distinctes, x et y. 
Deux plans successifs d'orientations perpendiculaires forment des doublets xy separes de seulement 
3 mm. Ainsi, on peut localiser en trois dimensions les points de passage des particules chargees en 
reperant les doublets xy qui ont detecte un signal, et leur position sur I'axe z, orthogonal au plan xy. 

Afin de maximiser la conversion de photons 7 en paires electron / positron, des couches de feuilles 
de tungstene (Z = 74) d'epaisseur plus ou moins importante ont ete intercalees dans les differents pla- 
teaux. Le schema de repartition du tungstene dans le detecteur repond au besoin de disposer a la fois 
d'une bonne reconstruction spatiale et d'une surface effective de collection des photons maximale. La 
precision de la reconstruction spatiale se heurte principalement a la diffusion multiple des particules 
de basse energie. A cause de cet effet, les faibles epaisseurs de materiau sont privilegiees. Cependant, 
une grande epaisseur de tungstene apporte plus de longueurs de radiatiorQ accroissant ainsi la capa- 
cite de conversion du materiau pour une surface donnee. La conception adoptee pour le trajectographe 

'La definition de la longueur de radiation Xq est la suivante : un electron traversant une epaisseur Xq d'un materiau 
donne voit son energie diminuer d'un facteur e par pertes radiatives (par Bremsstrahlung essentiellement). 
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est la suivante : dans les 12 premiers plateaux, les couches de tungstene sont fines representent 0,03 
longueurs de radiation. Pour les quatre plateaux qui suivent, elles sont epaisses et equivalent a 0,18 
longueurs de radiation. Enfin, les trois derniers plateaux ne contiennent pas de tungstene. 

IV.2.4 Le calorimetre 

Le calorimetre (CAL) a pour but de mesurer I'energie des particules chargees nees de la conversion 
du photon 7 initial dans le trajectographe, pour remonter a I'energie de ce photon. Chacune des 16 
tours du TKR possede un module de calorimetre dans sa partie inferieure. Un module comprend 
huit couches de 12 barreaux d'iodure de cesium (Csl) enrichi au thallium (Tl), chaque barreau ayant 
deux paires de photodiodes situees a ses extremites, mesurant la lumiere produite par scintillation du 
cristal. Grace a sa densite, le calorimetre a une epaisseur equivalente de 8,6 longueurs de radiation. 
A nouveau, les barreaux sont arranges perpendiculairement d'une couche a 1' autre, de fagon a ce que 
Ton puisse mesurer la position (x,y) du depot d'energie dans la couche consideree. La precision sur la 
mesure de la position est de quelques mm pour les faibles depots d'energie (typiquement 10 MeV), a 
moins d'un mm pour E > I GeV. L' information spatiale apportee par le CAL participe a I'algorithme 
de calcul des trajectoires des particules chargees dans le telescope, en particulier pour les photons y 
de haute energie, dont les cascades ont un developpement plus grand dans le detecteur. La resolution 
en energie (AE/E) du calorimetre est inferieure a 20%, voire a 10% dans certaines bandes en energie. 

IV.2.5 Reconstruction des evenements 

Les paragraphes qui precedent presentent les differents elements du LAT comme relativement 
independants. En realite les informations delivrees par les sous-systemes doivent le plus souvent etre 
mises en relation pour determiner la direction d' incidence et I'energie de la particule initiale, creant 
une gerbe electromagnetique dans le detecteur. 

La determination de la trace est centrale dans la reconstruction des evenements : elle foumit des a 
priori sur ce que devraient observer le CAL et 1' ACD. A partir des positions xy du TKR ou un signal 
a ete repere, des algorithmes generent des trajectoires de particules, et il s'agit alors de retenir la plus 
plausible. Le critere utilise le plus souvent fait intervenir le calorimetre : on fait I'hypothese que le 
barycentre du depot d'energie dans le CAL doit etre sur la trajectoire de la particule primaire. A haute 
energie ou la diffusion multiple est la moins importante, la gerbe est etroite : ainsi ce barycentre de 
depot dans le CAL fournit une contrainte de plus en plus forte sur la trajectoire. II se pent cependant 
que peu ou pas d'energie soit deposee dans le calorimetre. Dans ce cas I'algorithme de recherche de 
trajectoire optimale ne peut tenir compte que de la position des doublets xy ayant ete actives dans le 
TKR. 

La trajectoire ainsi determinee permet alors de raffiner la mesure d'energie par le calorimetre. 
L' estimation la plus rudimentaire se fait en effet en sommant les energies mesurees par tous ses 
cristaux. Selon le profil de la gerbe, on peut estimer I'energie perdue par echappement des particules 
sur les cotes, a I'arriere, ou encore dans les interstices separant les composants du calorimetre. A 
basse energie, une part importante de I'energie est laissee dans le trajectographe, et par consequent il 
est important de mesurer ce depot et de I'ajouter a la quantite mesuree par le calorimetre. Le nombre 
de pistes de silicium ayant detecte un signal foumit une estimation de I'energie ainsi deposee dans le 
TKR. Dans tous les cas, la forme de la cascade se developpant dans le telescope est ajustee de fagon 
a bien estimer I'energie. 
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IV.2.6 Declenchement et rejet du fond 

Le LAT est essentiellement traverse par des rayons cosmiques, les photons 7 ne representant 
qu'une part marginale des particules incidentes. Le rejet des particules chargees indesirables est done 
indispensable. Cependant, on ne peut pas simplement rejeter les evenements pour lesquels du signal 
est detecte dans I'ACD : le backsplash engendre par les photons 7 est souvent recupere par des tuiles 
de I'ACD. Les tuiles ou la detection de signal est autorisee sont voisines de la trajectoire hypothetique 
de la particule primaire, et I'energie reconstruite permet de fixer la taille de la region de I'ACD auto- 
risee. Plus la particule primaire est energetique, plus la gerbe est etroite, moins grande est la surface 
de I'ACD ou du signal doit etre detecte. Si des tuiles en dehors de la region permise sont activees, 
I'evenement est rejete. Ceci est un exemple des verifications faites pour le rejet du fond, faisant inter- 
venir la encore les differents sous-systemes pour declencher 1' acquisition des donnees. Sur les 3 kHz 
d'evenements consideres par le LAT en moyenne, seuls 400 Hz environ sont enregistres et transmis 
au sol. 

Les evenements transmis sont ensuite analyses, afin de determiner la qualite de la reconstruction : 
precision de la direction d' incidence, de la mesure de I'energie de la particule primaire, et probabilite 
que la particule corresponde a un photon. A ce stade, les algorithmes mettent en oeuvre des arbres 
decisionnels (Classification Trees, CT), entraines par des simulations de la reponse attendue du LAT. 
Les evenements attribues a des photons sont distribues en classes, selon la qualite de leur reconstruc- 
tion : Transient, Source et Dijfuse, constituant une nouvelle couche de rejet du fond et optimisees pour 
differents objectifs scientifique^ 

La classe Transient forme le jeu de coupure le plus lache, donnant le plus de statistique, au prix 
d'une qualite de reconstruction moindre et d'une contamination par des evenements de fond impor- 
tante. Le taux d'evenements Transient est d'environ 5 Hz dont 2 Hz de contamination par du fond. Les 
evenements de type Source forment une sous-classe des evenements Transient, de meilleure qualite 
et correspondant a des photons 7 avec une plus grande probabilite. Pour cette classe le taux d'evene- 
ments total est d'environ 1 Hz, avec un fond estime a environ 0,4 Hz. Enfin, les evenements de classe 
Diffuse, une sous-categorie de la classe Source, sont les mieux reconstruits et ont la plus grande pro- 
babilite de correspondre a de veritables photons, avec un taux de contamination de 0,1 Hz environ sur 
0,5 Hz au total. En pratique, c'est cette derniere classe qui est utilisee pour I'etude des pulsars avec 
le LAT, le rapport signal-sur-bruit y etant le meilleur. Les evenements n'appartenant a aucune de ces 
trois classes se repartissent en electrons, positrons, ions lourds, et dechets. 



IV.2.7 Performances 



Les performances du LAT et certaines de ses caracteristiques sont donnees dans la Table IV. 1 Les 
proprietes du LAT sont generalement determinees par les performances individuelles des differents 
instruments ainsi que par les algorithmes d' identification des rayons 7, a la fois dans la reconstruction 
de leur energie et de leur direction d'incidence, mais egalement dans le rejet des detections fortuites. 

L' observation de pulsars en rayons 7 avec le LAT regroupe trois aspects : spatial, spectral et 
temporel. Les performances de datation du LAT sont mentionnees dans la suite de ce chapitre. Le bon 
fonctionnement des horloges embarquees par Fermi est primordial pour 1' analyse des pulsars dans la 
mesure ou la dimension temporelle est ce qui fait leur specificite. 



^Avec r amelioration des algorithmes decisionnels, de nouvelles classes pourraient etre definies a I'avenir pour 1' ana- 
lyse des donnees du LAT. 
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Parametre Valeur 

Dimensions 1,8 x 1,8 x 0,72 m^ 

Masse totale 2789 kg 

Gamme d' energie 20 Me V - 300 Ge V 

Surface effective maximale a incidence normale 9500 cm^ 

Champ de vue 2,4 sr 

Temps mort 26,5 ^s 

Precision en temps < 1 ;tls 

Resolution en energie (incidence normale) : 

100 MeV - 1 GeV 9% - 15% 

1 GeV - 10 GeV 8% - 9% 

> 10 GeV 8,5% - 18% 

Resolution angulaire ( incidence normale, 68% de la PSF) : 

100 MeV 3,5° 

1 GeV 0,6° 

> 10 GeV < 0,15° 

Tableau IV.l - Caracteristiques et performances du LAT. 



A chaque classe d'evenement, Transient, Source ou Diffuse, correspond une « fonction de reponse 
instrumentale » (Instrument Response Function, IRF), que Ton peut formuler de la maniere suivantej^: 



IRF(0,£) = Aeff(0,£) X PSF(0,£) X AE{d,E) (IV.l) 

Dans cette expression, et E designent Tangle d'incidence par rapport a I'axe du telescope et 
I'energie reconstruits. Aeff(0,£') est la surface effective de collection. EUe est la plus importante dans 



la classe Transient, la plus faible pour la classe Diffuse. EUe est illustree en Figure pAOj pour la classe 
Source. 



La Point Spread Function, decrite par le terme PSF(0,£') dans I'equation IV.l represente la 
reponse spatiale du systeme pour une source ponctuelle. La resolution angulaire concerne tous les 
aspects de 1' analyse des donnees du LAT. EUe resulte essentiellement des performances du trajecto- 



graphe. Son comportement en fonction de I'energie des photons incidents est montre en Figure IV.4 



pour la classe Source. Comme deja mentionne, la reconstruction des directions est meilleure a haute 



energie ; en effet la diffusion multiple est importante a basse energie. La Figure IV.4 indique que 68% 
des photons de 100 MeV reellement emis par la source ponctuelle consideree sont detectes avec une 
separation angulaire inferieure ou egale a 3,5° (en incidence normale). A 1 GeV, la distance angulaire 
correspondante est de 0,6°, et continue de decroitre lorsque I'energie augmente. 



La Point Spread Function, decrite par le terme PSF(0,£') dans I'equation IV.l represente la 
reponse spatiale du systeme pour une source ponctuelle. La resolution angulaire concerne tous les 
aspects de 1' analyse des donnees du LAT. EUe resulte essentiellement des performances du trajecto- 



graphe. Son comportement en fonction de I'energie des photons incidents est montre en Figure IV.4 



pour la classe Source. Comme deja mentionne, la reconstruction des directions est meilleure a haute 



energie ; en effet la diffusion multiple est importante a basse energie. La Figure IV.4 indique que 68% 



^Notons que cette formulation est elementaire et ne rend pas compte de la topologie de I'instrument. Par exemple, 
I'energie et la trace des particules interagissant au niveau des interstices separant les tours du TKR sont moins bien 
reconstruites, ce qui introduit une dependance en (x,y) de la qualite de reconstruction. 
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Figure IV.3 - Surface effective de collection en fonction de I'energie, pour des photons en incidence 
normale (trait plein), et a 60° d'inclinaison (ligne pointillee), dans la classe d'evenements Source. Figure 
extraite de |Atwood et a/.| ( |2009l ). 



des photons de 100 MeV reellement emis par la source ponctuelle consideree sont detectes avec une 
separation angulaire inferieure ou egale a 3,5° (en incidence normale). A 1 GeV, la distance angulaire 
correspondante est de 0,6°, et continue de decroitre lorsque I'energie augmente. 

En resume, les performances du LAT sont tres superieures a celles de son predecesseur, EGRET. 
II est moins assujetti au back-splash, ainsi sa gamme d'energie s'etend de 20 MeV a 300 GeV, la ou 
EGRET n'observait que jusqu'a quelques dizaines de GeV. Le champ de vue est de 2,4 sr contre 0,4 
pour EGRET, et la surface effective maximale, qui est de 9500 cm^, est six fois superieure environ. 
Entre les deux telescopes le temps mort de I'electronique a ete reduit de fa9on spectaculaire : 100 ms 




Energy (MeV) 

Figure IV.4 - Resolution angulaire du LAT en fonction de I'energie dans la classe d'evenements Source, 
pour des photons en incidence normale (ligne pleine), et a 60° d'inclinaison (ligne pointillee), pour des 
conversions ayant lieu dans la region du trajectographe oil les couches de tungstene sont fines. Figure issue 
de |Atwood efoT] ( |2009l l. 
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Energy (MeV) 



Figure IV.5 - Resolution en energie du LAT en fonction de I'energie, pour des photons en incidence 
normale (ligne pleine) et a 60° d'inclinaison (ligne pointillee), dans la classe d'evenements Source. Figure 
extraite de [Atwood et'oL] p009| l . 



pour EGRET, centre 26,5 /is pour le LAT. Les resolutions angulaire et en energie ont egalement ete 
ameliorees. Globalement, la strategic d' observation (balayage plutot que pointe), la surface effective 
et la precision de localisation rendent le LAT 30 fois plus sensible que son predecesseur a des sources 
ponctuelles. 

IV.2.8 Precision de datation 

La precision de datation des evenements mesures par le telescope est importante pour I'etude des 
sursauts 7, necessaire pour I'observation de pulsars jeunes et cruciale pour les pulsars milliseconde. 
C'est pour cette raison qu'avant le lancement de Fermi des tests de la precision de mesure des dates 
ont ete conduits. Une paire de scintillateurs a ete disposee a proximite du satellite de fa^on a ce que 
des muons atmospheriques puissent traverser a la fois le TKR et les scintillateurs. Pour ces derniers, 
les dates etaient enregistrees a I'aide d'un systeme GPS utilise auparavant par le telescope CELESTE 
et dont la precision avait ete demontree grace a des observations du pulsar du Crabe dans le domaine 
optique ( |de Naurois et a/.|2002[ ). 



Les premiers essais de la datation du LAT, monte sur le satellite et utilisant son systeme GPS ont 
ete effectues les 22 et 23 fevrier 200'^ II s'agissait de tester deux configurations : 1) le LAT re9oit et 
interprete le signal GPS pour calculer les dates des muons, 2) on empeche le LAT d'utiliser le signal 
GPS pour calculer les dates, celles-ci sont derivees d'horloges internes. Pour le premier cas, les temps 
enregistres par les deux instruments different, et les ecarts forment un motif de dents de scie au cours 
du temps. Les ecarts s'echelonnent de a — 1 ms, avec une pente de quelques /is/s (un exemple est 



presente en Figure IV.6). Le deuxieme cas est plus critique encore : au lieu de revenir a lorsque 
Ar = — 1 ms, les ecarts continue de croitre, et atteignent quelques dizaines de ms. 

Une erreur de programmation dans les algorithmes de calcul des dates a ete diagnostiquee. Apres 
correction, les tests finaux ont revele que les deux systemes GPS etaient en accord a 300 ns pres (cf. 



Ces tests ont ete realises par D. A. Smith, J. E. Grove, D. Dumora, D. P. Sandora et E. J. Siskind. 
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Figure IV.6 - Tests de la precision de datation du LAT, par comparaison avec un systeme GPS de reference 
(ici, test 77014193 realise le 22 fevrier 2007). A gauche : histogramme des differences de temps entre le 
LAT et les scintillateurs. Le pic legerement en de9a de represente le signal de muons. A droite : ecarts 
entre les dates enregistres par les deux telescopes, en fonction du temps. On observe ici un comportement 
en dents de scie, dont la pente est de -3,4 /is/s. 



Figure IV.7). Le controle des donnees de telemetrie indique que la datation des evenements par le 
LAT sur orbite se comporte comme lorsque le telescope etait au sol (Abdo et al. 2009j[ ). La preci- 
sion de datation est done inferieure a 1 jUs, garantissant la possibilite d'etudier les pulsars, a fortiori 
milliseconde, avec le LAT. 
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Figure IV.7 - Resultats des tests de datation finaux. A gauche : ecarts entre les dates mesurees par les deux 
dispositifs en fonction du temps. A droite : histogramme des differences de temps. L'ecart moyen est de 0,3 
jj.s. Figure extraite de Abdo et al. (2009j i. 
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Observation de pulsars milliseconde 

avec le lat 
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Chapitre V 
Pulsars milliseconde du champ galactique 
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Resume 



DISPOSANT d'ephemerides contemporaines aux observations du LAT pour la plupart des pulsars 
milliseconde situes en dehors d'amas globulaires, nous avons recherche des pulsations en 
rayons 7, et detecte huit MSPs avec une significativite superieure a 5 a. Les huit objets detectes sont 
J0030+0451, J0218+4232, J0437-4715, J0613-0200, J0751+1807, J1614-2230, J1744-1134 et 
J2124— 3358. En plus de ces pulsars, deux MSPs montrent des signaux interessants, J0034— 0534 et 
J1713+0747, et la progression de leur significativite semble indiquer une detection probable dans les 
mois a venir. Nous presentons dans ce chapitre les courbes de lumiere 7 et les proprietes spectrales 
des huit MSPs detectes au cours de cette these. Nous montrons aussi une analyse de donnees XMM 
pour J0030+0451, oia contrairement a des resultats anterieurs, nous determinons la position des pics 
X par rapport aux pics radio. Les composantes radio et X sont en alignement proche, ce qui conforte 
I'idee que remission X de J0030+0451 est produite aux poles, par des processus thermiques. De 
maniere generale, les phasogrammes 7 des huit MSPs ressemblent a ce qui est observe par ailleurs 
pour les pulsars normaux, de meme que leurs spectres. En outre, on note que les MSPs detectes 
dominent cette classe d'objets en terme de E /d^, tout comme chez les pulsars ordinaires. Ces points 
communs laissent a penser que la production de rayonnement 7 dans la magnetosphere des pulsars 
normaux et des MSPs a sensiblement les memes causes et origines. En I'occurence, les modeles 
d'emission dans la magnetosphere exteme sont privilegies chez les MSPs, comme chez les pulsars 
normaux. Ce chapitre constitue une mise a jour des resultats presentes dans |Abdo et aL\ ( 12009f|al ), 
avec ici un lot de donnees plus large. 
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CHAPITRE V. PULSARS MILLISECONDE DU CHAMP GALACTIQUE 



V.l Detection de huit pulsars milliseconde 

Parmi les 72 pulsars milliseconde du champ galactique {cf. Annexe |A]), 66 ont ete suivis dans le 
cadre de la campagne de datation. Nous avons utilise les ephemerides ainsi obtenues et recherche 
des pulsations dans les donnees du LAT, de maniere analogue a la recherche effectuee une decennie 
auparavant dans les donnees d'EGRET ( [Fierro et a/.|1995[ ). 



V.1.1 Analyse des donnees 

Pour la recherche de pulsations des pulsars milliseconde du champ galactique, nous avons utilise 
les donnees enregistrees par le LAT depuis sa mise en service, le 30 juin 2008, jusqu'au 2 juin 2009, 
essentiellement en mode de balayage. Nous avons retenu les evenements d'energie reconstruite supe- 
rieure a 100 MeV, et appartenant a la classe Diffuse. Cette classe regroupe les evenements qui ont la 
plus grande probabilite d'etre de veritables photons 7, comme nous 1' avons evoque dans le chapitre 
precedant. Pour s'affranchir du maximum de photons yd' albedo terrestre (photons emanant des in- 
teractions des rayons cosmiques dans la haute atmosphere terrestre), nous avons exclu les evenements 
dont la direction reconstruite forme un angle de plus de 105° avec le zenith. A cause de remission 
diffuse galactique, causee par I'interaction des rayons cosmiques dans le gaz inter stellaire, les sources 
situees dans le plan galactique (done a faible latitude galactique b) baignent dans un fond de photons 
7 intense. Les sources a haute latitude sont contaminees par un fond extragalactique isotrope, moins 
intense. Pour cette raison, le rayon p des regions d'interet retenues autour de chaque pulsar a ete 
adapte a leur position en coordonnees galactiques. En I'occurence : 

p = 0,5° \b\<lO° 
p = r sinon 

Remarquons neanmoins que nous avons fait une exception pour deux pulsars, PSR J02 18+4232 
et J1614— 2230, dont la latitude galactique est superieure a 10° mais pour lesquels nous avons reduit 
p a 0,5°. J02 18+4232 est en effet situe a 1° du blazar 3C 66A, tres brillant en rayons 7. Quant a 
J1614— 2230, celui-ci est situe dans une region riche en sources 7. La consequence de ce changement 
de coupure angulaire vis-a-vis de la significativite des signaux est discutee dans la suite. 

Nous avons ensuite transfere les dates des evenements selectionnes au barycentre du systeme so- 



laire, selon 1' equation in.4 puis calcule les phases rotationnelles a partir des ephemerides foumies 
par les radiotelescopes. Certains pulsars n'ont pas ete converts, a savoir Jl 125— 6014, J1216— 6410, 
J1618-39, J1629-6902, J1757-5322 et J1933-6211. Pour ces objets nous avons utilise les para- 
metres de rotation foumis en ligne par le catalogue de pulsars de I'ATNF. Ces ephemerides etant 
anciennes et/ou approximatives, nous prendrons soin de considerer les resultats les concemant avec 
precaution. Pour la realisation des etapes decrites ici, nous avons utilise le logiciel TEMP02. Les 
outils standard d' analyse temporelle gtbary et gtpphase presentent en effet un certain nombre de li- 
mitations : ils ne permettent pas de corriger d'un eventuel mouvement propre ou effet de parallaxe de 
chronometrie, d' employer plus de deux derivees temporelles dans le developpement de Taylor de la 



frequence de rotation du pulsar {cf. equation IIL1[ ), ou encore de s'affranchir d' eventuel timing noise 



par un developpement en harmoniques. En pratique, un module de TEMP02 a ete developpe au cours 
de cette these, permettant I'analyse des donnees du LAT a I'aide des fonctionnalites de TEMP02. Le 
module fermijilug.C est livre dans la distribution de TEMPO2Q 

'http ://tempo2. sourceforge.net/ 
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V.1.2 Recherche de pulsations 

Un certain nombre de methodes statistiques sont disponibles pour determiner si un phasogramme 
donne est compatible avec un signal pulse, 1' absence de pulsation correspondant a une distribution 
uniforme en phase. Parmi ces tests d e periodicite, citons p ar exemple le test de de Pearson, le test 
de Rayleigh, ou encore le H-test ( de Jager et a/.||l989 |), ces methodes etant adaptees a differentes 
situations. Nous avons opte pour le H-test, qui presente I'avantage de ne dependre d'aucun parametre 
arbitraire introduisant des nombres d'essais, comme par exemple le nombre de subdivisions du pha- 
sogramme {binning) necessaire au test de x^, ou le nombre d'harmoniques du test Z,^. 

Le H-test est construit de la maniere suivante : soit / la fonction normalisee decrivant I'echantillon 
des A'^ phases mesurees. On peut representer / par une somme de fonctions d : 



(V.2) 



On peut egalement ecrire / comme un developpement en serie de Fourier d'ordre m : 



Les coefficients ocq, Ofyt et /3yt sont donnes par : 



(V.3) 



1 r^^ 1 
2.L f^^^'^ = Tn 



2n 



2it 



f{^)cos{k^)d(j) 
f{^)sm{k^)d^ 



(V.4) 

(V.5) 

(V.6) 
(V.7) 



Le developpement en fonctions 5 (cf. equation V.2) conduit a : 



CXk 



Pk 



1 ^ 
^cos(fc(^/ 



(=1 
1 ^ 

sin(^<^/) 



!=1 



(V.8) 

(V.9) 
(V.IO) 



Soit Z„j I'ecart entre la distribution de phase observee, representee par la fonction f,„, et le phaso- 
gramme uniforme (qui prend pour valeur sur 1' ensemble de la rotation) : 



lit 



In 



d(l> = 2N'£{a^ + Pi) 

k=l 



(V.ll) 
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Comme le montre 1' expression precedente, est un critere de periodicite, qui depend cependant 
du parametre m, a choisir en fonction de la forme de profil attendue. Le H-test remedie a ce probleme 
en recherchant I'harmonique m optimale : 



//= max (z2-4m + 4) (V.12) 

l<m<20 

La probabilite d'obtenir une valeur superieure a // de maniere fortuite a ete evaluee par simula- 
tions Monte Carlo. EUe est donnee par : 

/ 0, 9999755 e-0'39802 h H<23 
P{>H)^/ 1,210597 e-0'4590i//+0,00229 /f2 23 < // < 50 (V.13) 
\ < 4 X 10"^ sinon 

Le H-test s'avere etre parmi les plus efficace^pour des profils d'emission a un pic raisonnable- 
ment large, ou des profils a deux pics de type Vela, Crabe ou Geminga. Neanmoins, notons que le test 
de Rayleigh est plus efficace pour des profils sinusoidaux (type de profil d'emission jusqu'a present 
non rencontre en rayons 7) et que les tests et de sont plus performants pour des profils a plus 
de deux composantes, ou lorsque les pics sont tres etroits (typiquement, lorsque le cycle utile de la 
courbe de lumiere est inferieur a 5%). Ce choix de test de periodicite introduit par consequent un biais 
de selection : le H-test favorise des detections de profils tels que ceux observes auparavant en rayons 
7, par exemple par EGRET. Cependant, en 1' absence de detections fermes de MSPs avant cette etude, 
il etait raisonnable d'utiliser le type de profil des emetteurs 7connus comme hypothese a priori. 

La recherche de pulsations dans le lot de donnees decrit precedemment conduit a la distribu- 



tion de significativites montree en Figure V.l Rappelons que pour deux des pulsars, J02 18-1-4232 et 
Jl 6 14— 2230, un essai supplementaire a ete introduit du fait du changement de coupure angulaire. 
Pour ces deux pulsars, la probabilite de detection fortuite est inferieure a 4 x 10^^ d'apres le H-test 
pour des rayons d'integration p = 0,5°. Un essai supplementaire amene done a P < S x 10^^. Par 
consequent, la significativite est toujours superieure a 5 a, et le changement de jeu de coupure reste 
sans consequence. Huit pulsars milliseconde sont detectes avec certitude au moment de la redaction de 
ce manuscrit (soit un lot de donnees d'environ un an) : PSR J0030-I-0451, J0218-I-4232, J0437-4715, 
J0613-0200, J0751-I-1807, J1614-2230, J1744-1134 et J2124-3358. En outre, deux pulsars sont 
detectes marginalement, avec des niveaux de significativite de 3,7 et 4,2 o : PSR J 17 13-1-0747 et 
J0034-0534. 

Avant de presenter les courbes de lumiere et les proprietes spectrales des huit MSPs detectes avec 
certitude, ainsi que de discuter des candidats a la detection au cours des mois a venir, remarquons que 
la majorite des 72 MSPs traites ici presentent un signal dont la significativite est inferieure a 3 C7. Pour 
les huit pulsars detectes, on s'aper9oit que la valeur de H augmente de fa^on lineaire avec le temps. 



en premiere approximation (la Figure V.2 illustre quelques exemples). En supposant qu'un signal 
actuellement compris entre et 1 C7 (soit // < 2,9) voie sa valeur de H augmenter lineairement, on 
aboutit aH < 14, 5 au bout de cinq ans. Or une detection a 3 c correspond a une valeur de H d'environ 
15. On peut par consequent avancer que les pulsars dont la significativite n'a pas atteint 1 a au cours 



^L'efficacite d'un test de periodicite est definie comme etant la probabilite qu'une source periodique soit identifiee par 
ce test en tant que telle au-dessus d'un certain seuil de vraisemblance. Cette efficacite depend des proprietes de la source 
etudiee : nombre et largeur des pics, intensite, etc. Cf. \de Jager et a/.| ( |l989| l pour une comparaison de I'efficacite du H-test 
avec celle d'autres tests de periodicite. 
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Figure V.l - Resultats de la recherche de pulsations pour les 72 MSPs du champ galactique. Les noms des 
pulsars pour lesquels la significativite est d'au moins 3, 4, ou 5 (7 sont indiques. 



de cette premiere annee ne seront pas detectes dans les cinq premieres annees de fonctionnement du 
LAT, avec les jeux de coupures utilises ici. En revanche, pour les pulsars dont le signal est compris 
entre 1 et 2 cr (soit H < 7,76) le meme raisonnement conduit a des valeurs de H inferieures a 38,8 
apres cinq ans. Cette valeur correspond approximativement a une detection a 5 C7, par consequent le 
LAT pourrait detecter certains de ces objets dans les quatre annees a venir. 

Cependant, les 62 MSPs pour lesquels la significativite est inferieure a 3 C7 ne presentent pas 
de telles croissances lineaires de H en fonction du temps, indicatrices de detections a venir. Pour 
les futures recherches de pulsations y dans les donnees LAT, il faudra done envisager des jeux de 
coupures differents, en angle ou en energie, ou opter pour un test de periodicite different : rappelons 
en effet que le H-test est efficace pour des profils d'emission a un ou deux pics relativement larges, 
mais son efficacite chute pour les pics etroits ou nombreux. 



V.1.3 Observations X et 7 de PSR J0030+0451 

Le premier MSP pour lequel des pulsations en rayons 7 ont ete detectees au cours de cette these 
est PSR J0030+045 1 . Ce pulsar, de periode P = A, 865 ms, a ete decouvert independamment par deux 
recherches de pulsars dans le domaine radio, conduites respectivement a Arecibo ( |Somer|[2000 1 et a 



Bologne ( [D'Amico 2000). La chronometrie a conduit a une m esure de parallaxe de 3,3 ±0,9 mas, 
correspondant a une distance de 300^gg*^ pc ( Lommen et a/. [2006^ . Le mouvement apparent du pulsar 



est relativement faible. Les valeurs de P apparentes et intrinseques different de moins d'un pourcent. 
Les quantites derivees E, T, Bs et Bjjq sont donnees en Annexe [a} 

Ce pulsar n'etait pas encore connu lors de la recherche de MSPs dans les donnees d'EGRET 
( [Fierro et al.\\995) . Cependant il est interessant de noter que EGR J0028+0457, appartenant au cata- 
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Figure V.2 - Exemples de progression de la valeur de H en fonction du temps, pour quatre des huit MSPs 
detectes. La progression est lineaire, en premiere approximation. 



logue de sources d' EGRET revise, coincide avec la position du pulsar ( [Casandjian & Grenier|2008| ). 
Les perspectives de detection de J0030+0451 avec le LAT etaient done importantes. 

Par ailleurs, PSR J0030+0451 a egalement ete detecte en rayons X en tant que source pulsee, 
d'abord par ROSAT ( [Becker et aL\pmU^ , puis par XMM-Newton ( [Becker & Aschenbach] [20021 ). 
Les deux telescopes X ont observe un profil d'emission comprenant deux pics larges, separes d' en- 
viron 180°, mais n'ont pas permis de determiner le decalage en phase des composantes radio et X, 
les horloges de ces telescopes ne le permettant pas a I'epoque des observations. Depuis, des correc- 
tions ont ete apportees a la calibration des horloges embarquees par XMM-Newton, de sorte que 
celles-ci delivrent desormais des dates dont la precision absolue est 300 /is, soit 0,06 rotations de 
J0030-I-045 1[^ Parallelement a I'analyse des photons 7 du LAT, nous avons re-analyse les observa- 
tions de J0030+0451 par XMM, afin de mesurer I'alignement relatif entre les emissions radio, X et 

r- 

PSR J0030-I-0451 a ete observe par XMM les 19 et 20 juin 2001, pour une duree totale de 29 ks. 
Pour construire une ephemeride radio couvrant a la fois les observations recentes de Fermi ainsi que 
celles de XMM, nous avons utilise 700 TOAs radio enregistres au radiotelescope de Nan9ay entre 
juillet 1999 et aujourd'hui. L'essentiel des observations a ete realise a 1,4 GHz (1360 ± 72 MHz 
avant 2002, 1398 ± 32 MHz apres 2004 et 1398 ± 64 MHz depuis juillet 2008), en complement 
d' observations a 2 GHz afin de contraindre au mieux la mesure de dispersion. Cette derniere induit 
en effet un retard entre la composante radio et remission de haute energie, qu'il faut corriger pour 



pouvoir mesurer le decalage intrinseque (cf. equation III.5). L'incertitude moyenne sur les TOAs est 
de 3,6 jUs pour les donnees enregistrees apres 2004, et 8,6 /is pour I'ensemble du lot de donnees. Les 



residus de chronometrie sont montres en Figure V.3 La rms des residus est de 3,7 [is. Nous mesurons 



ainsi une parallaxe chronometrique de 4,1 ± 0,7 mas, comparable a celle de Lommen et al. (2006). 



cf. http ://xmm2.esac.esa.int/docs/documents/CAL-TN-0045-l-0.pdf 
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Figure V.3 - Residus en fonction du temps pour le pulsar milliseconde J0030+0451, observe a Nan9ay. Ces 
TOAs radio couvrent a la fois les observations du MSP par XMM, en juin 2001, et les recentes observations 
du Fermi-LAT. 



De meme, la valeur de DM obtenue dans cette etude, de 4,333 ± 0,001 cm^-^ pc, est en accorcQ 

Pour r analyse des donnees X, nous avons utilise le signal enregistre par I'instrument EPIC-pn 
en mode de datation. Nous avons selectionne les evenements d'energie comprise entre 300 eV et 
2,5 keV. Une precision temporelle optimale est obtenue en enregistrant le signal des CCD du pn en 
continu, ce qui a pour effet d'etaler I'information dans la direction y. Par consequent, les evenements 
sont conserves en ecrasant les donnees dans une seule dimension, oii chaque pixel de la direction 
X contient la somme des pixels de la direction y. Nous avons selectionne une region centree sur le 
pulsar, large de sept pixels dans la direction x. Le niveau de fond a ete estime en utilisant une region 
voisine libre de sources X et de meme dimension. Nous 1' avons ensuite moyenne sur le nombre de 
subdivisions du phasogramme X et soustrait au phasogramme. Enfin, les dates ont ete transferees au 
barycentre du systeme solaire grace a I'outil dedie barycei^ 



Les courbes de lumiere radio, X et 7 de J0030+0451 sont presentees en Figure V.4 Le profil 
d'emission ycomprend deux pics, separes de A = 0,45 ± 0,01. Comme on I'a vu dans le deuxieme 
chapitre, ce type de profil a deux pics separes de 0,4 a 0,5 rotations est relativement commun parmi les 



pulsars emetteurs 7 comme Vela, le Crabe, Geminga ou encore PSR J2021+3651 ( Thompsonl|2004 



Abdo et al. 2009e). Le premier pic 7 est decale de 5 = 0,16 du maximum de remission radio, ce 
qui est egalement une caracteristique commune : d est compris entre 0,11 et 0,16 pour Vela, PSR 
B 195 1+32 ou J2021+3651. La separation des pics et le decalage radio / 7indique que le rayonnement 
7 est vraisemblablement produit dans la magnetosphere exteme, au contraire de remission radio. 



La Figure V.4 montre que remission en rayons X est en alignement proche avec la composante 
radio (rappelons que la precision temporelle est de 300 /is, soit 0,06 fractions de tour), et done decalee 
de remission 7. Ce resultat indique que les emissions X et radio ont des origines communes dans la 
magnetosphere, et que le rayonnement 7 est produit dans une region differente. L' alignement radio/X 
conforte I'idee que PSR J0030+0451 fait partie des emetteurs de rayonnement X thermique, a I'instar 

^La construction de TOAs a partir des observations NBPP et BON de J0030+045 1 a ete realisee par I. Cognard. 
^La selection et la barycentrisation des donnees XMM ont ete realisees par N. Webb et B. Pancrazi. 
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Figure V.4 - Courbes de lumiere radio, X et yde PSR J0030+0451. Deux rotations du pulsar sont montrees. 



de J0437— 4715, pour lesquels remission est produite pres des calottes polaires, la ou la temperature 
est la plus importante ( Zavlin|2007 ). Par ailleurs, 1' analyse spectrale des donnees XMM a montre que 



remission est compatible avec un modele purement thermique ( Becker & Aschenbach|2002[ ) 



V.1.4 Confirmation de PSR J0218+4232, detection de six autres MSPs 

En plus de la detection de PSR J0030+0451, sept autres MSPs ont ete detectes avec le LAT au 



cours de cette these. Leurs courbes de lumiere 7 sont presentees en Figures V.5 et V.6 Remarquons 



neanmoins que ces courbes de lumiere sont basees sur des coupures angulaires fixes, ce qui a ten- 
dance a biaiser les profils 7 vers les hautes energies, compte tenu de la resolution angulaire du LAT. 
Des courbes de lumiere obtenues avec des coupures angulaires dependantes de I'energie sont presen- 
tees en Annexe |Cj Les phasogrammes globaux au-dessus de 100 MeV et de 1 GeV sont semblables 
qualitativement avec ceux presentes en Figures [V5]et|V.6[ avec le lot de donnees actuel. 
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Les sept pulsars sont les suivants : 

- J0218+4232 : ce MSP appartenant a un systeme binaire a ete decouvert a Westerbork puis 
confirme au radiotelescope de Jodrell Bank ( [Navarro et a/.|[T995] ). Comme nous I'avons evo- 
que dans le second chapitre, des pulsations marginales ont ete decouvertes dans les donnees 
d'EGRET ( |Kuiper et ar][200^ . La Figure |V5] montre que le profil 7 de J0218+4232 est tres 
large, une caracteristique que Ton retrouve dans le domaine radio. Le pic principal radio, ici 
centre autour de 0, possede deux composantes symetriques et la structure globale est large (~ 
0,4). On constate que les profils radio et 7 sont similaires : remission 7 forme un faisceau large, 
et a haute energie (E > 1 GeV) deux pics semblent emerger, bien que le nombre de photons ne 
permet pas pour le moment de I'affirmer. Les profils radio et 7 sont decales d'environ 0,5. En 
raison du faible nombre de photons et de la largeur du profil 7, il est difficile pour le moment 
de determiner le decalage radio / 7 avec plus de precision. Nous avons done confirme la detec- 
tion de J02 18+4232 par EGRET, bien que le phasogramme 7 obtenu dans cette etude differe 
qualitativement de celui d'EGRET (cf. Figure IL6 1. 

- J0437— 4715 : decouvert avec le radiotelescope de Parkes (Johnston et al. 1993), ce MSP est 
le plus proche connu a present. Malgre sa faible distance, la recherche de pulsations 7 dans les 
donnees d'EGRET n'a pas abouti ( [Fierro et a/.|1995[ ). La contribution de I'effet Shklovskii pour 
ce pulsar est importante, si bien que la valeur de E intrinseque est environ quatre fois inferieure 
a la valeur apparente. Le pulsar presente un seul pic 7, separe de 5 = 0,45 avec le pic radio. 
Comme chez PSR J0030+045 1, remission X est alignee avec la composante radio, et le spectre 
est compatible avec un rayonnement X thermique ( [Zavhn 2007), tandis que remission 7 est 
decalee, implicant des regions de production separees. 

- J0613— 0200 : ce pulsar situe dans un systeme binaire a ete decouvert au radiotelescope de 
Parkes ( Lorimer et al. [1995 ). Le profil d'emission 7est analogue a celui de J0437— 4715. Le pic 
7 est ici separe du maximum de remission radio de 5 = 0,42. En supposant que la similitude 
avec J0437— 4715 persiste aux autres longueurs d'onde, on pourrait s'attendre a ce que des 
observations en X revelent un pic en phase avec remission radio et un spectre thermique. 

- J0751+1807 : parmi les objets presentes ici, ce pulsar est a part : il a en effet ete decouvert en 
radio dans la boite d'erreur d'une source 7non identifiee du premier catalogue d'EGRET, GRO 
J0749+17 (Lundgren et al. 1995). La luminosite de la source d'EGRET etait cependant trop 
importante pour etre due a la seule contribution de PSR J0751+1807. Par la suite, la recherche 
de pulsations dans les donnees d'EGRET n'a pas permis la detection du pulsar (Fierro et al 
1995). L' analyse de donnees XMM a revele des pulsations marginales ( Webb et q/.||2004 ). Le 
decalage des composantes radio et X n'a cependant pas pu etre mesure, les horloges ne le 
permettant pas. Le phasogramme 7de J0751+1807 vu par le LAT possede des caracteristiques 
communes avec celui de J0218+4232 {cf. Figure V.6). En effet pour ce pulsar le faisceau radio 
presente deux sous-pics formant une largeur totale d'environ 0,2, largeur que I'on retrouve dans 
le pic d'emission 7 ici centre autour de 0,5. La structure du faisceau 7 semble presenter deux 
pics, en particulier a haute energie, la separation n'etant pas significative avec le lot de donnees 
actuel. Comme pour J0218-I-4232, le decalage 5 est d'environ 0,5. Avec plus de photons, il sera 
possible de mesurer le decalage plus precisement et de conclure sur la presence d'un eventuel 
deuxieme pic. 

- J1614— 2230 : ce MSP a egalement ete decouvert au cours d'une recherche de pulsars dans les 
boites d'erreur de sources EGRET non identifiees ( [Crawford et a/. [2006] ). Cette fois il s'agissait 
de la source du troisieme catalogue 3EG J1616— 2221. La courbe de lumiere 7 presente deux 
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Figure V.5 - Courbes de lumiere 7 au-dessus de 1 GeV et de 100 MeV, et profil radio pour quatre des huit 
MSPs detectes par le LAT. Deux rotations sont montrees. 
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pics, separes de A = 0,5 et dont le premier succede au pic radio de 5 = 0,2. Comme pour 
J0030+0451, le decalage radio / 7 et la separation des pics sont semblables a ceux des pulsars 
jeunes. En revanche, si un bridge est clairement observe chez J0030+0451, la statistique est 
trop faible a present pour determiner si J1614— 2230 presente une emission inter-pics. 
- J1744— 1134 : ce pulsar milliseconde isole a ete decouvert par le radiotelescope de Parkes 
( [Bailes et al]|1997| ). ROSAT a montre que la source RX J1744.4-1134 est vraisemblable- 



ment associee au pulsar milliseconde, mais le nombre de photons enregistres etait insuffisant 
pour la recherche de pulsations (Becker & Triimper 1999| ). La complexite de remission 7 de 



J1744— 1134 a haute energie est frappante. Un pic principal se degage, separe de 5 = 0,85 
du pic radio. II est pour I'instant difficile de denombrer les pics. Cependant le phasogramme 
au-dessus de 1 GeV suggere qu'un certain nombre de composantes secondaires existent. 
- J2124— 3358 : tout comme J1744— 1134, ce pulsar est isole et a ete detecte en radio a Parkes 
( Bailes et q/.||1997| ). Des pulsations ont ete decouvertes en rayons X grace aux telescopes RO- 



SAT (I Becker & Trumper| [T999) et XMM (Zavlin 2006]). Ces instruments ont mis en evidence 



un profil d'emission X comportant un pic principal asymetrique, mais n'ont malheureusement 
pas pu mesurer le decalage radio/X. La complexite de remission radio de PSR J2124— 3358 
est remarquable, suggerant que le pulsar emet tout au long de sa rotation. Un pic principal est 
neanmoins observable en radio, que Ton trouve en rayons 7 decale de 5 = 0,85 comme chez 
J1744— 1134. Un seul pic 7 se degage dans les donnees actuelles. Etant donnee la structure 
complexe de remission radio, la courbe de lumiere 7pourrait egalement exhiber de multiples 
composantes. 



La Table V.l liste les valeurs de 5 et A. En resume, J0030-I-0451 et J1614— 2230 possedent des 
profils 7 a deux pics separes de A ~ 0,45 et retardant de 5 ~ 0,15 par rapport a remission radio, des 
caracteristiques identiques a celles de beaucoup des pulsars 7ordinaires. J0437— 4715 et J0613— 0200 
ne presentent qu'un pic, decale d'environ 0,4 par rapport a la radio, ce qui les rend qualitativement 
similaires a J2229-I-61 14, pulsar 7jeune ( Abdo et a/.|2009i| ). Pour les quatre MSPs restants, le nombre 



de photons est encore insuffisant pour conclure sur la multiplicite des pics et leur position. 

V.l. 5 Detections a venir ? 



Outre les huit detections de MSPs mentionnees precedemment, la Figure V.l indique que PSR 
J0034— 0534 et J1713-I-0747 sont detectes marginalement, avec des significativites de 4,2 o et 3,7 a, 
respectivement. On a vu que pour les detections certaines, la valeur de H croit de fa9on lineaire (cf. 



Figure V.2 1. Le meme graphe de progression de H en fonction du temps constitue par consequent un 
indicateur des detections a venir : on s'attend en effet a ce qu'un pulsar detectable voie sa significati- 
vite augmenter regulierement avec le temps, et que la significativite varie de fagon erratique pour un 
pulsar non emetteur. 



La Figure V.7 montre, pour les deux pulsars marginalement detectes, la progression de H en 
fonction du temps. Pour les deux pulsars, H est globalement croissant (en particulier pour PSR 
J0034— 0534, oil la croissance est quasi-monotone). L'allure des courbes amene a penser que ces 
pulsars sont bel et bien emetteurs 7. Des ajustements lineaires sont montres par des lignes pointillees 
obliques. lis atteignent H = 39,5 (equivalent k 5 o) entre 2,85 et 2,95 x 10^ MET s. Les pulsars de- 
vraient done etre detectes par le LAT au cours des mois a venir en utilisant le meme jeu de coupure, 
en supposant que ces regimes de progression se maintiennent. Nous verrons par la suite que ces deux 
pulsars sont associes a des sources d'emission continue detectees par le LAT, renfor§ant I'hypothese 
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Figure V.6 - Courbes de lumiere 7 au-dessus de 1 GeV et de 100 MeV, et profil radio pour les quatre autres 
MSPs detectes par le LAT. Deux rotations sont montrees. 



92 



y. 1 . DETECTION DE HUIT PULSARS MILLISECONDE 



selon laquelle les deux MSPs sont emetteurs de rayons 7. 
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Figure V.7 - Progression de la valeur de H en fonction du temps, pour les pulsars milliseconde J0034— 0534 
et J171 3+0747. La valeur de H equivalant a une detection a 5 C7 est indiquee par une Ugne pointillee 
horizontale, et un ajustement lineaire est montre par une ligne pointillee oblique. 



Les courbes de lumiere actuelles pour PSR J0034-0534 et J1713+0747 au-dessus de 100 MeV 



et de 1 GeV sont presentees en Figure V.8 ainsi que leur profil d'emission radio a 1,4 GHz. Pour 
J0034— 0534, il semble que le profil d'emission ycomprenne deux pics, separes d'environ 0,3, et en 
alignement avec les pics radio. Dans la mesure oii les huit MSPs detectes a present ont leurs pics 7 
et radio en decalage, J0034— 0534 pourrait etre le premier MSP a avoir ses pics alignes dans les deux 
gammes d'energie, tout comme le Crabe. II pourrait egalement etre interessant d' observer le pulsar 
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en rayons X. Jusqu'a present les recherches en X se sont averees infructueuses ( |Zavlin|2006 ). 

En ce qui concerne PSR J17 13+0747, la description est moins aisee, bien qu'une structure semble 
se degager vers 0,4 en phase. Si le MSP est confirme en tant qu'emetteur 7 pulse, il semble a present 
qu'il entrerait dans la categoric des pulsars a un seul pic radio et 7, separes d'environ 0,4 en phase, 
comme J0437-4715, J0613-0200 ou encore le pulsar jeune J2229+61 14 ( ,Abdog? a/.|2009i| ). 
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Figure V.8 - Phasogrammes yde PSR J0034-0534 et J1713+0747, probables futures detections de MSPs 
par le LAT. Les cadrans superieurs montrent les courbes de lumiere au-dessus de 100 MeV et 1 GeV. Les 
parties inferieures montrent les profits radio a 1,4 GHz, acquis avec le radiotelescope de Nanfay. 



V.1.6 Analyse spectrale 



L' analyse spectrale de s huit MSPs detectes dans cette etude a ete realisee a I'aide des methodes 



decrites dans 



Abdo et al. 



{2009 gj^ De fa9on schematique, remission 7 est ajustee par un modele 
comprenant des fonds galactique, extragalactique et instrumental (ces deux demieres sources de fond 
etant isotropes), ainsi que les sources avoisinantes et le pulsar lui-meme, dont on souhaite mesurer les 
proprietes spectrales. La fonctionnelle utilisee pour I'ajustement du spectre des pulsars milliseconde 
est une loi de puissance a coupure exponentielle, conformement a I'equation IL2[ avec le parametre 
de forme b egal a 1 : 



dN 
dE 



No 



1 GeV 



e 



(V.14) 



^Dans cette partie, les resultats spectraux presentes sont ceux presentes dans Abdo et al. ( 2009a 1, correspondant a un 
lot de donnees s'etendant jusqu'au 15 mars 2009. 
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Pulsar 


P 


8 


A 


/ (> 100 MeV) 


h (> 100 MeV) 


r 


Ec 


Ly 




(ms) 






(10-'* 
ph/cm^/s) 


(10-" 
erg/cm^/s) 




(GeV) 


(10^2 gj.g/s) 


J0030+0451 


4,865 


0,16 


0,45 


5,5 ± 0,7 


4,9 ± 0,3 


1,3 ± 0,2 


1,9 ±0,4 




J02 18+4232 


2,323 


0,5 


- 


5,6 ± 1,3 


3,5 ± 0,5 


2,0 ± 0,2 


7 ±4 


(3,0 ± 0,9) X 10^ 




^,151 


0,45 




4,4 ± i,U 


1,9 ± 0,3 


2,1 ± 0,3 


2,1 ± 1,1 


n _i_ r\ r\c\ 
0,50 ± 0,09 


J0613-0200 


3,062 


0,42 




3,1 ± 0,7 


3,1 ±0,3 


1,4 ± 0,2 


2,9 ± 0,7 


8 4+''-** 


J0751+1807 


3,479 


0,5 




2,0 ± 0,7 


1,7 ±0,2 


1,6 ± 0,2 


3,4 ± 1,2 


T 7+10,9 
'''-3,8 


JI614-2230 


3,151 


0,2 


0,5 


2,3 ± 2,1 


2,5 ± 0,8 


1,0 ± 0,3 


1,2 ±0,5 


50 ±22 


J1744-1134 


4,075 


0,85 




7,1 ± 1,4 


4,0 ± 1,0 


1,5 ± 0,2 


1,1 ±0,2 


nil 


J2124-3358 


4,931 


0,85 




2,9 ± 0,5 


3,4 ± 0,3 


1,3 ± 0,2 


2,9 ± 0,9 


2,5!?:^ 



Tableau V.l - Proprietes des MSPs detectes avec le Fermi-LAT. Le parametre 5 correspond au decalage en 
phase entre le maximum de remission radio et le pic y le plus proche, A est la separation des pics y pour 
les profils a deux pics. Les parametres / et h donnent les flux en photons et en energie, au-dessus de 100 
MeV. Enfin, F, et Ly — sont I'indice spectral, I'energie de coupure et la luminosite y au-dessus 

de 100 MeV. Les incertitudes citees ici sont statistiques. Les erreurs systematiques, liees a I'incertitude sur 
la reponse de I'instrument et sur le fond diffus sont (-0,1 ; +0,3) pour F, (-10% ; +20%) pour Ec, (-10% ; 
+30%) pour / et (-10% ; +20%) pour h. 



Rappelons que dans cette expression E represente I'energie, T est I'indice spectral et Ec est I'ener- 
gie de coupure du spectre. Le flux dN/dE correspond au nombre de photons provenant de la source 
et detecte par unite de temps, d'energie, et de surface du detecteur. La reponse de I'instrument (reso- 
lution angulaire, surface effective, resolution en energie) est variable en fonction de I'energie, comme 
on I'a evoque dans le chapitre precedent ; elle doit done etre prise en compte par I'estimateur lors du 
calcul de Nq, T et Ec. La fonction de reponse de I'instrument utilisee pour cette analyse spectrale est 
P6_V1_DIFFUSE. Ces estimations de la reponse du LAT ont ete etablies avant lancement. II semble 
a present que ces IRFs sous-estiment les flux mesures, comme indique par les erreurs systematiques 



citees en legende de la Table V.l Par integration de dN/dE sur I'energie, on peut calculer le flux de 
photons / au-dessus de 100 MeV, ainsi que le flux en energie h, par : 



J 100 MeV \dt J 



100 MeV V"^ 



(V.15) 
(V.16) 



La Table V.l donne les valeurs de F, Ec, f &t h pour les 8 MSPs. Deux exemples de spectres 



en energie sont presentes en Figure V.9 PSR J0030-I-0451 est le pulsar milliseconde le plus brillant 
vu de la Terre; son analyse spectrale est relativement aisee. PSR J1614— 2230 est un exemple de 
pulsar moins lumineux pour lequel les parametres spectraux sont plus difficilement contraints, la 



consequence etant ici une erreur statistique importante sur le flux de photons / {cf. Table V.l I, en 



particulier a cause du mauvais ajustement a basse energie, oil la m ajorite des photons sont emis. 
Les valeurs de luminosite y Ly = Aitf^hd^ {cf. equation IL3 1 sont egalement indiquees dans la 



Table V.l En I'absence d'hypothese sur la geometric de remission y des MSPs, nous avons utilise 
= 1 . Notons neanmoins que fa est certainement variable d'un pulsar a 1' autre, aussi une incertitude 
importante est introduite dans le calcul de Ly. Pour six des huit MSPs, la distance est fiable, provenant 



de mesures de parallaxes trigonometriques ou chronometriques (Lommen et al. 2006, Hotan et al. 
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Figure V.9 - Distributions spectrales en energie pour les pulsars J0030+0451 et J1614— 2230, avec leurs 
ajustements par des lois de puissance a coupure exponentielle. Les parametres de ces ajustements sont 
donnes en Table |V.l] Credit : T. J. Johnson, M. Kerr. 
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Erreur (dB) 



Figure V.IO - Histogramme des erreurs sur les distances derivees du modele NE2001, pour les pulsars 
ayant une mesure de parallaxe. Les erreurs sont nuUes en moyenne, mais ont une deviation standard de 2 



dB et s'etendent jusqu'a 6 dB. Figure extraite de Deller (2009 1. 



20061 IDeller et al. |2008D . Les distances de PSRJ02 18+4232 etJ16 14-2230 sontbasees sur le modele 
NE2001 de densite d'electrons dans la Galaxie ( |Cordes & Lazio|[2002 l). Si les distances derivees du 
modele NE2001 sont fiables en moyenne, des comparaisons particulieres avec des distances issues 
de mesures de parallaxes ont montre que NE2001 produit des resultats errones dans de nombreux 
cas (cf. Figure V.IO). L'efficacite d'emission y au-dessus de 100 MeV, rj = Lj/E, varie de 2% pour 



J0437-4715 a environ 30% pour J1744-1134, a I'exception de PSR J1614-2230 pour lequel 7] 
est d'environ 100%. Ceci indique que la distance du pulsar est vraisemblablement surestimee par 
le modele NE2001. Une distance plus faible reduirait I'effet Shklovskii dont la contribution a P est 
Pshk = l^^d/ (cP), augmentant ainsi E ; tout en diminuant Ly. La valeur de l'efficacite s'en trouverait 
ainsi reduite. Si Ton ne tient pas compte de la variation de E due a I'effet Shklovskii, une efficacite 
de 10% requiert une distance inferieure d'un facteur ~ 3, amenant J1614— 2230 a J ~ 0,4 kpc. 

Bien que le cas de figure ne soit pas prevu par les differents modeles theoriques, remarquons 
qu'il sera judicieux a I'avenir d'ajuster les spectres des MSPs en introduisant le parametre de forme 



b (cf. equation II. 2 1, et ainsi s' assurer que celui-ci est compatible avec 1. Les donnees actuelles ne 
permettent pas a present d'effectuer cette verification et d'en tirer des conclusions significatives. 



V.1.7 Emission continue 

La recherche de sources d'emission continue parmi les 72 MSPs galactiques est complementaire 
de la recherche de pulsations. En effet elle fournit un a priori fort sur les pulsars pour lesquels la 
detection de pulsations est envisageable ou ne Test pas. 
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Tableau V.2 - Table des quinze MSPs galactiques associes a des sources continues, dont la significativite 
est superieure a 5 CT. Les objets pour lesquels un signal pulse est detecte sont indiques par les lettres p 
(detection certaine) ou m (detection marginale). 



On peut resumer la methode mise en place dans le cadre de la construction des catalogues de 
sources continues de Fermi de la fa9on suivante : dans un premier temps on recherche des regions 
du ciel oil la densite de photons detectes devie d'un modele de fond galactique et extragalactique 
(un seuil est fixe sur la deviation pour limiter le nombre de sources - par exemple 4 c). La seconde 
etape consiste a affiner la position de chaque source candidate et determiner 1' incertitude sur celle- 
ci. L' etape finale est I'estimation de la significativite de remission continue observee, c'est-a-dire 
la probabilite que cette emission provienne de fluctuations statistiques. Pour cela un algorithme de 
vraisemblance compare les probabilites des situations {fond uniquement} et {fond + source hypo- 
thetique} ou la source est representee par une loi de puissance (soit deux parametres : flux et indice 
spectral). La quantite TS = 2 Aln(vrai semblance) est une mesure de la significativite de la source. La 
probabilite d'obtenir au moins TS est environ egale a 0,5 fois I'esperance associee a une distribution 
de ^ deux degres de liberte, de sorte qu'une significativite de 5 a correspond a ~ 30. Pour plus 
de details, le lecteur pourra se reporter a Abdo et al. (2009dl. La Table V.2| donne les valeurs de TS 



pour les MSPs associes a des sources d'emission continue dont la significativite est superieure a 5 a. 

On note que les huit detections pulsees de MSPs sont egalement observes comme sources conti- 
nues. De meme, J0034— 0534 et J1713+0747 sont associes a des sources continues, ce qui renforce la 
perspective de detection pulsee pour ces deux pulsars, comme presente precedemment. 

La situation est moins claire pour ce qui est des cinq pulsars restants, J0610— 2100, J062 1+1002, 
J1600-3053, J1939+2134 et J1959+2048. Aucun de ces cinq MSPs, pour lesquels des ephemerides 
contemporaines etaient disponibles lors de cette etude, n'est detecte comme source 7 pulsee, y com- 
pris marginalement. Davantage de photons et un jeu de coupures different pourraient permettre la de- 
tection de ces pulsars. En outre, la nebuleuse a vent de pulsar (Pulsar Wind Nebula, PWN) G59.2— 4.7, 
emettrice de rayons X, est associee a J1959+2048 ( [Stappers et a/.|2003] ). Cette nebuleuse pourrait etre 



responsable de remission continue en 7. Une analyse detaillee des proprietes spectrales et spatiale de 
remission pourrait reveler la nature de la source. 



98 



W.2. DISCUSSION 



V.2 Discussion 



V.2.1 Proprietes des MSPs detectes 



Les courbes de lumiere 7 des pulsars milliseconde et ralignement relatif entre les composantes 
radio et 7 sont similaires a ce qui est observe par ailleurs pour les pulsars normaux. Par exemple, 
les profils a deux pics de J0030+0451 et de J1614— 2230, separes de A ~ 0,45 et dont le premier pic 
retarde de d ~ 0,15 par rapport a la radio sont relativement communs parmi les pulsars ordinaires 
( |Abdo et Ql]|2009g|h|C,e) . Un plus grand nombre de MSPs detectes presente un seul pic, separe de 
5 ~ 0,4 - 0,5 de remission radio, a I'image de J2229+6114 ( |Abdog? a/.|2009il ). On en conclut que la 
geometric de remission 7 dans la magnetosphere des MSPs d'une part et des pulsars jeunes d' autre 
part est relativement similaire. 



L' analyse spectrale, dont les resultats sont donnes dans la Table V.l a montre que les indices sont 



generalement durs, les valeurs de F etant typiquement inferieures a 2. Les indices spectraux sont done 
comparables a ceux des pulsars normaux. Le spectre le plus dur est celui de PSR J1614— 2230, dont 
I'indice a ete mesure a 1,0 ± 0,3. Le durcissement des spectres en fonction de I'age caracteristique 



avait ete observe pour les pulsars d'EGRET (Fierro et al. 1993 1. II est interessant de remarquer que 
les valeurs de F mesurees pour les MSPs semblent decorrelees de Bs, B^c et E. La Figure V.l 1 (A) 
montre I'age caracteristique x = P / {IP) en fonction de I'indice spectral F pour les MSPs detectes. 
Aucune correlation ne semble se degager entre ces deux grandeurs. Cependant, pour ces pulsars dont 
la rotation a ete acceleree par un processus de recyclage, I'age caracteristique n'est pas representatif 
de I'age veritable de ces objets, comme on I'a evoque precedemment. Dans ces conditions, une corre- 
lation eventuelle entre I'age et I'indice spectral est probablement indecelable sans une connaissance 
de I'histoire d' accretion. 

Les energies de coupure exponentielle Ec sont comprises entre 1 et 4 GeV, si Ton fait exception 
de J0218+4232 pour lequel I'incertitude sur Ec est importante {Ec = 1±4 GeV). A nouveau, ces 
caracteristiques sont comparables a celles des pulsars emetteurs 7 normaux. Ces energies de coupure 
ne presentent pas de correlation particuliere avec B^c et E. Comme pour I'indice spectral des pulsars 
vus par EGRET, il est apparu que I'energie de coupure des pulsars de CGRO semble diminuer lorsque 



^5 augmente {cf. Figure 4 de Thompson (2008b I). II faut cependant mettre un bemol, car la tendance 
repose encore une fois sur les valeurs minimale et maximale, celles de B1509— 59 et B 1951+32. Pour 
ce dernier, 1' analyse des donnees du LAT revele que Ec a ete surestimee par EGRElQ La Figure 
V.l 1| (B) montre les valeurs de Ec en fonction du champ magnetique surfacique, Bs- La dispersion 



est grande, principalement a cause des incertitudes sur Ec- Neanmoins, bien qu'aucune correlation 
claire ne se degage d'apres le graphe, il semble que la tendance vue pour les pulsars de CGRO ne soit 
pas confirmee, a savoir, I'energie de coupure exponentielle Ec ne diminue pas en fonction du champ 
magnetique surfacique Bs- 

La luminosite 7 au-dessus de 100 MeV, Ly = Anfohd^ est representee en fonction de E pour les 



MSPs et les pulsars normaux detectes jusqu'a present dans la Figure V.12 Pour les MSPs, une valeur 
de 1 a ete utilisee pour le facteur de correction du flux /n. Les valeurs de /f^ predites par le modele 
Outer Gap ont ete adoptees ici pour les pulsars normaux (1' utilisation des facteurs de correction issus 
du modele Slot Gap ne modifierait pas qualitativement la Figure V.12[ ), d'apres Watters et al. (2009 1. 
Une ligne pointillee, representant un ajustement de la luminosite 7 des pulsars normaux par une loi en 
est egalement tracee. La luminosite de J02 18+4232, un MSP tres energetique, est en bon accord 



^Information apportee par T. Reposeur. L' analyse de B 195 1+32 par le LAT sera publiee prochainement. 



99 



CHAPITRE V. PULSARS MILLISECONDE DU CHAMP GALACTIQUE 




T (yrs) 

(a) r en fonction de T 



Bs (G) 

(b) Ec en fonction de Bs 



Figure V.ll - Indice spectral en fonction de I'age caracteristique T et energie de coupure Ec en fonction 
du champ magnetique surfacique Bs pour les huit MSPs detectes, oil T et Bs ont ete corriges de I'effet 
Shklovskii. 



avec cette loi. Les sept autres MSPs sont distincts de J0218+4232 de par leur faible E. Pour ces 
objets (a I'exception de J1614— 2230 pour lequel la luminosite 7 est certainement surestimee comme 
on I'a deja evoque), Ly est inferieur a la valeur attendue selon V£. En particulier, la luminosite 7 
de J0437— 4715 n'est pas conciliable avec la loi empirique, a moins d'un facteur de correction /fj 
irrealiste. On en conclut qu'il doit exister une valeur de E seuil, en de9a de laquelle la relation entre 
Ly et E s'inflechit, afin d'expliquer les efficacites 7 de ces MSPs. D' autres detections de pulsars, 
milliseconde ou normaux, dans ce regime de faible E pourraient aider a comprendre la relation entre 
luminosite et perte d' energie par freinage electromagnetique. 

Dans cette etude nous avons recherche des pulsations a partir d'ephemerides contemporaines pour 
tous les MSPs, a six exceptions pres. Par consequent, il est possible de distinguer les objets pour les- 
quels une emission pulsee a ete detectee, parmi la population de MSPs connus ou plus generalement 
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Figure V.12 - Ly en fonction de E pour les huit pulsars milliseconde detectes (points et traits pleins) et 
les pulsars normaux (carres vides et traits pointilles). Les valeurs de Ly pour les MSPs sont donnees dans 
la Table 



V.l 



Elles sont calculees en utilisant /n = 1- Les barres d'erreur sur Ly pour les MSPs ne tiennent 
pas compte de la contribution de 1' incertitude sur /n. Pour les pulsars normaux, les valeurs de fa et les 



incertitudes associees proviennent du modele Outer Gap (Table 1 de Walters et al. ( 2009 1 et K. Watters, 
communication privee), a savoir/n = 1 sauf pour CTAI (0,6 ± 0,3), Geminga (0,125 ± 0,025), J0205+6449 
(0,95 ± 0,05), J1028-5819 (1,1), J1057-5226 (0,55), J1709-4429 (0,85 ± 0,15), J1952+3252 (0,925 ± 
0,175) etJ2021+3651 (1,05). 



parmi la population des pulsars. La Figure V.l 3 montre les 72 MSPs galactiques connus dans un dia- 
gramme E - distance. On constate que les pulsars pour lesquels des pulsations sont detectees sont 
energetiques et generalement proches. J02 18+4232 est lui plus eloigne, mais a un £ plus fort. PSR 
J0034— 0534, dont les perspectives de detection sont importantes, est egalement un pulsar a grand E 
et a distance faible. En resume, les MSPs detectes se distinguent par un grand E/d^. 

montre E/d^ en fonction de P pour la population de pulsars milliseconde ainsi 



La Figure 



V.14 



que pour les pulsars normaux du champ galactique. Conformement a ce qui est observe dans la Fi- 
gure V. 1 3 les MSPs detectes dominent en terme deE/d^,de meme que les pulsars normaux emetteurs 
7 sont distribues dans la partie haute. On remarque que tous les pulsars detectes sont au-dessus d'un 
seuil situe a environ 5 x 10^-^ erg/s/kpc^, un autre point commun entre les deux populations de pulsars. 
Un certain nombre d'objets non detectes sont neanmoins situes au-dessus de cette valeur {cf. Table 
V.3[ ). Differentes causes peuvent I'expliquer : d'abord les distances peuvent etre erronees, notamment 



lorsqu'elles sont calculees a partir d'un modele de distribution galactique d'electrons, ce qui est le 
cas pour PSR J 101 2-1-5307, J1843-1113, J1911-1114, J1933-6211 et J2129-5721. Rappelons ce- 
pendant que les parametres de rotation utilises pour PSR J1933— 6211 sont anciens, aussi le resultat 
de la recherche de pulsation pour cet objet est a considerer avec precaution. Neanmoins ce pulsar 
n'est pas detecte en tant que source d' emission continue, contrairement aux huit MSPs pour lesquels 
des pulsations sont observees, ce qui signifie que le pulsar n'est vraisemblablement pas detectable a 
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Figure V.13 - E en fonction de la distance pour les 72 MSPs du champ galactique connus. Les huit pulsars 
detectes comme sources pulsees sont indiquees par des points pleins. Les triangles indiquent les pulsars 
associes a des sources / continues. Les carres vides sont les MSPs pour lesquels des ephemerides contem- 
poraines n'etaient pas disponibles. Les autres pulsars sont indiques par des cercles. 



present. En revanche, on dispose pour PSR J1024— 0719 et J1909— 3744 de distances fiables basees 
sur la mesure de leur parallaxe chronometrique ( Hotan et a/. [2006 ). Une autre raison possible de ces 
non-detections est I'orientation de remission : le faisceau 7, s'il existe, peut etre oriente de maniere 
defavorable par rapport a la ligne de visee, empechant ainsi la detection. Enfin, nous avons vu dans 
le premier chapitre que I'incertitude sur le moment d'inertie / des etoiles a neutrons est importante : 
celui-ci pourrait done prendre une valeur plus faible que les 10"^^ g cm^ generalement employes ; de 
sorte qu'a un ralentissement P donne corresponde un taux de perte d'energie E plus faible. 



Pulsar 


P 


E 


Distance 


Reference 




(ms) 


(10" erg/s) 


(kpc) 




J1012+5307 


5,256 


3,7 


0,41 


NE2001 I'Cordes & Lazio'2002| 
Hotan ^f«/. (20061 


J1024-0719 


5,162 


5,318 


0,52 


J1843-1113 


1,846 


60,2 


1,69 


NE2001 (Cordes & Lazio 2002) 
Hotan et al. ( 2006 1 


J1909-3744 


2.947 


4,1 


1,14 


J1911-1114 


3,626 


7 


1,22 


NE2001 (Cordes & Lazio 2002 1 


J1933-6211 


3,543 


3,28 


0,52 


NE2001 (Cordes & Lazio 2002 1 


J2129-5721 


3,726 


22,1 


1,36 


NE2001 (Cordes <fe Lazio|2002| 



Tableau V.3 - Table des MSPs non-detectes en tant que sources d' emission pulsee ou continue, et dont le 
flux E/c/^ est superieur a 5 x 10^^ erg/s/kpc^. 
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Figure V.14 - E divise par la distance au carre en fonction de la periode de rotation pour les pulsars situes 
en dehors des amas globulaires. Les MSPs et pulsars normaux emetteurs 7 pulses sont indiques par des 
points pleins (pulsars normaux : detections de COMPTEL, EGRET et AGILE, et detections recentes du 
LAT ( Abdo et a/.|2009cp 1). Les MSPs probablement associes a des sources 7 continues sont montres par 
des triangles. Les MSPs pour lesquels des ephemerides contemporaines n'etaient pas disponibles pour cette 
etude sont representes par des carres vides. Les pulsars non detectes jusqu'a present sont indiques par des 
cercles pour les MSPs, et par des points fins pour les pulsars normaux. Figure issue de ( |Abdo et al.\2009a\ . 



V.2.2 Conclusion 

Les similitudes des courbes de lumiere, des proprietes spectrales et du role de E dans remission 
7 entre les MSPs et les pulsars normaux detectes suggerent fortement que les memes mecanismes 
de production de rayonnement 7 operent chez ces deux classes d'objets, et ce malgre les ecarts de 
periode de rotation et de taux de ralentissement. Les champs magnetiques surfaciques Bs different 
de plusieurs ordres de grandeur entre les MSPs et les pulsars ordinaires. En revanche, les valeurs du 
champ magnetique au cylindre de lumiere B^c pour les pulsars emetteurs 7, milliseconde ou ordi- 
naires, sont tout a fait comparables. De plus ces pulsars sont generalement distribues aux grandes 
valeurs de B^c^ comme le montre la Figure [VT5 Ceci suggere que les proprietes electromagnetiques 



de la magnetosphere exteme ont un role important dans le processus de rayonnement 7. 

Pour les modeles Polar Cap, le rayonnement prend son origine au-dessus des calottes polaires. 
Par consequent ces modeles predisent des faisceaux de rayons 7 en alignement proche avec remis- 
sion radio. Au contraire, les modeles Slot Gap et Outer Gap prevoient une production de photons 7 en 
altitude dans la magnetosphere, formant des faisceaux larges generalement non alignes avec remis- 



sion radio. Les phasogrammes presentes dans les Figures V.5 et V.6 montrent que les pics 7 sont 
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Figure V.15 - Champ magnetique au cylindre de lumiere Bic en fonction de la periode de rotation pour les 



pulsars situes en dehors des amas globulaires (cf. Figure V. 14 pour la signification des differents symboles) 



toujours en decalage par rapport a la radio, ce qui suggere a nouveau une origine du rayonnement 
dans la magneto sphere externe. Enfin, les energies de coupure mesurees sont de quelques GeV, ce qui 
indique que 1' acceleration de particules ne se produit pas au-dessus des poles ou le champ electrique 
est plus intense, et oil les energies de coupure atteindraient 10 GeV ou plus ( Bulik et a/.|2000[ ). 

En resume, les modeles theoriques pour lesquels 1' acceleration de particules se produit dans des 
regions situees en altitude, comme Slot Gap ou Outer Gap, sont favorises pour les MSPs comme ils 
le sont pour les pulsars normaux ( Abdo et a/.|2009h|e|c|g ). Cependant, les observations ne permettent 
pas a present de discriminer parmi ces deux types de modele. Plus de detections de MSPs, plus de 
details dans les courbes de lumiere et davantage de precision sur les proprietes spectrales pourraient 
permettre des comparaisons particulieres entre modeles theoriques. 
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Chapitre VI 
Pulsars milliseconde des amas globulaires 



Sommaire 
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VI.2 Recherche de pulsations 11061 

VI.2.1 Analyse des donnees 11061 
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VI.3 Conclusion [I09l 



Resume 

PLUS de la moitie des MSPs recenses appartiennent a des amas globulaires, tels que 47 Tucanae 
ou Terzan 5. Nous avons vu precedemment que certains MSPs emettent un rayonnement 
7 detectable par le LAT. De plus, des amas globulaires sont detectes en rayons y en tant que 
sources d'emission continue. Les perspectives de detection de MSPs au sein de ces amas sont done 
importantes, motivant la recherche de pulsations pour ces pulsars. Dans ce dernier chapitre nous 
presentons les resultats de cette recherche de pulsations. Aucun pulsar n'est fermement detecte en 
tant que source de rayonnement 7 pulse, bien que pour la plupart d'entre eux, ces pulsars n'ont pas 
ete couverts lors de la campagne de chronometrie. Un autre facteur redhibitoire est certainement leur 
grande distance, par comparaison aux huit MSPs galactiques detectes. Nous nous penchons enfin sur 
le cas de J1824— 2452A dans I'amas M28, marginalement detecte par AGILE et de prime abord non 
detecte dans notre etude. Bien que la courbe de lumiere vue par le LAT ne corresponde pas a une 
detection ferme, elle est qualitativement dilferente de celle observee par AGILE. 



VI.l Interet 



Les amas globulaires sont des groupes d'etoiles tres denses, situes dans le halo galactique. Ce sont 
des formations anciennes (typiquement, quelques milliards d'annees), interessantes du point de vue de 
revolution stellaire. En raison de la grande densite des amas globulaires, les collisions entre etoiles 



y sont plus frequentes que dans le champ galactique (Verbunt & Hut 1987). Une consequence des 
interactions frequentes entre etoiles est le grand nombre de systemes multiples CClark 1975| ), ce qui 
est propice au recyclage des etoiles a neutrons en pulsars milliseconde. Ainsi, les amas globulaires 
sont riches en MSPs. A ce jour, pres de 150 pulsars ont ete detectes dans 26 amas globulaires {cf. 
Annexe |B} et voir |Camilo & Rasiol ( |2005| ) pour une revue recente). 

Avant I'arrivee du LAT, les amas globulaires avaient ete observes et detectes dans toutes les lon- 
gueurs d'onde, a I'exception du domaine 7. Depuis, le LAT a detecte I'amas 47 Tucanae en tant que 
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Amas globulaire 


TS 


NmSPs 


47 Tucanae 


395,2 


23 


Terzan 5 


310,6 


33 


M62 


112,6 


6 


M28 


106,5 


12 



Tableau VI. 1 - Table des amas globulaires coincidant avec une source continue detectee par le LAT avec 
une significativite superieure a 5 C7, et comprenant au moins un pulsar milliseconde. 



source d'emission continue ( |Abdo et a/. 1200%] ) . Le spectre d'emission de 47 Tucanae est bien repre- 
sente par une loi de puissance a coupure exponentielle, avec un indice spectral F = 1,3 ± 0,3 et une 
energie de coupure Ec =2,5^q g GeV. Par ailleurs, le catalogue de sources continues de Fermi apres 
six mois d'observation comprend trois autres amas globulaires dans lesquels des MSPs sont connus 
(cf. Table VL1[ ). Pour expliquer remission yde 47 Tucanae, un scenario propose est la production de 
vents de particules relativistes par les MSPs, interagissant avec le vent des autres etoiles ou d' autres 
MSPs de I'amas, creant ainsi des chocs capables d'accelerer des particules chargees jusqu'au TeV 
( Bednarek & Sitarekl|2007 ). Les auteurs predisent un flux en energie comparable a celui observe par 
le LAT, mais en se basant sur une valeur de E moyenne pour les MSPs de 47 Tucanae d'un facteur 30 
trop importante ( |Abdo et a/.|2009b[ ). Ce scenario semble a priori ne pas pouvoir expliquer remission 
7 observee. 

Comme le montre la Table V.l[ les proprietes spectrales des huit MSPs galactiques detectes par 
le LAT sont tout a fait comparables a ce qui est observe pour 47 Tucanae. Compte tenu du fait que 
les pulsars milliseconde peuvent emettre du rayonnement 7 pulse, comme nous I'avons vu dans le 
chapitre precedent, remission continue detectee par le LAT dans la direction de I'amas globulaire 
pourrait etre due a la contribution ajoutee des pulsars milliseconde de I'amas, dont 23 sont connus 
mais dont le nombre total pourrait s'elever a 60 objets ( Camilo & Rasio|2005 1. Pour cette raison, nous 
avons egalement recherche des pulsations pour les MSPs des amas globulaires. 



VI.2 Recherche de pulsations 
VI.2.1 Analyse des donnees 

Le lot de donnees considere ici est similaire a celui du chapitre precedent (donnees du LAT en- 
registrees entre le 30 juin 2008 et le 2 juin 2009). De meme, nous avons retenu les evenements de 
classe Diffuse, d'energie superieure a 100 MeV, et compris dans un rayon de 0,5° autour du pulsar si 
\b\ < 10°, ou dans un rayon de 1° sinon. 

Contrairement aux MSPs du champ galactique, I'essentiel des pulsars des amas globulaires n'a 



pas ete suivi par la campagne de datation (cf. chapitre II). La Table VI.2 indique les pulsars pour 
lesquels des ephemerides contemporaines ou recentes ont ete obtenues. 

Pour les autres objets, les ephemerides utilisees pour la recherche de pulsations sont issues de 
la base de donnees de I'ATNF, a I'exception des parametres de datation des pulsars de 47 Tucanae, 



mises a jour des resultats de Freire et al. (2003) mais non contemporains. En tout, des ephemerides 
ont ete coUectees pour 106 pulsars. Pour la majorite d'entre eux, les ephemerides sont anciennes ou 
minimales. Par consequent les resultats de la recherche de pulsations sont a considerer avec prudence. 
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Tableau VI.2 - Pulsars des amas globulaires pour lesquels des ephemerides recentes ou contemporaines 
aux observations du LAT ont ete obtenues. 



La Figure VI.l montre rhistogramme des significativites des courbes de lumiere obtenues pour 
les pulsars des amas globulaires, d'apres le H-test. Aucun MSP n'est detecte a present avec une 
significativite de 3 C7 au moins. La grande distance de ces objets est certainement une cause importante 
des non-detections : les MSPs etudies ici ont des distances comprises entre 4 et 10 kpc, tandis que 
les MSPs detectes par le LAT sont proches, avec des distances inferieures a 1 kpc pour la majorite 
d' entre eux. II faut cependant considerer separement les cas de PSR J02 18+4232 et des sept autres 
pulsars milliseconde emetteurs 7, nettement en de9a en terme de E. Un MSP de luminosite 7 egale 
a celle de J0218+4232 et place a J > 4 kpc produirait un flux en energie h = Ljo2i8/ {^^fo.d^)^ 
soit h < 1,6 X 10^^^ erg/cm^/s, en prenant = L Cette limite superieure est tout a fait comparable 



aux flux en energie vus par le LAT, pour les huit MSPs emetteurs 7 galactiques (voir Table V.ip . Un 



pulsar tel que PSR J02 18+4232 semble done detectable, a faible distance (J ~ 4 kpc). Neanmoins le 
flux attendu decroit rapidement vers des niveaux pour le moment non detectables par le LAT au-dela 
de 4 kpc. Les sept autres MSPs galactiques emetteurs 7 sont moins lumineux et plus proches. Leur 
luminosite 7 etant de Ly ~ 1,2 x 10-^^ erg/s en moyenne, le flux en energie attendu au-dela de 4 kpc 
est /z < 6 X 10^^^ erg/cm^/s, ce qui tres inferieur aux valeurs mesurees actuellement. 

D' autre part, les ephemerides utilisees ici sont generalement approximatives, et empechent cer- 
tainement le calcul precis des phases pour certains pulsars. Neanmoins, remarquons que les pulsars 



correctement couverts, dont les noms sont donnes en Table VI.2 ne sont pas detectes a present. En 
particulier, le cas de PSR J1824— 2452A est interessant, dans la mesure oii une detection marginale 
de ce MSP a ete rapportee par le telescope AGILE au-dessus de 100 MeV. 

VI.2.2 Le cas de PSR J1824-2452A dans M28 

Le telescope 7 AGILE a rapporte la detection du pulsar milliseconde J1824— 2452A, dans I'amas 



globulaire M28, avec un niveau de confiance de 4,2 o (Pellizzoni et al. 2009a). Cependant, cette 
valeur de significativite n'est atteinte que pour une partie seulement des observations d'AGILE, alors 
que le lot de donnees complet ne mene a aucune detection. Pour I'intervalle de temps favorable, la 
courbe de lumiere vue par AGILE est semblable a une sinusoide. Le pic principal radio a 1,4 GHz 
coincide avec le pic 7. 
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Figure VI. 1 - Resultats de la recherche de pulsations pour 106 pulsars d'amas globulaires. Les noms des 
pulsars pour lesquels la significativite est d'au moins 2 a sont indiques. 



La recherche de pulsations dans les donnees du LAT avec le jeu de coupures decrit dans le para- 
graphe precedent conduit, pour J1824— 2452A, a une significativite inferieure a 1 a. Bien que I'article 
d'AGILE ne detaille pas les coupures angulaires utilisees pour ce pulsar en particulier, il est indique 
que la region d' integration utilisee respecte la resolution angulaire de 1' instrument, avec un rayon 
d'ouverture maximal de 2°, pour les regions bruitees. Nous avons utilise ce schema pour produire la 
courbe de lumiere de J1824-2452A, a savoir : p = min(68% PSF, 2°) (rappelons que la PSF du LAT 
est donnee en Figure IV.4). Le phasogramme est montre en Figure [VL2 La significativite correspon- 
dante est de 3,1 o, si Ton ne tient pas compte de I'essai supplementaire. 



Au-dessus de 100 MeV, le profil est en disaccord avec les resultats d'AGILE (jPellizzoni et al. 



2009a) : on n'observe pas de modulation sinusoidale, et le pic principal radio a 1,4 GHz n'est aligne 



avec aucune structure en y. En revanche, il est interessant de constater que les pics radio secondaire 
et tertiaire (ici a 0,75 et 0,25 en phase, respectivement) sont alignes avec des structures emergeant du 
niveau de fond dans la courbe de lumiere 7 au-dessus de 100 MeV. Au contraire, le pic radio principal 
semble ne pas etre accompagne d'une emission 7. Or, ce pulsar milliseconde est emetteur de GRPs, 
principalement en phase avec le pic radio secondaire, les GRPs restants s'alignant avec le pic tertiaire 
( Knight et a/.||2006a| ). Ce parallele entre emission de GRPs et rayonnement 7 est interessant, mais il 
faut neanmoins attendre davantage de photons pour confirmer la detection de J 1824— 2452 A. 
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Figure VI.2 - Courbe de lumiere ypour le pulsar milliseconde J1824— 2452A, dans I'amas globulaire M28. 
Un profil radio de reference, enregistxe au radiotelescope de Nan9ay a 1,4 GHz est egalement montre. 



VI.3 Conclusion 

Parmi les pulsars des amas globulaires, aucun n'est done detecte a present. Outre le fait que les 
ephemerides utilisees dans cette etude n'etaient pas contemporaines des observations ypour la plupart 
d'entre elles, la non-detection de pulsars s'explique vraisemblablement par leur eloignement. De plus, 
rappelons qu'a cause de la resolution angulaire du LAT, qui est de quelques degres a basse energie, 
plusieurs sources ysont potentiellement presentes dans le champ de vue, de sorte que chacune d'entre 
elles participe au bruit de fond des autres objets. Dans I'hypothese oia un pulsar est detecte dans un 
amas, une strategic a envisager pour la recherche des autres pulsars est done d'isoler les photons 
appartenant a ses pics, puis de les eliminer. Ceci diminuerait le niveau de fond pour les autres sources 
presentes dans le champ de vue, ce qui pourrait favoriser la detection de pulsations. 

La detection de 47 Tucanae par le LAT laisse neanmoins a penser que celui-ci pourrait detecter 
des pulsations pour des MSPs des amas globulaires au cours des annees a venir. II serait done interes- 
sant de disposer d' observations radio de ces pulsars afin de construire des ephemerides couvrant les 
prochains mois d'observation y. Bien que la plupart des pulsars milliseconde soient stables, certains 
presentent des orbites complexes, et leur chronometrie requiert un soin particulier. La detection d'au 
moins un MSP emetteur y dans un amas globulaire permettrait vraisemblablement de lever I'ambi- 
guite sur le mecanisme d'emission 7 dans ces formations stellaires. 

Enfin, dans I'hypothese ou des pulsars des amas globulaires soient detectes, I'etude du spectre 
d'emission de ces objets permettrait la comparaison de leurs proprietes spectrales avec celles des 
pulsars galactiques. La mesure des luminosites yfournirait des contraintes sur les valeurs intrinseques 
de E pour ces pulsars, en supposant une efficacite d'emission typique des MSPs connus par ailleurs 
(t] ~ 10%). A cause des mouvements locaux, les taux de perte d'energie rotationnelle intrinseques 
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sont generalement inaccessibles par chronometrie. Bien que les incertitudes soient importantes, en 
particulier sur rj ou sur la distance, des estimations de E apporteraient des informations sur le budget 
energetique des pulsars des amas globulaires, dont le scenario d' evolution et I'environnement sont 
tres differents de ceux des pulsars galactiques. 
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Le satellite Fermi est en orbite autour de la Terre depuis le 1 1 juin 2008. A son bord, le Large 
Area Telescope (LAT) observe depuis le ciel yentre 20 MeV et plus de 300 GeV, avec une sensibilite 
bien superieure a celle de son predecesseur EGRET, ou d' AGILE, autre telescope / operant a I'heure 
actuelle. Avant les observations du ciel par le LAT, seuls neuf pulsars avaient ete detectes dans le 
domaine des rayons 7, grace notamment au telescope EGRET. De fagon remarquable, le LAT a aug- 
mente le nombre de pulsars yconnus d'un facteur cinq : des pulsations ont en effet ete observees pour 
une cinquantaine de pulsars a I'heure actuelle. Ce nombre continuera de croitre, la mission Fermi 
n'achevant que sa premiere annee d'activite. A I'origine de ces decouvertes, le gain en sensibilite du 
LAT par rapport aux instruments anterieurs a naturellement occupe une place importante. Pour ce qui 
est des pulsars, deux autres facteurs ont joue des roles cles : 

- La precision temporelle, testee avant le lancement de Fermi. Ces tests ont justement mis au 
jour des problemes importants dans la datation des evenements, si bien que moins de pulsars 
7 auraient vraisemblablement ete detectes sans ces corrections. En particulier, on pent penser 
qu'aucun pulsar milliseconde n'aurait ete decouvert sans cela. On salt desormais que les pho- 
tons enregistres par le LAT sont dates avec une precision inferieure a la jUs, ce qui conduira a 
un niveau de details sans precedent dans les courbes de lumiere, y compris chez les pulsars les 
plus rapides. 

- la campagne de chronometrie des pulsars emetteurs radio et/ou X, dont le succes a eu diffe- 
rentes consequences. D'une part la collaboration avant lancement a permis, via I'utilisation de 
donnees reelles, la verification et la validation des outils d' analyse temporelle de Fermi. D' autre 
part, le bon fonctionnement de la coordination des efforts de chronometrie pour le LAT a per- 
mis la couverture de presque tons les pulsars definis comme prioritaires (ceux aE > 10^'^ erg/s). 
L'enthousiasme cree autour des resultats de Fermi permet desormais de chercher des pulsations 
pour de nombreux pulsars initialement non prioritaires, done a caracteristiques differentes, ce 
qui pourrait conduire a des resultats inattendus. Cette these a directement beneficie du bon de- 
roulement de cette campagne de datation, puisque I'essentiel des 72 MSPs galactiques connus 
a ete suivi, de sorte que cette population d'objets a quasiment ete etudiee dans son ensemble. 

La recherche de pulsations en rayons 7 pour les pulsars milliseconde, presentee dans cette these, a 
conduit a huit detections fermes. Ainsi, le debat de remission des MSPs dans le domaine 7 est clos : 
les pulsars milliseconde sont de puissants accelerateurs de particules, creant une emission 7 detectable 
a grande distance. L' analyse temporelle, basee sur des ephemerides tres precises construites a partir 
d' observations des radiotelescopes de Nan9ay, Parkes et Green Bank, a revele que les huit MSPs pre- 
sentent un ou deux pics 7, decales de remission radio. Alors que les MSPs sont differents des pulsars 
normaux a priori, de par leur cylindre de lumiere beaucoup plus etroit et leur champ magnetique plus 
faible, les phasogrammes sont tout a fait similaires a ceux des pulsars ordinaires. Les spectres d' emis- 
sion des MSPs ont des indices relativement durs (F < 2) et des energies de coupure de I'ordre du GeV. 
A nouveau, ces proprietes different peu de ce qu'on observe pour les autres pulsars. Enfin, on constate 
que pour les deux classes d'objets, normaux et milliseconde, le LAT detecte les pulsars dont le flux de 
freinage electromagnetique (E / cf-) est important. Le seuil de detection actuel parait commun entre les 
deux populations. Ces differents elements amenent a la conclusion que le mecanisme de production 
de rayonnement 7 est commun entre les pulsars normaux et milliseconde. Les modeles d'emission 7 
de type Polar Cap semblent exclus par les observations, au profit des modeles Slot Gap ou Outer Gap. 
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La prochaine etape dans la comprehension de remission 7 des pulsars est done la discrimina- 
tion des modeles situant 1' acceleration des particules chargees dans la magnetosphere exteme. Pour 
la premiere fois nous disposons de donnees concretes sur remission 7 des MSPs, tant du point de 
vue temporel que spectral. Les differentes observables sent maintenant a confronter aux predictions 
theoriques, afin d'affiner ou d'infirmer les modeles. Les prochaines arniees d'observation du ciel 7par 
le LAT ameneront plus de photons, done plus de details dans les courbes de lumiere, et en particulier 
la possibilite d'etudier remission dans des bandes d'energie etroites. Au bout de cinq a dix ans d'ac- 
tivite, le seuil de detection devrait diminuer d'un facteur deux a trois, augmentant ainsi le nombre de 
pulsars observables en 7et done la variete des caracteristiques. En particulier, le LAT pourrait detecter 
des MSPs appartenant a des amas globulaires, ce qui est pour le moment premature. Pour les MSPs les 
plus brillants, comme PSR J0030+0451, une statistique accrue permettra I'analyse spectrale resolue 
en phase, donnee importante pour la comparaison des modeles theoriques. Par ailleurs, on ne dispose 
d' observations en rayons X avec datation absolue par rapport a remission radio que pour un petit 
nombre d'objets parmi les huit MSPs detectes dans cette these. Ces huit detections motivent I'obser- 
vation dans le domaine des rayons X, afin de cartographier remission au sein de leur magnetosphere 
dans les differentes gammes d'energie. Une fois ameliores, les modeles theoriques vont permettre des 
syntheses de population de MSPs plus realistes, et ainsi foumir des predictions plus precises au sujet 
des MSPs de la Galaxie, leur nombre, leurs proprietes et leur contribution a remission 7 diffuse. 

Finalement, la detection de MSPs emetteurs 7 en bon nombre amene a penser que certaines 
sources non identifiees detectees par le LAT, notamment a haute latitude galactique, sont des MSPs 
pour le moment inconnus. Le faisceau d'emission des MSPs dans le domaine radio est plus large 
que celui des pulsars normaux, a cause du cylindre du lumiere, plus etroit chez les MSPs. On attend 
done moins de MSPs muets dans le domaine radio que chez les pulsars normaux. Par consequent la 
recherche de pulsations pour ces sources non identifiees pourrait etre efficace en radio comme en 7. 
Les perspectives de detection de MSPs emetteurs radio ou non par le LAT au cours des annees a venir 
sont done grandes. Le LAT a deja permis une avancee sans precedent dans I'etude de remission des 
pulsars a haute energie. Tout porte a croire que lorsque la mission Fermi viendra a son terme, le LAT 
aura egalement constitue un allie de poids pour les radiotelescopes, dans la comprehension globale 
de ces objets extremes et de leur environnement. 
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Annexe A 
Table des pulsars milliseconde galactiques 



Les valeurs de la derivee de la periode P, de la luminosite de ralentissement E, des champs ma- 
gnetiques surfacique Bs et au cylindre de lumiere Bijc et de I'age caracteristique T donnees dans les 



tableaux, ainsi que les incertitudes associees, tiennent compte de I'effet Shklovskii ( jShklovskii|1970[ ) 
lorsque possible. Les differentes grandeurs donnees dans les tableaux qui suivent sont obtenues a 
partir des formules ci-dessous. 
- Derivee de la periode : 



^ = ^obs-i^Shk (A.l) 



^p^\|^p^+^pl^ (A.2) 



ou Pshk — 2.43 X 10 ^^P D est la correction de I'effet Shklovskii, avec P en s, D en kpc et 
jl en mas/yr. L' incertitude associee est done : 




APshk = i^ShkX W( — 1 +2^^^ +(^] (A.3) 



La composante due a 1' incertitude sur la periode est negligeable, les contributions de la distance 

D et du mouvement propre transverse /i sont typiquement de I'ordre de 10%. 

Les valeurs de P fournies dans les tableaux sont corrigees de I'effet Shklovskii si la distance 

D et le mouvement propre transverse jl sont connus. De plus on interdit P <0, ce qui peut se 

produire si la distance D est surestimee. Si les conditions precedentes ne sont pas realisees, on 

utilise la valeur de P apparente. 

Luminosite de ralentissement : 



4n^I-^ (A.4) 




AE = £x W3 — + -r (A.5) 



On utilise si possible la valeur de P corrigee de I'effet Shklovskii. Sinon, les valeurs de E et de 
son incertitude sont celles du catalogue ATNF ( Manchester et al. \2005 ), qui n 'emploient que la 
valeur de P apparente. 
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- Champ magnetique surfacique : 



, 3/c3 . 
Bs= \ — ^PP 



(A.6) 




Detail du calcul : cf. E. 

Champ magnetique au cylindre de lumiere : 



(A.7) 



5.c = 4;r%/^/'-V2/.i/2 



(A.8) 



Detail du calcul : cf. E. 
Age caracteristique : 




(A.9) 



P 

2P 



(A.10) 



At 



T X 




(A.ll) 



Detail du calcul : cf. E. 
Remarques : 

- Les valeurs de E, Bs et B^c dependent du moment d'inertie de I'etoile a neutrons /, et les 
champs magnetiques Bs et B^c dependent egalement de son rayon R. Les incertitudes sur / et 
R devraient done intervenir dans les equations |A.5[|A.7]et|A.9[ De plus, on a fait I'hypothese 



de rotateurs orthogonaux dans les equations A.6 et A.8 II faudrait en principe tenir compte 



de I'inclinaison magnetique a et son incertitude. Les contributions de I,R et a sont malheu- 
reusement mal connues, aussi nous n'en tenons pas compte dans les tableaux qui suivent. 



Pour les pulsars dont la distance est calculee a partir du modele NE2001 (Cordes & Lazio 



2002 1, une incertitude sur D de 20% est utilisee. 
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Ci 
KJ 


Jlo, lu 




98 6'^5 




n fio-i- n 19* 

U,OUIlI U, IZ 


fl 9'^R('9"l 
u,zjov^z^ 


■^9 nr'9'i 


90 A(^\ 
ZU,f v^l J 


7 QtA^ 


T14S'^-i-lQn9 




IN 




^iO Rl 


5 7Q9 


1 n f\(fi\ 

1U,DV^D^ 


1 1 5-1- n 9*^* 


9 9^1\ 


87^5\ 


9 5^1"\ 


1 17^5\ 
11 /p; 


J itj J — jjj\J 








99 5fi 


7 087 






1 A(A\ 


7i'9^ 

'l.z; 


3,6(1) 




TlfiOO— 3nS3 







344 09 


16 45 


3 598 


8 8ni 


1 fi^n-l- 326* 


7 5fn 


65n^ 


i,0(,z; 


352(4) 


J1603-7202 







316,63 


-14,50 


14,842 


13,9(5) 


1,170± 0,234* 


0,168(7) 


169(7) 


4,6(1) 


12,8(4) 


J1614-2230 







352,54 


20,30 


3,151 


4(2) 


1,290± 0,258* 


5(2) 


12(6) 


1,1(4) 


3(1) 


J16I8-39 


_ 





340,77 


7,89 


11,987 




2,730± 0,546* 










J1629-6902 




N 


320,37 


-13,93 


6,001 


10,0(3) 


0,960± 0,192* 


1,83(5) 


95(3) 


2,48(5) 


104(2) 


Jl 640+2224 







41,05 


38,27 


3,163 


1,7(2) 


1,160± 0,232* 


2,1(3) 


30(4) 


0,73(7) 


21(2) 


J1643-1224 







5,67 


21,22 


4,622 


16(2) 


2,410± 0,482* 


6,5(7) 


45(5) 


2,8(2) 


26(2) 


J 1709+23 13 


- 





44,52 


32,21 


4,631 


2,0(4) 


1,410± 0,282* 


0,8(2) 


37(7) 


1,0(1) 


9(1) 


J171 3+0747 


_ 





28,75 


25,22 


4,570 


8,07(3) 


1 05+°-'"* 

' — 07 


3,34(1) 


89,8(3) 


1,943(5) 


185,1(5) 


J1721-2457 




N 


0,39 


6,75 


3,497 


6(1) 


1,290± 0,258* 


5,4(9) 


9(2) 


1,5(2) 


31(4) 


J1730-2304 




N 


3,14 


6,02 


8,123 


20,21(1) 


0,530± 0,106* 


1,4887(7) 


63,68(3) 


4,1(1) 


69,55(2) 


J1732-5049 







340,03 


-9,45 


5,313 


13,8(9) 


1,410± 0,282* 


3,6(2) 


61(4) 


2,7(1) 


166(8) 


J1738+0333 







27,72 


17,74 


5,850 


23,1(6) 


1,430± 0,286* 


4,6(1) 


40(1) 


3,72(7) 


169(3) 


J1741+1354 







37,89 


21,64 


3,747 




0,90± 0,18* 










J1744— 1134 




N 


14,79 


9,18 


4,075 


6 9f51 


47+!!-'- 

"■^'-0.07 


4 Of31 


93('71 


1 7('91 




J 1 / — UyjZ 








Q SQ 






i,OJUIt U,jDO 


n S9(^9^ 
u,Jz^^zl 


jzm 


1 ^ liA'\ 
1 J, 


1 7 9f'^^ 

1 /,zpi 


T17S1 9857 
J 1 / J i — Zo J / 






65 


-1 12 


1 Ql 5 


1 1 Of>{A\ 


i , ium u,zz 


7 A^ (^'\ 
1 ,^iPl 


J J, H^z ) 


9 1 9S/^S^ 

Z, iZJt^J f 


jjLi,yyo ) 


J1756-2251 







6,50 


0,95 


28,462 


1017,1(2) 


2,480± 0,496* 


1,7416(3) 


4,4337(9) 


54 445r81 


21,469(3) 


J1757-5322 







339,64 


-13,98 


8,870 


26,3(4) 


0,960± 0,192* 


1,49(2) 


53,4(8) 


4,89(5) 


63,7(7) 


J1801-1417 




N 


14,55 


4,16 


3,625 


5,3(1) 


1,520± 0,304* 


4,4(9) 


108(2) 


1,4(2) 


268(4) 


J1802-2124 







8,38 


0,61 


12,648 


72(1) 


2,940± 0,588* 


1,4(2) 


27,8(4) 


9,66(9) 


43,4(4) 


J1804-2717 







3,50 


-2,74 


9,343 


40,89(8) 


0,780± 0,156* 


1,979(4) 


36,2(7) 


6,255(9) 


69,7(1) 


J1841+0130 







33,12 


2,94 


29,773 


817(5) 


3,590± 0,718* 


12,22(7) 


0,577(4) 


157,8(7) 


54,4(2) 


J1843-1113 




N 


22,05 


-3,40 


1,846 


9,59(5) 


1,690± 0,338* 


60,2(3) 


30,5(2) 


1,346(5) 


1947(7) 


J1853+1303 


_ 





44,88 


5,37 


4,092 


8,9(1) 


2,090± 0,418* 


5,1(6) 


73,3(8) 


1,93(2) 


256(2) 


Jl 857+0943 


Bl 851+09 





42,29 


3,06 


5,362 


17 yxs 




4 42('41 


49 2<S\ 


3,08(2) 


182(1) 


11 90^+0^27 




o 


37 32 


-0 99 


2 150 


18,79(2) 


f. "Jfin-I- 1 272* 


74 65C81 


18,13(2) 


2 034('2"l 


1861(1) 


J1905+0400 


I 


N 


38^09 


-1,29 


3,784 


4,86(6) 


1,710± 0,342* 


3,54(4) 


123(2) 


1,37(1) 


230(2) 


J1909-3744 







359,73 


-19,60 


2,947 


2,7(5) 


1,14±0,05 


4,1(8) 


17(3) 


0,9(1) 


32(4) 


J1910+1256 







46,56 


1,79 


4,984 


9,77(7) 


2,330± 0,466* 


3,12(2) 


80,8(6) 


2,23(1) 


164,0(8) 


J1911-1114 




o 


25,14 


-9,58 


3,626 


8(5) 


1,220± 0,244* 


7(4) 


7(4) 


1,7(7) 


3(1) 


J1911+1347 




N 


47,52 


1,81 


4,626 


17,1(2) 


2,070± 0,414* 


6,83(8) 


42,8(5) 


2,85(2) 


262(2) 


J1918-0642 




o 


30,03 


-9,12 


7,646 


24(3) 


1,240± 0,248* 


2,1(3) 


50(6) 


4,3(4) 


88(8) 


J1933-6211 




o 


334,43 


-28,63 


3,543 


3,7(1) 


0,520± 0,104* 


3,28(9) 


152(4) 


1,16(2) 


237(5) 


J1939+2134 


B 1937+21 


N 


57,51 


-0,29 


1,558 


105,1(4) 


o 01 + 17,0 
»PJ_3,30 


1097,6(4) 


2,3484(9) 


4,095(1) 


9847(3) 
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ANNEXE A. TABLE DES PULSARS MLLLISECONDE GALACTIQUES 



JName 


BName 


Binaire ? 


1 


b 


P 


P 


Distance 


E 


T 


Bs 










(°) 


D 


(ms) 


(10-21) 


(kpc) 


(10^3 erg/s) 


(10** yr) 


{Uf G) 


(10^ G) 


J1944+0907 


- 


N 


47,16 


-7,36 


5,185 


6(3) 


1,790± 0,358* 


1,8(9) 


13(6) 


1,8(6) 


12(4) 


J 1955+2908 


B 1953+29 





65,84 


0,44 


6,133 


28,7(2) 


4,640± 0,928* 


4,92(4) 


33,8(3) 


4,25(2) 


167,4(9) 


J 1959+2048 


B 1957+20 





59,20 


-4,70 


1,607 


8(2) 


2,490± 0,498* 


7(2) 


32(7) 


1,1(2) 


25(4) 






N 


29,45 


-23,54 


< oil 
5,223 


4,82(7) 


1,U2U± U,2U4^ 


1 '2 /I /OA 

1,34(2) 


172(2) 


1,61(2) 


102(1) 








DM-, / J 


-D,DZ. 






1 40ft-l- 9QS* 




oo, I jy^} 













60 86 


-13 12 


5 949 




2 Ozb 4* 




86(8) 


2,6(2) 











39 19 


-30 41 


4 509 




1 040-1- 908* 


5 35('5"l 


57,5(5) 


2,39(2) 


238(2) 


T91 94— '^'^ '^8 




N 


10 93 


-45 44 


4 931 


1 9('8'i 


25+"--^ 


4('3'l 


ov+^ 


3M ■! 




J2129-5721 







338,00 


-43,57 


3,726 


29,0(5) 


1,360± 0,272* 


22,1(3) 


20,3(3) 


3,33(4) 


585(6) 


J2145-0750 




o 


47,78 


-42,08 


16,052 


26(2) 




0,25(2) 


98(6) 


6,5(3) 


14,3(6) 


J2229+2643 




o 


87,69 


-26,28 


2,978 


1,46(2) 


1,45±0,29* 


2,18(3) 


323(4) 


0,667(6) 


230(2) 


J2317+1439 







91,36 


-42,36 


3,445 


2,0(3) 


0,830± 0,166* 


1,9(3) 


28(5) 


0,8(1) 


19(2) 


J2322+2057 




N 


96,52 


-37,31 


4,808 


4(2) 


0,80± 0,16* 


1,4(5) 


20(8) 


1,4(4) 


11(3) 



Tableau A.1 - Table des pulsars miUiseconde du champ galactique 



Les deux premieres colonnes donnent les noms J2000 et eventuellement B1950. La troisieme in- 
dique les pulsars appartenant a des systemes binaires. On donne egalement leur position en coordon- 
nees galactiques, leur periode de rotation P et leur taux de ralentissement P, intrinseque si possible, 
apparent sinon. Les colonnes suivantes donnent les valeurs de d, E, T, Bs et Bic- Les chiffres places 
entre parentheses indiquent I'incertitude sur la demiere decimale citee. Les distances calculees avec 
le modele NE2001 sont indiquees avec des asterisques. Les valeurs de P, P ainsi que les quantites 
derivees pour le pulsar J1614— 2230 sont differentes de celles du catalogue ATNF. EUes sont issues 
de I'ephemeride utilisee pour la detection avec le LAT. Un article traitant de la decouverte et de la 
chronometrie de PSR J1614— 2230 est par ailleurs en cours de redaction (M. Roberts, S. Ransom, 
communication privee). 
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Annexe B 

Table des pulsars milliseconde des amas globulaires 



La Table qui suit donne les noms des pulsars milliseconde repertories (non necessairement confir- 
mes) si tues dans des amas globulaires. Les distances citees ont ete obtenues a partir du catalogue de 
Harris ('Harris 1996), regulierement mis a jouij^ Lorsque possible les valeurs de P et des quantites 
derivees tiennent compte de I'effet Shklovskii (pour le detail, cf. Annexe |A]). Neanmoins les valeurs 
citees sont a manipuler avec prudence, car elles ne tiennent pas compte des effets locaux d' accelera- 
tion gravitationnelle, interferant sur la valeur de P observee. 



JName 


BName 


Amas 


Binaire ? 


1 


b 


P 


P 


Distance 


E 


Notes 










D 


D 


(ms) 


(10-21) 


(kpc) 


(10^3 erg/s) 




J0024-7204C 


B0021-72C 


47Tuc (NGC104) 


N 


305,92 


-44,89 


5,757 




4,5 






J0024-7204D 


B0021-72D 


47Tuc (NGC104) 


N 


305,88 


-44,89 


5,358 




4,5 






J0024-7204E 


B0021-72E 


47Tuc (NGC104) 


O 


305,88 


-44,88 


3,536 


96,7(2) 


4,5 


86,3(2) 




J0024-7204F 


B0021-72F 


47Tuc (NGC104) 


N 


305,90 


-44,89 


2,624 


63,7(1) 


4,5 


139,2(3) 




J0024-7204G 


B0021-72G 


47Tuc (NGC104) 


N 


305,89 


-44,89 


4,040 




4,5 




=1: 


J0024-7204H 


B0021-72H 


47Tuc (NGC104) 


O 


305,90 


-44,90 


3,210 




4,5 




* 


J0024-7204I 


B0()2 1-721 


47Tuc (NGC104) 


O 


305,89 


-44,89 


3,485 




4,5 




* 


J0024-7204J 


B0021-72J 


47Tuc (NGC104) 


o 


305,91 


-44,90 


2,101 




4,5 




* 


J0024-7204L 


B0021-72L 


47Tuc (NGC104) 


N 


305,90 


-44,89 


4,346 




4,5 




* 


J0024-7204M 


B0021-72M 


47Tuc (NGC104) 


N 


305,91 


-44,88 


3,677 




4,5 




* 


J0024-7204N 


B0021-72N 


47Tuc (NGC104) 


N 


305,89 


-44,90 


3,054 




4,5 




* 


J0024-7204O 


- 


47Tuc (NGC104) 


O 


305,90 


-44,89 


2,643 


29,5(2) 


4,5 


63,1(5) 




J0024-7204P 


- 


47Tuc (NGC104) 


o 


305,90 


-44,90 


3,643 




4,5 






J0024-7204Q 


- 


47Tuc (NGC104) 


o 


305,88 


-44,90 


4,033 


34,01(2) 


4,5 


20,47(1) 




J0024-7204R 




47Tuc (NGC104) 


o 


305,90 


-44,90 


3,480 




4,5 






J0024-7204S 




47Tuc (NGC104) 


o 


305,90 


-44,89 


2,830 




4,5 






J0024-7204T 




47Tuc (NGC104) 


o 


305,89 


-44,89 


7,588 


293,7(1) 


4,5 


26,54(1) 




J0024-7204U 




47Tuc (NGC104) 


o 


305,89 


-44,91 


4,343 


93,0(3) 


4,5 


44,8(2) 




J0024-7204V 




47Tuc (NGC104) 


N 


305,90 


-44,90 


4,810 




4,5 






J0024-7204W 




47Tuc (NGC104) 


o 


305,90 


-44,90 


2,352 




4,5 




* 


J0024-7204X 




47Tuc (NGC104) 


N 


305,90 


-44,90 


4,771 




4,5 




* 


J0024-7204Y 




47Tuc (NGC104) 


O 


305,90 


-44,90 


2,197 




4,5 






J0024-7204Z 




47Tuc (NGC104) 


N 


305,90 


-44,90 


4,554 




4,5 




Non publie 


J0514-4002A 




NGC1851 


O 


244,51 


-35,04 


4,991 


1,2(1) 


12,1 


0,37(4) 




J1342+2822A 




M3 (NGC5272) 


N 


42,21 


78,71 


2,545 




10,4 






J1342+2822B 




M3 (NGC5272) 


O 


42,22 


78,71 


2,389 


18,58(4) 


10,4 


53,8(1) 




J1342+2822C 




M3 (NGC5272) 


N 


42,21 


78,71 


2,166 




10,4 




Incertain 


J1342+2822D 




M3 (NGC5272) 


O 


42,22 


78,71 


5,443 




10,4 






J1518+02()4A 


B 1516+02 A 


M5 (NGC5904) 


N 


3,87 


46,80 


5,554 


41,22(6) 


7,5 


9,5(1) 




J1518+0204B 


B1516+02B 


M5 (NGC5904) 


O 


3,86 


46,81 


7,947 




7,5 






J1518+0204C 




M5 (NGC5904) 


O 


3,86 


46,80 


2,484 




7,5 




* 


J1518+02()4D 




M5 (NGC5904) 


o 


3,86 


46,80 


2,988 




7,5 




* 


J1518+0204E 




M5 (NGC5904) 


o 


3,86 


46,80 


3,182 




7,5 






J 1623 -2631 


B 1620-26 


M4(NGC6121) 


o 


350,98 


15,96 


11,076 


62(1) 


2,2 


18,1(3) 




J1641+3627A 


B1639+36A 


M13 (NGC6205) 


N 


59,00 


40,91 


10,378 




7,7 






J1641+3627B 


B1639+36B 


M13 (NGC6205) 


o 


59,01 


40,91 


3,528 




7,7 




* 


J1641+3627C 




M13 (NGC6205) 


N 


59,01 


40,91 


3,722 




7,7 




* 


J1641+3627D 




M13 (NGC6205) 


O 


59,01 


40,91 


3,118 




7,7 




* 


J1641+3627E 




M13 (NGC6205) 


o 


59,01 


40,91 


2,487 




7,7 




* 


J1701-3006A 




M62 (NGC6266) 


o 


353,58 


7,32 


5,242 




6,9 




* 


J1701-3006B 




M62 (NGC6266) 


o 


353,57 


7,32 


3,594 




6,9 







cf. http ://physwww.mcmaster.ca/ harris/mwgc.dat 
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ANNEXE B. TABLE DES PULSARS MILLISECONDE DES AMAS GLOBULAIRES 



JName 


BName 


Amas 


Binaire? 


1 

n 


b 

n 


P 
(ms) 


P 

(10-21) 


Distance 
(kpc) 


E 

(10^3 erg/s) 


Notes 


J1701-3006C 






M62 (NGC6266) 


O 


353,S7 


7,32 


7,613 




6,9 




* 


J1701-3006D 






M62 (NGC6266) 


O 


353,S8 


7,33 


3,418 




6,9 




* 


J1701-3006E 






M62 (NGC6266) 


O 


353,58 


7,33 


3,234 




6,9 




* 


J1701-3006F 


_ 




M62 (NGC6266) 





353,58 


7,33 


2,295 




6,9 




* 


J1740-5340 






NGC6397 


O 


338,17 


-11,97 


3,650 


168(7) 


2,3 


136(6) 


* 


J1748-2021B 


_ 




NGC6440 


O 


7,73 


3,80 


16,760 




8,4 




* 


J1748-2021C 






NGC6440 


N 


7,72 


3,80 


6,227 




8,4 




* 


J1748-2021D 


_ 




NGC6440 


O 


7,73 


3,81 


13,496 


586,8(2) 


8,4 


9,424(3) 


* 


J1748-2021E 






NGC6440 


N 


7,73 


3,80 


16,264 


312,4(4) 


8,4 


2,867(4) 


* 


J1748-2021F 






NGC6440 


O 


7,73 


3,80 


3,794 




8,4 




* 


J1748-2446A 


B1744- 


■24A 


TerzanS 





3,84 


1,70 


11,563 




10,3 




* 


J1748-2446aa 






Terzan5 


N 


3,84 


1,69 


5,788 




10,3 




Non publie 


J1748-2446ab 






TerzanS 


N 


3,84 


1,69 


5,120 




10,3 




Non publie 


J1748-2446ac 






Terzan5 


N 


3,84 


1,69 


5,087 




10,3 




Non public 


J1748-2446ad 






TerzanS 





3,84 


1,69 


1,396 




10,3 




* 


J1748-2446ae 






TerzanS 





3,84 


1,69 


3,659 




10,3 




Non publie 


J1748-2446af 






TerzanS 


N 


3,84 


1,69 


3,304 




10,3 




Non publie 


J1748-2446ag 






TerzanS 


N 


3,84 


1,69 


4,448 




10,3 




Non publie 


J1748-2446ah 






TerzanS 


N 


3,84 


1,69 


4,965 




10,3 




Non publie 


J1748-2446C 






TerzanS 


N 


3,84 


1,69 


8,436 




10,3 




* 


J1748-2446D 






TerzanS 


N 


3,84 


1,69 


4,714 




10,3 




* 


J1748-2446E 






TerzanS 





3,84 


1,69 


2,198 




10,3 




* 


J1748-2446F 






TerzanS 


N 


3,84 


1,69 


5,540 




10,3 




* 


J1748-2446G 






TerzanS 


N 


3,84 


1,69 


21,672 




10,3 




* 


J1748-2446H 






TerzanS 


N 


3,84 


1,69 


4,926 




10,3 




* 


J1748-2446I 






TerzanS 


O 


3,84 


1,69 


9,570 




10,3 




* 


J1748-2446K 






TerzanS 


N 


3,84 


1,69 


2,970 




10,3 




* 


J1748-2446L 






TerzanS 


N 


3,84 


1,69 


2,245 




10,3 




* 


J1748-2446M 






TerzanS 


O 


3,84 


1,69 


3,570 




10,3 




* 


J1748-2446N 






TerzanS 





3,84 


1,69 


8,667 




10,3 




* 


J1748-24460 






TerzanS 





3,84 


1,69 


1,677 




10,3 




* 


J1748-2446P 






TerzanS 





3,84 


1,69 


1,729 




10,3 




* 


J1748-2446Q 






TerzanS 





3,84 


1,69 


2,812 




10,3 




* 


J1748-2446R 






TerzanS 


N 


3,84 


1,69 


5,029 




10,3 




* 


J1748-2446S 






TerzanS 


N 


3,84 


1,69 


6,117 




10,3 




* 


J1748-2446T 






TerzanS 


N 


3,84 


1,69 


7,085 




10,3 




* 


J1748-2446U 






TerzanS 





3,84 


1,69 


3,289 




10,3 




* 


J1748-2446V 


_ 




TerzanS 





3,84 


1,69 


2,073 




10,3 




* 


J1748-2446W 






TerzanS 


o 


3,84 


1,69 


4,205 




10,3 




* 


J1748-2446X 


_ 




TerzanS 





3,84 


1,69 


2,999 




10,3 




* 


J1748-2446Y 






TerzanS 





3,84 


1,69 


2,048 




10,3 




* 


J1748-2446Z 






TerzanS 





3,84 


1,69 


2,463 




10,3 




* 


J1750-3703B 






NGC6441 


o 


353,53 


-5,01 


6,075 


19,2(5) 


11,7 


3,38(9) 


* 


J1750-3703C 






NGC6441 


N 


353,53 


-5,01 


26,569 




11,7 




* 


J1750-3703D 


_ 




NGC6441 


N 


353,53 


-5,01 


5,140 


492,8(3) 


11,7 


143,27(9) 


* 


J1801-0857A 






NGC6S17 


N 


19,23 


6,76 


7,176 




10,8 




Non public 


J1801-0857B 


_ 




NGC6S17 


N 


19,23 


6,76 


28,961 




10,8 




Non public 


J1801-0857C 






NGC6S17 


N 


19,23 


6,76 


3,739 




10,8 




Non public 


J1803-3002B 






NGC6S22 


N 


1,02 


-3,93 


4,397 




7,8 




Non public 


J1803-3002C 






NGC6S22 


N 


1,02 


-3,93 


5,840 




7,8 




Non public 


J1803-30A 






NGC6S22 


N 


1,02 


-3,93 


7,101 




7,8 




* 


J1804-0735 


B1802- 


-07 


NGC6S39 





20,79 


6,77 


23,101 


467,1(3) 


8,4 


1,496(1) 


* 


J1807-2459A 


_ 




NGC6S44 


o 


5,84 


-2,20 


3,059 




2,7 




* 


J1807-2459B 






NGC6S44 


N 


5,84 


-2,20 


4,186 




2,7 




* 


J1823-3021A 


B1820- 


-30A 


NGC6624 


N 


2,79 


-7,91 


5,440 


3384,14(3) 


7,9 


829,871(7) 




J1823-3021D 






NGC6624 


N 


2,79 


-7,91 


3,020 




7,9 




Non publie 


J1823-3021E 






NGC6624 


N 


2,79 


-7,91 


4,394 




7,9 




Non publie 


J1823-3()21F 






NGC6624 


N 


2,79 


-7,91 


4,850 




7,9 




Non publie 


J1824-2452 


B1821- 


-24 


M28 (NGC6626) 


N 


7,80 


-5,58 


3,054 


1617,8(5) 


5,6 


2241,6(7) 




J1824-2452B 






M28 (NGC6626) 


N 


7,80 


-5,58 


6,547 




5,6 




Non publie 


J1824-2452C 






M28 (NGC6626) 





7,80 


-5,58 


4,159 




5,6 




Non publie 
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JName 


BName 


Amas 


Binaire ? 


1 


b 


P 


P 


Distance 


E 


Notes 










n 


D 


(ms) 


(10-21) 


(kpc) 


(10" erg/s) 




J1824-2452E 




M28 (NGC6626) 


N 


7,80 


-5,58 


5,420 




5,6 




Non publie 


J1824-2452F 




M28 (NGC6626) 


N 


7,80 


-5,58 


2,451 




5,6 




Non publie 


J1824-2452G 




M28 (NGC6626) 


O 


7,80 


-5,58 


5,909 




5,6 




Non publie 


J1824-2452H 




M28 (NGC6626) 


O 


7,80 


-5,58 


4,629 




5,6 




Non publie 


J1824-2452I 




M28 (NGC6626) 


o 


7,80 


-5,58 


3,932 




5,6 




Non publie 


J1824-2452J 




M28 (NGC6626) 


o 


7,80 


-5,58 


4,039 




5,6 




Non publie 


J1824-2452K 




M28 (NGC6626) 


o 


7,80 


-5,58 


4,461 




5,6 




Non publie 


J1824-2452L 




M28 (NGC6626) 


o 


7,80 


-5,58 


4,100 




5,6 




Non publie 


J1836-2354A 




M22 (NGC6656) 


o 


9,89 


-7,55 


3,354 




3,2 




Non publie 


J1836-2354B 




M22 (NGC6656) 


N 


9,89 


-7,55 


3,232 




3,2 




Non publie 


J1905+0154A 




NGC6749 


o 


36,21 


-2,20 


3,193 




7,9 






J1905+0154B 




NGC6749 


N 


36,21 


-2,20 


4,968 




7,9 




Incertain 


J1910-5959A 




NGC6752 


O 


336,52 


-25,73 


3,266 


2,18(7) 


4 


2,48(7) 


* 


J1910-5959B 




NGC6752 


N 


336,49 


-25,63 


8,358 




4 






J1910-5959C 




NGC6752 


N 


336,46 


-25,66 


5,277 




4 






J191()-5959D 


_ 


NGC6752 


N 


336,49 


-25,63 


9,035 


964,31(6) 


4 


51,612(3) 




J1910-5959E 


_ 


NGC6752 


N 


336,49 


-25,63 


4,572 




4 






J1911+0101A 


B 1908+00 


NGC6760 


O 


36,11 


-3,92 


3,619 




7,4 




* 


J1911+0101B 




NGC6760 


N 


36,11 


-3,93 


5,384 




7,4 




* 


J1953+1846A 


- 


M71 (NGC6838) 


O 


56,74 


-4,56 


4,888 




4 




* 


J2129+1210D 


B2127+11D 


Ml 5 (NGC7078) 


N 


65,01 


-27,31 


4,803 




10,3 






J2129+1210E 


B2127+11E 


Ml 5 (NGC7078) 


N 


65,01 


-27,31 


4,651 


178(7) 


10,3 


70(3) 




J2129+1210F 


B2127+11F 


Ml 5 (NGC7078) 


N 


65,01 


-27,31 


4,027 


32(8) 


10,3 


19(5) 




J2129+1210H 


B2127+11H 


M15 (NGC7078) 


N 


65,01 


-27,31 


6,743 


2(1) 


10,3 


3(2) 




J2140-2310A 




M30 (NGC7099) 


O 


27,18 


-46,84 


11,019 




8 




* 


J2140-2310B 




M30 (NGC7099) 


O 


27,18 


-46,84 


12,986 




8 




* 



Tableau B.l - Table des pulsars milliseconde des amas globulaires 



Les deux premieres colonnes donnent les noms J2000 et eventuellement B1950. La troisieme co- 
lonne donne le nom de I'amas globulaire hote. La quatrieme indique les pulsars appartenant a des 
systemes binaires. Puis les positions en coordonnees galactiques sont donnees, ainsi que P et P. La 
distance de I'amas (issue de Harris] ( fl 996) ) est donnee ensuite, ainsi que la valeur de E. Enfin la co- 
lonne « Notes » apporte des renseignements complementaires ; indiquant par exemple les pulsars non 
publics jusqu'a present. Les chiffres places entre parentheses indiquent I'incertitude sur la derniere 
decimale citee. 
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Annexe C 
Courbes de lumiere des huit MSPs detectes 



V.6| ) correspondent a 



Les courbes de lumiere presentees dans le chapitre V {cf. Figures |V.5| et 
des coupures angulaires independantes de I'energie : p = 0,5° ou 1°. Or la resolution angulaire de 
r instrument est variable en fonction de I'energie : elle est de plusieurs degres a 100 MeV et elle 
decroit lorsque I'energie augmente. Ainsi, la selection d'une region etroite autour d'une source y a 
pour consequence 1' elimination d'une grande partie des photons de basse energie. Ces phasogrammes 
sont done biaises vers les hautes energies. 

Dans cette annexe nous montrons les courbes de lumiere des huit MSPs pour lesquels des pulsa- 
tions ont ete detectees, en appliquant cette fois des coupures angulaires adaptees au profil de la PSF en 
fonction de I'energie. Cette methode presente I'avantage de ne privilegier aucune partie du spectre des 
pulsars, ce qui est utile pour 1' interpretation de remission des pulsars a differentes gammes d'energie. 

Nous avons retenu les photons enregistres depuis la mise en service du LAT, le 30 juin 2008, 
jusqu'au 2 juin 2009. Les evenements appartenant a la classe Diffuse ont ete selectionnes, et ceux 
dont la direction forme un angle de plus de 105° avec le zenith ont ete exclus. 

La coupure angulaire est determinee de la maniere suivante : soit un photon d'energie E, d'in- 
clinaison et d'azimut (0,0) par rapport a I'axe du telescope et de direction reconstruite formant un 
angle p avec la source consideree. Pour cette configuration (E, 0, 0, p), nous calculons la fraction 
de la PSF qui inclut ce photon (la fonction de reponse de I'instrument utilisee ici est P6_V3_dijfuse). 
Enfin nous Axons une limite sur la fraction de PSF desiree : 68% ou 33% pour les pulsars situes dans 
des regions bruitees ou lorsque des sources brillantes sont a proximite, comme pour J02 18+4232 {cf. 
Table C.l I. La PSF de I'instrument a 68% est representee dans la Figure IV.4 Pour une fraction de 
33% et en incidence normale, elle est de 1,8° a 100 MeV, 0,75° a 300 MeV et 0,3° a 1 GeV 

Les courbes de lumiere y au-dessus de 100 MeV, entre 100 et 300 MeV, entre 300 MeV et 1 
GeV et au-dessus de 1 GeV sont montrees en Figures ICTj |C3l |C3l |C7l |C9l [CTTj |CT3] et [CT5 



Les lignes pointillees horizontales indiquent les niveaux de fond, pour les bandes en energie et les 
coupures angulaires considerees. Ces niveaux ont ete calcules en utilisant 16 regions voisines de 



Pulsar 


PSF 


Note 


> 100 MeV 


Niveaux de fond 
1 00 - 300 MeV 300 - 1 000 MeV 


> 1 GeV 


J0030+0451 


68% 




984(19,7) 


784(15,7) 


193 (3,9) 


4(0,1) 


J0218+4232 


33% 


Blazar 3C 66A a 1° 


575 (25) 


420 (16,8) 


125 (5) 


5 (0,2) 


J0437-4715 


68% 




760 (30,4) 


638 (25,5) 


120 (4,8) 


3 (0,1) 


J0613-0200 


68% 




2632 (105,3) 


2044 (81,8) 


570 (22,8) 


26(1) 


J0751+1807 


33% 


BlazarPKS0735+I7a3'' 


322 (12,9) 


274(11,0) 


56 (2,2) 


0(0) 


J1614-2230 


33% 


Region bruitee 


747 (29,9) 


587 (23,5) 


153 (6,1) 


7 (0,3) 


J1744-1134 


33% 


Region bruitee 


1623 (64,9) 


1265 (50,6) 


347(13,9) 


18 (0,7) 


J2124-3358 


68% 




1227 (49,1) 


937 (37,5) 


285 (11,4) 


13 (0,5) 



Tableau C.l - Fractions de PSF retenues pour les Figures [CT] [C3] [Q] [C^ [C^ [CTT] [CT^ Les 
quatre dernieres colonnes indiquent les niveaux de fond attendus dans les differentes gammes d'energie. 
La premiere valeur est le nombre de photons, la valeur entre parentheses est le nombre de photons de fond 
moyenne sur le nombre de subdivisions du phasogramme. 
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dimensions identiques, situees a 1° du pulsar (les points centraux de ces regions forment un cercle 
autour du pulsar). Pour ces photons, on calcule la phase rotationnelle du pulsar, et on rejette les 
photons appartenant a la region pulsee du phasogramme. On obtient ainsi 16 valeurs dont on retient 
la mediane. Contrairement a la valeur moyenne, la mediane est moins sensible aux perturbations des 
valeurs extremales, introduites par exemple lorsque des sources brillantes sont situees a proximite. La 
Table C. 1 donne les valeurs calculees pour le niveau de fond dans chaque bande d'energie. Les Figures 
C:2|[C4l[C3|[O||O0l[CT2|[CT4l i 



et 



C.16 montrent des cartes du ciel yau-dessus de 100 MeV dans 



un rayon de 5° autour des pulsars, et I'energie des photons yen fonction de la phase rotationnelle. 
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Figure C.2 - Carte du ciel autour de J0030+045 1 et energie 
des photons en fonction de la phase rotationnelle 



Figure C.l - Courbes de lumiere pour J0030+0451 (rayon 
dependant de F energie) 



127 



ANNEXE C. COURBES DE LUMIERE DES HUIT MSPS DETECTES 




des photons en fonction de la phase rotationnelle 



Figure C.3 - Courbes de lumiere pour J0218+4232 (rayon 
dependant de I'energie) 
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Figure C.6 - Carte du ciel autour de J0437— 4715 et energie 
des photons en fonction de la phase rotationnelle 



Figure C.5 - Courbes de lumiere pour J0437— 4715 (rayon 
dependant de F energie) 
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des photons en fonction de la phase rotationnelle 



Figure C.7 - Courbes de lumiere pour J0613— 0200 (rayon 
dependant de I'energie) 
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Figure C.IO - Carte du ciel autour de J0751+1807 et ener- 
gie des photons en fonction de la phase rotationnelle 



Figure C.9 - Courbes de lumiere pour J0751+0751 (rayon 
dependant de I'energie) 
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Figure C.12 - Carte du ciel autour de J1614— 2230 et ener- 
gie des photons en fonction de la phase rotationnelle 



Figure C.ll - Courbes de lumiere pour J 16 14— 2230 
(rayon dependant de I'energie) 
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Figure C.14 - Carte du del autour de J 1744— 11 34 et ener- 
gie des photons en fonction de la phase rotationnelle 



Figure C.13 - Courbes de lumiere pour J1744— 1134 
(rayon dependant de I'energie) 
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Pulse phase 

Figure C.15 - Courbes de lumiere pour J2124— 3358 
(rayon dependant de I'energie) 
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Pulse phase 



Figure C.16 - Carte du ciel autour de J2124— 3358 et ener- 
gie des photons en fonction de la phase rotationnelle 
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Resume 

Le satellite Fermi a ete lance le 11 juin 2008, avec a son bord le Large Area Telescope (LAT). 
Le LAT est un telescope sensible au rayonnement 7 de 20 MeV a plus de 300 GeV. Au debut de 
I'activite de Fenni, neuf pulsars jeunes et energetiques etaient connus dans le domaine 7. Le nombre 
de detections de pulsars par le LAT predit avant lancement etait de plusieurs dizaines au moins. Le 
LAT permettait egalement I'etude des pulsars milliseconde (MSPs), jamais detectes avec certitude a 
tres haute energie jusqu'alors. 

Cette these aborde dans un premier temps la campagne de chronometrie des pulsars emetteurs 
radio et/ou X, candidats a la detection par le LAT, en collaboration avec les grands radiotelescopes et 
telescopes X. Cette campagne a permis la recherche de signaux 7 pulses avec une grande sensibiUte. 
En outre, la plupart des MSPs galactiques ont ete suivis dans le cadre de cette campagne, sans biais 
de selection a priori sur cette population d'etoiles. 

Pour la premiere fois, des pulsations ont ete detectees pour huit MSPs galactiques au-dessus de 
100 MeV. Quelques bons candidats a une detection prochaine apparaissent. Une recherche similaire 
a ete conduite pour des MSPs d'amas globulaires, sans succes a present. L' analyse des courbes de 
lumiere et des proprietes spectrales des huit MSPs detectes revele que leur rayonnement 7 est relative- 
ment similaire a celui des pulsars ordinaires, et est vraisemblablement produit dans la magnetosphere 
externe. Cette decouverte suggere que certaines sources non identifiees sont des MSPs, pour 1' instant 
inconnus. 

Mots-cles : Astronomic gamma, Fermi, Large Area Telescope (LAT), chronometrie des pulsars, 
pulsars milliseconde. 

Abstract 

The Fermi observatory was launched on June 11, 2008. It hosts the Large Area Telescope (LAT), 
sensitive to 7-ray photons from 20 MeV to over 300 GeV. When the LAT began its activity, nine 
young and energetic pulsars were known in 7 rays. At least several tens of pulsar detections by the 
LAT were predicted before launch. The LAT also allowed the study of millisecond pulsars (MSPs), 
never firmly detected in 7 rays before Fermi. 

This thesis first presents the pulsar timing campaign for the LAT, in collaboration with large 
radiotelescopes and X-ray telescopes, allowing for high sensitivity pulsed searches. Furthermore, it 
lead to quasi-homogeneous coverage of the galactic MSPs, so that the search for pulsations in LAT 
data for this population of stars was not affected by an a priori bias. 

We present a search for pulsations from these objects in LAT data. For the first time, eight galactic 
MSPs have been detected as sources of pulsed 7-ray emission over 100 MeV. In addition, a couple 
of good candidates for future detection are seen. A similar search for globular cluster MSPs has not 
succeeded so far. Comparison of the phase-aligned 7-ray and radio light curves, as well as the spectral 
shapes, leads to the conclusion that their 7-ray emission is similar to that of normal pulsars, and is 
probably produced in the outer-magnetosphere. This discovery suggests that many unresolved 7-ray 
sources are unknown MSPs. 

Keywords : Gamma-ray astronomy, Fermi, Large Area Telescope (LAT), pulsar timing, millisecond 
pulsars. 



